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SECTION 1 


INTRODUCTION 


THe purp:)3e of .his manual is to present the capabilities of 
the LATDYN (Large Angle Transient DYNaraics) computer program and 
the instructions necessary for its efficient use. Detailed user 
input and output specifications are provided and discussed , and 
example runstreams, including selected output are given. The 
theoretical basis for the computer program may be found in the 
attached reference document at the end of the manual. 

The manual is structured so that it may be used for both 
tutorial and reference purposes. For tutorial purposes, a new 
user is advised to read Sections 1 through ^ and then proceed 
directly to Section 7, to examine'and duplicate on their computer 
the sample problems of that Section. The sample problems are 
arranged such that in the early problems the more commonly used 
features of the program are demonstrated and more specialized 
commands are shown in the later problems. In this way, the user 
can rapidly learn the commands required for solving most basic 
problems and then proceed to more complex ones. As a reference 
manual, a summary of all input commands and their function is 
provided in Table 4.1-1. 

The LATDYN program consists of the following three 
components; 

(1) a pre-processor to verify and interpret input data 

(2) a set of computational routines 

(3) a post-processor to plot and analyze results 

The entire program may be run in one execution or the components 
may be executed separately. The structure of these three 
components is more fully described in Appendix A. , 

TO date, LATDYN has been configured for CDC/NOS and the DEG 
VAX/VMS computers. Operating instructions for executing LATDYN on 
these computers are provided in Appendices B and C. 

All LATDYN coding is done in ANSII-77 FORTRAN. This provides 
a large measure of portability between computers. Nevertheless, 
small differences may still exist between versions for different 
computers. Detailed instructions regarding the Interfaces with 
particular computer operating systems and file structures are 
contained in Appendices B and C of this manual. 
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SECTION 2 - CAPABILITIES 

The LATDYN program is a tool for analyzing the controlled or 
uncontrolled dynamic transient behavior of interconnected 
deformable multi-body systems which can undergo large angular 
motions of each body relative to other bodies. The program can 
accommodate large structural deformation as well as large rigid 
body rotations and is applicable but not limited to the following 
areas : 

(1) large flexible space structures such as lattice or truss 
type structures which are deployed through unfolding of , 
the lattice deformable beam members from a packaged state 

(2) slewing of large space structure components 

(3) mechanisms which can be modeled as a collection of 
springs, masses, viscous dampers and rigid or deformable 
beam members 

(4) robotic manipulations of beam members 

The current program is limited to two-dimensional problems, 
but in many cases, three dimensional problems can be exactly or 
approximately reduced to two dimensions. 

The program is finite element based and utilizes convected 
coordinates each finite element to affect the large angular 
motions involved in the analysis. Because of its finite element 
nature, the program permits general geometry and users having 
familiarity with other more general finite element programs should 
find LATDYN easy to use. The use of finite elements eliminates 
the user concern with open or closed loop topologies common to 
other multi-body analysis procedures and therefore this 
terminology will not be used in the remainder of this document. 
Those users already familiar with finite element theory may 
comprehend the analytical manipulations behind the user 
instructions, but knowledge of finite element theory is not a 
prerequiste to successful program execution. 

One of the complexities encountered in the solution of multi- 
body transient analysis problems are changes in the structural 
stiffness and mass properties which occur during the anaylsis. 

Many of these changes occur due to contact problems, such as 
encountered in docking and berthing of structural systems, and 
lock-up of deployable Joints during structural deployment. To 
enable the user to direct such structural changes LATDYN uses a 
unique flexible control language. It allows the user external 
access to internal portions of the programs coding, thereby 
enhancing the program’s capabilities. 
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The beginning user is not required to learn the full 
capablities of the control language' in order to use the program 
Rather, it is envisioned that the user will learn the control 
language progressi vely . 
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SECTION 3 - CLASSIFICATIONS OF INPUT DATA 


As an aid in discussion, the input is divided into nineteen 
different classifications as follows: 

1 . Title * 

2. Shorthand Definition 

3- Analysis Control * 

4. Material Properties * 

5. Grid Points » 

6. Beam Members * 

7. Scaling 

8. Springs 

9. Viscous Dampers 

10. Lumped Masses and Additional Inertial Coupling 

1 1 . Preloads 

12. Applied Loads 

1 3 • Reactions 

1 4 . Cons trai nts 

15. Initialization 

16. Logical Conditions 

17. Variable Assignments 

18. Lock-Up of Joints 

19. Output Control * 

Classifications marked by an asterisk (*) indicate that at 
least some input data within this classification is always 
required. Each of these classif cations will be discussed in 
detail in the next section. 

No specific order is required for input data, however, a 
suggested order is employed in the sample problems of Section 7 
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SECTION 4 - GENERAL FORMAT OF INPUT DATA 
Section 4.1 - .Command Names 

Input is supplied to the program in the form of commands. 

Each command may consist of many lines of data with each line 
being up to 80 columns in length. Each command has an associated 
alpha-numeric name which the user generally provides in the first 
line of the command. The user selects from the list of admissible 
command names in Table 4.1-1. The function of each of these 
commands is described later;- In addition Table 4.1-1 provides a 
brief format description of each command for quick reference. 

Commands which are placed in square brackets, [ ], are 
not available to the user in this release. Commands which are 
placed in braces, { }, are post-prbcessor commands. This is 
explained at greater length on page 5.19-8 

The command name commences in the first column of the first 
line of a command and is followed by a delimeter which may be 
either a comma or a colon. Following the delimeter, is a sequence 
of alpha-numeric data values placed in a specific order and 
separated by commas. 

In general the first line of a command takes the form, 

command name: sequenced data 
0 r 

command name , sequenced data 
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If a name is absent from the first line of a command, it is 
assumed that the command name is a repeat of the command name 
which proceeded it. In such a case a delimeter is not required 
and only the sequenced data is given, such as, 

command name : sequenced data 
sequenced data 

Blank spaces between data entries are generally ignored, (except 
that imbedded blanks in ah alpha-numeric name or a numerical data 
value are not permitted), so that it is also permissible to write, 

command name : se quenced data 
sequenced data 

which has the same effect as the preceeding example and is useful 
for visual inspection of data input. 

A list of admissable command names is given in Table 4.1-1. 
The list is broken down into the various classifications of 
Section 3- Required commands are flagged. Also Table 4.1-2 
provides an alphabetical listing of commands with the associated 
data format. Section 5 provides full documentation and examples 
for each command. 
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Table H.1-V Command Name Liisting by Classification 

■ ■ i i 

t 

Cl assification i Z£8£_No^ 

Re q ui red TITLE. . .j 5.1-1 

Requi re d PORDER . ..5.3-1 

Required TIMS PAN 5.3-6 


Shorthand D ef i ni t i'on ...... DEFINE '....' 5 

LIST 5 

Analysis Control REST;ART 5, 

INCMASS 5 

STIFF .5 


FREQ 


..STOPj 15 

Time Integration INTGTYP 5 

. .TIMS'TEP 5 

Material Properties.; MAT PR OP * 5 

Grid Points GRID .^ 5 

Beam Members I MEMBER 5 

Scaling SCALE 5 

Springs SPRING 5 

Dampers . . . . ; DAMPER ...... .5 

Lumped Mass 1 . . ADMASS 5 

Additional Inertial 

Cojupling ADMASS .5 

Preloads ... i PRELOAD ..... 5 

Applied Loads... APPLOAD...... .' .5 

; LOAD i . 5 

I : LOADSIN 5 

' LOADTIM 5 

; GRAVITY 5 

Constraints SDFC 5 

i ....... ^ MDFC 5 

i i PIN 5 

! • ;....... INDDOF 5 

In'itialization....; INIT 5 

. . . . i DIS 5 

. . . . i VEL .5 

Condition Statements 

and Labels. CLj or Cj 5 

Variable Assignments SET 5 

RESET 5 

; KEEP 5 

LOCKUP 5 


2-1 

2-2 

3-2 

3-^ 

3-7 

3-8 

3-10 

3-3 

3- 6 

4- 1, 

5- 1 

6 - 1 

7- 1 

8 - 1 

9- 1 

10 - 1 

10-1 
11-1 
1 2-1 

12- 5 
1 2-6 
12-8 
12-10 
1 3-1 

13- 3 
13-6 
1 3-8 
1 4-1 

14- 1 

14- ^1 

15- 1 

16 - 2 
1 6-4 
16-8, 
17-1 


4.1-3 Table of Command Names by ,C 1 ass i f i ca t i on 



Table ^4.1-1 Command Name Listing by Classification (Continued) 


: Com ma n d Name Z£8® 

Printing Output PRINT 5.18-^ 

GRDPRT 5.18-7 

MEMPRT 5.18-8 

REACT 5.18-9 

Plotting Output PLOTINC ...5.18-10 

. PLOTLIMIT 5. 18-10 

PLOTXLOC ..5.18-13 

PLOTYLOC 5.18-13 

PLOTROT ; 5.18-13 

PLOTMAG 5.18-14 

PLOTX 5.18-14 

. PLOTY .5:18-14 

PLOTPHI 5.18-14 

PLOTXP . . : 5.18-15 

PLOTYP 5.18-15 

PLOTANGLE 5.18-15 

PLOTSTRAIN 5.1 8-1 5 

PLOTFEXT 5.18-15 

PLOTC .5.18-16 

PLOTQ 5 . 18-16 

PLOTAX ; 5.18-16 

PLOTSHA 5.18-16 

PLOTSHB 5.18-16 

PLOTMOMA 5.18-16 

PLOTMOMB 5.18-16 

PLOTQVQ 5.18-17 

SNAPSHOT 5.18-18 

PSD .5.18-19 

PSL ....5.18-19 

TIMEWINDOW 5 . 18-20 

FREQWINDOW 5.18-20 

HELP 5.18-21 


4,1-4 Table of Command Names by Classification 
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Table 4.1-2 Command Names and Data Format 

C ommand_N ame Z£.8£_No_i 

ADMASS: grid point number and direction for two 

degrees-of -freedom which have an additional 

mass coupling, followed by value of mass 

coupling ? condition label .5.10-1 

APPLOAD: identification i nteger , angul ar orientation, 


scale factor, time values at which load 

comes on and off, grid points at which 

load is to be applied 5.12-1 

CL j : or Cj defines j th condition label using a 

FORTRAN expression 5.15-1 

DAMPER: grid point numbers at ends of damper, 

integer specifying or i ent at i on , damper 
constant? condition label 5.9-1 

DEFINE: @n = constant (where n is an integer) 5.2-1 

DIS: grid point number, direction, initial 

displacement value 5.14-1 


FREQ: integer indicating how often frequencies are 

ca 1 cul at ed , i nt eger specifying how many 
frequencies are desired, integer specifying 


how many mode shapes are desired 

? condition label 5.3"8 

FREQWINDOW: two frequency values specifying frequency 

window for PSD command 5.1S-20 

GRAVITY: acceleration due to gravity at sea level, 

X and Y coordinates of earth’s center 5.12-8 

GRDPRT: string of grid points at which motions 

are printed 5.18-7 


GRID: starting grid point number of string, number 

of grid points in string, x and y global 
coordinates of first and second grid points 
in s t r i ng . ( r equi r ed command) 5.5~1 

Unless otherwise noted, the use of a command name is optional. 
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Table Command Names and Data Format (Continued) 


Command 

INCMASS 

INIT : 
KEEP : 

INTGTYP 

LOAD: 

LOADSIN 

. LOADTIM 

LOCKUP: 

MATPROP 

MDFC: 

Unless 


Name Data Format Page No. 


time step increment for updating mass 
matrix, angular rotation for updating 
mass matrix ? condition label.... 5.3~1 

scaling value for initial displacements, 

scaling value for initial velocities 5.H-1 

Qn values to be saved, integer number of 

previous time steps at which designated Qn 

values are to be saved I 5.16-8 

integer specifying time integration 

algorithm 5 . 3 "3 

identification integer, load value 5.12-11 

identification integer, sinusoidal 

frequency, phase shift and amplitude 5.12-6 

identification integer, relative or 
absolute time, pairs of time, load values for 
piecewise linear definition of load 5.12-8 

a series of coefficients describing the 

locking joint model ? condition label ...5.18-1 

material name or number, E, A, I , density 
? condition label 5.11-1 


number of degress-of-f reedom linearly 
constrained together, grid point numbers, 
directions and multipliers for each 
degree-of-freedom constrained, right-hand-side 
of linear constraint relationship 

? condition label 5.13-3 

otherwise noted, the use of a command name is optional. 
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Table 4.1-2 Command Names and Data Format (Continued) 


Command Name Data Format Page No. 


MEMBER: grid point number of ends A and B of 

finite element, number of finite elements 
in string, material name or number 5.6-1 

MEMPRT: String of member numbers at which internal 

forces are printed.. ..5.18-8 

PIN: two grid point numbers which have a 

connecting pin between them allowing full 
rotational freedom 5 . 13"6 

PLOTAX: member numbers in which axial loads are 

plotted ...5.18-16 

PLOTC: list of condition labels whose TRUE-FALSE 

values are plotted 5.18-16 

PLOTFEXT: M or A for magnitude or orientation angle 

of plotted external forces, load set 

identification integers 5.18-15 

PLOTINC: time step for saving data for plotting 

? condition label 5-.18-11 

PLOTLIMIT: lower and upper time limits for plotting 5.18-11 

PLOTMOMA: member numbers in which moments at end A 

are plotted 5.18-16 

PLOTMOMB: member numbers in which moments at end B 

are plotted 5.18-16 

PLOTQ: list of user defined Q-variables whose 

values are plotted 5.18-16 

PLOTQVQ list of user defined Q-variabies which are plotted 

aganist the first Q-variable of the list 5.18-17 

Unless otherwise noted, the use of a command name is optional. 
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Table Command Names and Data Format (Continued) 


Command Name Data Format Page No. 


PLOTSHA: member numbers in which shear forces 

at end A are plotted 5.18-16 

PLOTSHB: member numbers in which shear forces 

at end B are plotted 5.18-16 

BORDER: number of finite elements, number of grid 

points (required command) 5.3~1 

PRELOAD: member number, preload value 5.11-1 


PRINT: time step increment for printing grid point 

motions, coordinate system for motions, time 
step increment for printing internal forces, 


coordinate system for forces, time step 
increments for printing mass and stiffness 
matrices ? condition label 5.18-M 

PSD: list of parameters to creat power spectral 

density plots 5.18-19 

PSL: list of parameters to create power spectral 

level calculations 5.18-19 

RESET: Qn = FORTRAN expression (n is an integer) 5.16-11 

RESTART: integer indicating restart 5.3-2 

ROTLOCK: two grid point numbers that rotationally 

lock? condition label 5.17-5 


Unless otherwise noted, the use of a command name is optional. 


Table of Commands & Data Format 
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Table 4.1-2 Command Names and Data Format (Continued) 


C ommand_N arae 
SCALE: 

SDFC: 

SET: 

SNAPSHOT : 
SPRING: 

STIFF : 

STOP : 
SUBTITLE : 
TIMEWINDOW: 

TIMSPAN : 
TIMSTEP: 
TITLE : 

VEL : 


Data Format 


Page No . 


scale for x direction, scale for y 

direction, rotation ....5.7-1 

grid point number, direction ? condition 

label 5.13-1 

Qn = FORTRAN expression (n is an integer) 5.l6"-2 

list of parameters for plotting a snapshot 
of structure 5.18-18 


grid point numbers at ends of spring, 
integer specifying spring orientation, 
spring cons t ant , i ni t i al spring stretch 
? condition label .....5.8-1 

integer specifying how often tangent 
stiffness is to be calculated, Y if 
incremental stiffness is to be Included 


or N if not ? condition label 5.3-7 

? condition label 5.3-10 

data case subtitle (up to 71 characters)..... 5. 2-1 
two time values specifying plot time 

limits 5.18-20 

start time, stop time ..5.3-6 

time step, size ? condition label 5.3-6 

data case title (up to 74 characters) 5.1-1 

grid point number, direction, initial 

velocity value 5.14-1 
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Section 4.2 - Continuation Lines 

Commands of more than one line in length may be constructed 
through the use of continuation lines. A continuation line do es 
not commence with the command name, but continues with the data 
input as though only one line were being used. The line 
preceeding a continuation line must terminate with an & symbol. 

The limit on the number of continuation lines that can be used is 
ten, with the exception that TITLE and SUBTITLE commands do not 
permit any continuation lines. 

In general, a command with continuation lines takes the form, 

command name: sequenced data & 
sequenced data 

Note that an input line without a command name is either a 
continuation line if the last entry of the preceeding line is an & 
symbol or is the first line of a command with the same name as the 
preceeding command if the last entry of the preceeding line is not 
an & symbol. 
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Section M.3 ~ User Supplied Comments 

Users are encouraged to make liberal, use of comments to 
annotate their input data. A sign at any location on a line 

of data indicates that whatever follows on that line is a comment 
and will be ignored by the program. Commented commands take the 
general form, 

command name: sequenced data $comment 

Continuation lines solely for comments require a in the first 

column of the continuation line. 

The user may wish to separate data types with blank lines by 
placing a sign in the first column on the line separating the 

data types as , 

command name : sequenced data $comment 

$ 

command name : sequenced data $ comment 


Similarly, the user may wish to tabulate the data associated with 
a command for ease of visual inspection as. 


$ 

heading A 

heading B 

heading C 

command name: 

data 1 A 

data IB 

data 

1 C 


data 2A 

data 2B 

data 

2C 


data 3A 

data 3B 

data 

3C 


NOTE The user is cautioned to avoid placing a continuation 


& symbol after a $ symbol, since it will be 
interpreted as part of a comment. 


% 
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Section 4 . M 


Conditional Commands 


Many of the admissible commands can be made conditional. 

That is, the user can specify conditions under which the command 
will be active or inactive. This is accomplished through the use 
of a condition label which is placed at the end of the command 
following the last data entry of the command. The condition label 
is separated from, the last data entry by a question mark, "?" such 
as , 


command name: sequenced data ? condition label 

A condition label identifies a user supplied TRUE or FALSE logic 
statement which makes the data preceedihg it conditional; that is, 
the data is used only when the logic statement is TRUE. The logic 
statement is written in FORTRAN V and can make use of variables 
calculated during program execution such as displacement, 
velocity, acceleration, member strain, etc. A complete list of 
variables available to the user is provided in Table 4.4-1. 

The employment of condition labels gives the user the 
flexibility of specifying such things as contact conditions during 
impact phenomena, lock-up conditions for joints in deployable 
structures, control of parameters used in data output and a host 
of other functions. Section 5.15 provides further information on 
the use of condition labels. 
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Table 4.4-1 Admissible Variables For Condition Labels 
and User Defined Variables 

Variable Name Explanation 


XLOC(n) for the current location of grid point n parallel 

to the global x axis 

YLOC(n) for the current location of grid point n parallel 

to the global y axis 

ROT(m) for the current rotated position of member m 

measured positive in a direction from the positive 
global X to the positive global y axis 

X(p,n) for component of motion quantity (displacement, if 

p=Dj velocity, if p=V; acceleration, if p=A) at 

grid point n parallel to the global x axis 

Y(p,n) for component of motion quantity (displacement, if 

p=D; velocity, if p=V; acceleration, if p=A) at 

grid point n parallel to the global y axis 

XP(p,n,l) for component of motion quantity (displacement, if 

p=D; velocity, if p=V; acceleration, if p=A) at 

grid point n parallel to a local x'axis (parallel 
to a designated member, 1) 

YP(p,n,l) for component of motion quantity (displacement, if 

p=D; velocity, if p=V; acceleration, if p=A) at 

grid point n parallel to a local y'axis (transverse 
to designated member, 1) 

where , 

p is a parameter indicating the type of motion desired as, 
D, for displacement 
V , for veloci ty 
A, for acceleration 

1 designates the member number associated with local axes 
m is a member number 
n is a grid point number 



11-1 5-85 


Table of Admissible Variables 


Table 4 

Va riable Name 
PHI (p ,n) 

MAG(p,n) 

ANGLECp ,n) 

PI 

PIH 

PIQ 

PI2 

where , 
p is a 


n i s a 


4-1 Admissible Variables For Condition Labels 
and User Defined Variables (continued) 

Explanation 


for rotation at grid point n in the positive 
angular direction; from positive global x axis to 
positive global y axis 

for magnitude of motion quantity (displacement, if 
p=D; velocity, if p=V; acceleration, if p=A) at 
grid point n 

for angular orientation of motion quantity 
(displacement, if p=D; velocity, if p=V; 
acceleration, if p=A) at grid point, n; positive ; 
angles are measured from positive global x axis 
towards positive global y axis 

for.. IT 
foriT/2 
for it/ 4 
for 2it 


parameter indicating the type of motion desired as, 
D, for displacement 
V , f or vel oc i ty 
A, for acceleration 
grid point number ' 


-15-85 
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Table ^.^*-1 Admissible Variables For Condition Labels 
and User Defined Variables (continued) 

Variable Name Explanation 


AX(m) for axial force in designated member m (positive 

for tension) 

SHA(m) for shear force at end A of beam member m 

SHB(m) for shear force at end B of beam member m 

BM(m,sbar) for bending moment in ra'' beam member at a 

dimensionless distance sbar (0 < sbar < 1) from end 
A of the member 

FEXTX(j) for component of j external force parallel to the 

global X axis 

t h 

FEXTY(j) for component of ,j external force parallel to 

global y axis 

t Yx 

FEXTM(j) for the magnitude of the j external force 

FEXTA(j) for angular orientation of the external force; 

positive angles are measured from the positive 
global X to positive global y axis 

where , 

m is a member number 

th 

j refers to the applied load specified by the j 
APPLOAD command 


sbar designates location at which bending moment is 
desired as the distance from end A of the beam 
ratioed to the beam length 



Table ^.4-1 Admissible Variables For Condition Labels 
and User Defined Variables (continued) 



REACTX(n) for reaction force at grid point n in the global x- 

direction 

REACTY(n) for reaction force at grid point n in the global x- 

direction 


REACTK n) 
REACTM(n) 
REACTA( n) 
REACTC ( n , 1 ) 

STRAIN (m) 

T . 

NTSTEP 
0MEGA( i ) 

CLj or Cj 

QH ( Qk , r ) 

where , 


for reaction torque at grid point n in the global 
x-direction 

for the magnitude of the reaction force at grid 
point n 

for the angular orientation of the .angular force at 
gridpointn 

for the component of the reaction force at grid 
point n parallel to the local coordinate system of 
member 1 

for axial strain in member m 
for time 

for time step number 

for frequency in radians per unit time of 
designated mode i, where the absolute value of i 
designates the i th lowest mode if i is negative 
and the i th highest mode if i is positive 

t 

for the j logical condition; provides logical 
TRUE or FALSE for condition label Cj 

for the. value of user defined variable Qk at a time 
r time steps prior to the current value of time 


m is a member number 
i is a mode number 

r is number of time steps prior to current time 


4.4-5 Table of Admissible Variable; 
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Section 4.5 - Multiple Occurrences of a Command Name 

The user may supply multiples of most commands. For example, 
there may be any number of GRID commands in the input data. 
However, for the following command names, only one command can be 
active at a time: TIMSTEP, INCMASS, STIFF, FREQ, MATPROP, PRINT 
and PLOTINC. To determine which command of a multiple set is 
active, the program uses a condition-order priority as explained 
in NOTE 5.15-1 of Section 5 . 1 5 When multiples of these commands 
are employed, the first of these Commands which appears in the 
data set cannot have an attached condition label, while the 
subsequent ones must each have an attached condition label. The 
preprocessor will search for violations of this rule and if any 
violations are found, execution will terminate. 
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Section 


Units for Input Data 


Any consistent set of units may be used with the exception 
that angular measurements must be in radians and frequencies in 
Hertz. 
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Units 



SECTION 5 - INPUT SPECIFICATIONS 


Sect! on 5. 1-Title 

Titles for data case identification are specified by using 
the commands, TITLE and SUBTITLE. Character strings following 
either of these commands may contain any ASCII standard characters 
in any order, but may not exceed one line in length (80 
characters) including command name. The designated title will 
appear at the beginning of the printed and plotted output. Use of 
^ the SUBTITLE command is 
optiona^. Data takes the form, 

TITLE, free field descriptor 
SUBTITLE, free field descriptor 

TITLE, DEPLOYMENT OF 10 BAY TRUSS ,, 

SUBTITLE, UNIFORM PROPERTIES 
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5.1-1 


TITLE & SUBTITLE 



Section 5 . 2 


Shorthand Definition 


The DEFINE and LIST commands are op t 1 onal and may be used to 
define constants or lists of integers which often appear in the 
input data. Other commands can use these values in a shorthand 
fashion. 

DEFINE 

Shorthand definitions begin with the name DEFINE, This is 
followed by a colon or comma an sign, an Integer, an equal 

sign and the specified constant value. Defined constants may then 
be used in other commands by using the sign followed by the 

appropriate integer in place of a numerical value. Data takes the 
form, 

DEFINE , §n=cons t ant 
where , 

ri is an i nteger 


Example_5^2- Def ^ni ng_co ns t ants ■ 

DEFINE, §1 =3. 1.^1 59265359$. DEFINE PI 

DEFINE, @99= 1 . 41 i|21 36$ DEFINE SQUARE ROOT OF 2 

A condition label ca nnot be us.ed in conjunction with a DEFINE 
command. 

N0TE_5^22i2. The user may define up to 99 constants. 

NOTE 5.2- 2: The @n may appear in place of any numerical value in 

nearly all commands. Forms of @n requiring a sign 
change are used as "~@n". 
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DEFINE 



LIST 


The list command is an option al command used to def i ne ’ a 
sequence of integers which can be used in a shorthand fashion in 
certain commands. At present, only the APPLOAD , GRDPRT, MEMPRT 
and LIST commands are able to use the defined lists created by a 
LIST command. Data takes the form, 

LIST,n,m1,m2,m3,... 

A list command may also use the symbols T and B. The T 
implies "through" and is used to define a string of quantities, 
and the B implies "by" and is used to designate the spacing 
between subsequent items in the list. 

Ex ampl e_5j^222j D ef ining a list of int eg e rs 

LIST, 1 , 1 , 2, i<T1 0 

In this example a list of integers 
1,2,4,5,6,7,8,9,10 is defined 

Exampl e_5^22^3 j Defining a list of equally spac ed 

integers 

LIST,20T30B 

In this example, a list of integers is defined with 
spacing of two, namely, 20,22,24,26,28 and 30. 

NOTE 5.2-^: Floating point values used in a LIST command are 

' truncated to integer values. 

NOTE 5.2-4: A LIST may not specify more then fifty integers. 
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Section 5.3 - Analysis Control 

Several commands are used to specify parameters controlling 
program execution. Some of these commands are required; others 
are optional. The following commands are employed: 

PORDER 

This is a required co mm a nd which defines the maximum 
structural model size as the number of finite elements in the 
model and the number of grid points in the model. Data takes 
the form, 

PORDER, maximum number of elements, maximum number of grid points 

Ex am pl e 5 . 3~ 1 : Specifying the ma ximum problem s iz e # 

PORDER , 50, 1 20$ A MAXIMUM OF 50 FINITE ELEMENTS AND 
$120 GRID POINTS ARE ALLOCATED IN THIS EXAMPLE 
In this example, computer storage requirements for 
the problem are allocated to allow for a maximum of 
50 finite beam elements and 120 grid points. 

A condition label ca nnot be used with this command. 

NOT E 5 . 3 - 1 : The values specified by the user must not be exceeded 

in the geometry of the problem. Since the computer 
will be instructed to set aside memory sufficient to 
handle the problem dimensions, it is efficient, but 
not necessary, to set the maximum problem dimensions 
equal to those actually existing in the structural 
model . 
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PORDER 



RESTART 


This is an optional command which indicates if a restart case 
is being executed or if a restart tape is to be generated. Data 
takes the form, 

RESTART, IN - indicates that a restart case is being 

executed 

RESTART, OUT - indicates that a restart tape is to be 

generated 

RESTART, BOTH - indicates that a restart case is being 
executed and that a restart tape is to be generated. 

A condition label cannot be used in conjunction with a 
RESTART command. If a RESTART command does not appear in the 
input data, it is assumed that neither is a restart case 
being executed nor a subsequent restart envisioned. 
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INTGTYP 


This is an op ti onal co mmand which indicates the time 
integration algorithm to be employed in the analysis . Data 
takes the form, 

INTGTYP, 1 - indicates the use of central differences 


INTGTYP, 2 - indicates the use of the Newmark Beta method 
in an explicit fashion with beta equal to 
unity. 

INTGTYP, 3 ■ indicates Adams-Bashf or th predictor in 

conjunction with Adams-Moulton corrector and 
automatic variable step size 

INTGTYP, 4 - indicates use of, fifth order Runge Kutta 
with automatic variable step size 

INTGTYP, 5, B - indicates the use of the Newmark Beta 
method in an implicit fashion. If the 
value of beta is left blank, a default 
value of 1/3 is used. 


If the INTGTYP command is omitted from the data set, th e 
de fault is the Newmark Beta m ethod i n an explicit fashion with 

j condition labels ca nnot be used 

with this command. 

NOTE 5 . 3~2 : INTGTYP values of 2 through 5 are not available in 

the current version of LATDYN. 
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INCMASS 


This is an optional com mand used to designate updating of the 
system mass matrix. Data takes the form, 

INCMASS, n, angle ? Cj 
where, 

n is a positive integer specifying that the mass matrix 
is to be periodically updated with a periodicity of n 
time steps; if zero or blank periodic updating of the 
mass matrix will not occur 

angle is the maximum magnitude of angular motion, in 
radians, that any rigid member, (see MEMBER 
command in section 5.*<), will be permitted to 
rotate without there having occured an updating of 
the mass matrix; if "angle" is zero or left blank, 
a default value of .0872 radians (5 degrees) is 
assumed 

Cj is the j th opt i o nal condition label defined with a 
CJ command; see Section 5.15 for further details. 

Example 5. 3~2: Specifying mass matrix updates 

INCMASS, 1 00 , 0 . 2$ MASS MATRIX UPDATED EVERY 100 TIME 
$STEPS 

In this example, the mass matrix is updated every 
100 time steps; however, it will be updated sooner 
if any rigid member of the structure rotates plus or 
minus 0.2 radians since the previousmass matrix 
update . 


11-15-85 5.3-^ INCMASS 



NOT E 5 . 3 ~3 •• Updating of the mass matrix will occur every n time 
steps, where the length of the time step is that 
value which is currently being used by the program. 

In addition, updating of the mass matrix Will occur 
when the increment of angular motion of a rigid 
member, since the most recent update, (independent of 
what caused the update), reaches the specified 
"angle" value. 

N OT E 5 . 3 -^ : Multiple INCMASS commands may be used, but they have 

a CO nd i t i on- or der priority as discussed in NOTE b.15- 
1 of Section 5.15. Thus, the first appearance of an 
INCMASS command in the data set cannot have an 
attached condition label. See also Section M.5. 

N0TE_5^325: If the INCMASS command is omitted from the data set, 

the mass matrix will not be updated on a periodic 
basis. However, events whi.ch can lead to large 
changes in the mass matrix such as imp<^ct or lock-up 
of joints are sensed by the program and result in 
automatic updating of the mass matrix. The printed 
output will indicate this occurrence. (If rigid 
elements are not employed in the modeling, periodic 
mass matrix updating is usually not necessary and 
this command may be omitted if there is no other 
reason to update the mass matrix. For further 
guidelines on updating the mass matrix, see Section 
6 . 2 ) 

NOTE 5 . 3 ~6 : Mass matrix updating affects the velocities of the 

structure. Structural velocities will be 
automatically altered so as to conserve the momentum 
of the system. See section 6.2 for further 
discussion. 
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TIMSPAN 


This is a required comma nd which supplies the starting and 
terminating times of the transient response. Data takes the 
form, 

TIMSPAN , tl , t2 
where, 

tl is the starting time for the transient response 
t2 is the final time for the transi ent . response 

Exampl e_5 j Setting the t ime interval to be a na ly zed 

TIMSPAN , 0 . 0, 200 , $ 200 SECONDS OF RESPONSE STARTING 
$AT 0. 

A condition label cannot be used with the TIMSPAN command. 


TIMSTEP 

This is a required command which supplies the time 
integration step and takes the form, 

TIMSTEP, delt ? Cj 

where,' 

delt Is a time increment 

Cjis the j th optional condition label defined 

with a Cj command; see section 5.15 for further 
details 


Example _5 ^3 ~ 4 j S £ 1 1 ^ng_t he _t; i_m e _s t e p 

TIMSTEP , 1 . E-^$ TIME INCREMENT IS 0.0001 TIME UNITS 

NOTE 5_^32§.1 Multiple TIMSTEP commands may be used, but they have 
a condi ti on-order priority as discussed in NOTE 5.15- 
1 of Section 5.15. Thus, the first appearance of a 
TIMSTEP command in the data set cannot have an 
attached condition label. See also Section 4.5. 
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STIFF 


This is an optional comma nd indicating the calculation of^the 
system tangent stiffness matrix. The program does not use the V, ' 
stiffness matrix to predict the time history of response, but it 
may be calculated for a variety of other reasons depending on user 
preference. For example, the user may wish to determine the 
system frequencies or may wish to use the stiffness to examine 
structural or control stability. Data takes the form, 

STIFF, n,p ? Cj 
where, 

n is an integer specifying that the tangent 

stiffness is to be up-dated every n time steps; 
if n=0 then the stiffness matrix is up-dated 
automatically when it is required for calculation 
of system frequencies ■ f 

p is a -parameter specifying that the incremental 
stiffness, arising from the axial load carried by 
the member, is to be included or excluded from 
the stiffness matrix. If p is Y (for yes) the 
incremental stiffness is included in the 
stiffness matrix calculation; if p is N (for no), 
the incremental stiffness is excluded from the 
stiffness matrix calculation. 

Cj is the j th optional condition label defined 
with a Cj command; see section 5.15 for further 
details 

Exampl e_5^3;^5.J Creating the system tangent stif fn ess, 

matrix 

STIFF , 1 000 , Y$ STIFFNESS CALCULATED EVERY 1000 STEPS 

In this example, the tangent stiffness including the 
effect of incremental stiffness is calculated every 
1 000 time steps . 

NOT E_5 Multiple STIFF commands may be used, but they have a 
condition-order priority as discussed in NOTE 5.15-1 
of Section 5.15. Thus, the first appearance of a 
STIFF command in the data set cannot have an attached 
condition label. See also Section . 5 . 
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FREQ 


This is an optional co mm a nd for calculating instantaneous 
frequencies and associated mode shapes of the system. (The method 
of bisection is employed to extract frequencies.) These 
calculations require the mass and tangent stiffness matrices of 
the system. The mass matrix is automatically calculated, but 
since the transient algorithm used in the program does not use the 
tangent stiffness matrix, it is not calculated unless the user 
requests its calculation through the STIFF command. Thus wh en 

program recalculates the matrices whenever frequencies are 
required. Data takes the form, 

FREQ,k,m,n ? Cj 
where , 

k is an integer specifying frequency and mode shape 
calculation and automatic output every k time 
steps 

|m| is an integer specifying the number of lowest 
frequencies or highest frequencies desired; if 
m is negative the lowest lm| frequencies are 
calculated and if m is positive, the highest m 
frequencies are calculated; of m=0, no 
frequencies are calculated, but if m is 'blank, 
then all frequencies are calculated 

n is an integer < lm| specifying the number o,f mode 
shapes desired; n mode shapes associated with 
the lowest frequencies, if m is negative, the 
highest frequencies if m is prositive and all 
mode shapes if n is blank 

Cj is the j th opti onal condition label defin^ed 

with a Cj command; see section 5.15 for further 
details 
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FREQ 


NOTE 5.3 


Exam£l e_5j_3 "ii 


FREQ, 1 500,-1 0, 1 0$ 10 LOWEST FREQUENCIES CALCULATED 

In this example, the 10 lowest frequencies and mode 
shapes are calculated every 1500 time steps. In 
performing this calculation, the program will use 
the most recently up-dated mass and stiffness 
matrices. 

8 ; Multiple FREQ commands may be used, but they have a 
condi t i on- or der priority as discussed in NOTE 5.15-1 
of Section 5.15. Thus, the first appearance of an 
FREQ command in the data set cannot have an attached 
condition label. See also Section 4.5. 
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FREQ 



STOP 


This is an which if used always appears wi th 

a condition lab el . It instructs the termination of the problem 
when the condition label is true and thus may be used to terminate 
job execution before maximum time in TIMSPAN command is reached. 
Data takes the form, 

STOP ? CJ 

where, 

Cj is the j th condition label defined with a Cj 
command; see section 5.15 for further details 

NOTE_5^32l®- Section 5.15 on condition labels provides additional 
information concerning this command and appropriate 
examples. 
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Section 5 . 4 


Beam Material and Structural Properties 


At present only spring, damper and flexible beam members are 
available. Spring and damper properties are described under their 
own subheading. Beam material and structural properties are 
specified using the command MATPROP. This command is re quired if 
beam members are used. Data takes the form, 


MATPROP : matname ,ratype , E, A, I , rho ? Cj 


I 

\ 

( 

i 


where , 

matname is a unique property identification integer 
or alpha-numeric name 

mtype is the type of beam member with four beam 
types being permitted 

mtype=0, flexible member 
mtype=1 , extensionally rigid, but 
flexible in bending 

mtype=2, flexurally rigid, extensionally 
flexible 

mtype=3i fully rigid 

(Presently only mtype =0 and 3 are 
operational . ) 

If mtype entry is blank, the default 
val ue_^3_zero . 

E is Young's modulus 

A is the beam cross sectional area 

I is the beam areamoment of inertia 

rho is the mass density per unit volume 

Cj is the j th optional condition label defined 

with a Cj command; see section 5.15 for further 
details 


NOT E_5 to 99 different beam properties can be specified. 

Any consistent set of units may be used in the 
MATPROP command. 
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e2epoxy_proper t i^es 

MATPROP, ALUMINUM, 0, 10. 5E6, .i7.01,2.588E-J4$ DEFINES 
$ ALUMINUM MATERIAL PROPERTIES ' 
GRAPHITE/EP0XY,0,110.E6, .2,.01,2.E-^$ DEFINES 
$GRAPHITE-EPOXY MATERIAL PROPERTIES 

In this example, two beam material properties are 
defined. The first material corresponds to aluminum 
beam members and the presence of zero in the second 
numerical entry indicates that beam members which 
refer to this MATPROP command are fully flexible. 

The third line of this example does not commence 
with a command name. Since the numerical data of 
the first line does ' not terminate with a & sign, the 
program recognizes that the third line constitutes a 
new command which is of the same type as the most 
previous command, namely, a MATPROP command. (Had 
the previous line of data ended with a & sign, the 
program would have assumed the third line of the 
example to be a continuation line of the first.) 

The second MATPROP command would be referred to by 
beam members which are fully flexible and made of 
graph! t e- epoxy 

Example 5.4-2: Material Properties of Rigid Members 

MATPROP , 3 , 3 , , , , 2 . E-4 $ MATERIAL 3 REFERS TO FULLY 
$ RIGID MEMBERS 

4 , 3, , , , 5. E-4$ MATERIAL 4 REFERS TO FULLY RIGID 
$ MEMBERS 

In this example, two material properties are 
defined. Beam members which refer to these property 
commands are fully rigid. 

N0TE_5^422: There may be only eight MATPROP commands with the 

same identification number and with attached 
condition labels and there must be at least one 
MATPROP command without a condition label for each 
material identification number appearing in the 
input data. , 

I 

NOTE 5.4- 3: Multiple INCMASS commands with the same 

identification number or name may be used, but they 
have a condi tion- or der priority as discussed in NOTE 
5.15-1 of Section 5.15. Thus, the first appearance 
in the data set of an INCMASS command, with a 
particular identification number or name, cannot 
have an attached condition label. See also Section 
4.5. 
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MATPROP 


NOTE 5. M-4; 


When rigid members are employed, MTYPE greater than 
zero, it is advisable to periodically update the 
mass matrix. (See INCMASS command in Section 5.2 
and user guidelines in Section 6.2 for further 
details . ) 

N0TE_5 When the mass property of a member is changed during 
exectuion, the program will automatically update the 
mass matrix and adjust the velocitities of the 
structure to conserve ^system momentum. Energy, 
however, may not be conserved. See section 6.2 for 
“ further discussion. 
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Section 5.5 - Grid Point Locations and Stationary Global 
Coordinate System 

The stationary global cartesian coordinate system (X,Y) is 
depicted in figure 5.5-1. Locations of grid points in this system 
are specified with the GRID command. A single grid point or a 
continuous string of grid points may be specified with one GRID 
command. 

The GRID command automatically generates grid points in an 
equally spaced pattern between two user specified grid points. 

Grid point numbers are automatically assigned sequentially with a 
user supplied increment. Data takes the form, 

GRI D , ns t , nn , X 1 , y 1 , x2 , y 2 , i nc 
where, 

nst is the starting grid point number 

nn is the number of grid points defined by this 
string 

xl and y1 are the x and y coordinates of the first 
grid point in the string 

x2 and y2 are the x and y coordinates of the last 
grid point in the string 

inc is the increment in the numbering of the grid 
points 

To define only one grid point, nn is set to 1, inc is set to zero 
or blank, and x2 and y2 are left blank. If values for x2 and y2 
are given, they are ignored by the program when nn is 1., 

Exampl e_5_^52lj 

GRID, 10;1 1,1 00., 100., 1 000 . , 1 000 . , 1 0$ FIRST 

$STRING 

110, 3,100., 100., 1000., 100. ,10$ SECOND STRING 

In this example, two strings of grid points are 
defined. The first GRID command defines the first 

string with grid point numbers 10, 20, 30 1 00 

and coordinates (100. ,100.), (200. ,200.), 

300 ., 300 .),...,( 1 000 ., 1 000 .) , respectively. The 
second line is also a GRID command and defines the 
second string with grid point numbers 110, 120, 130, 

and coordinates (100.,1Q0.), (500. ,100.) 

( 1 000 ., 100 .). 
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NO TE 5.5-1: An unlimited number of GRID commands may be used, 

but the total number of grid points defined by all 
GRID commands cannot exceed the maximum value 
designated with the PORDER command. At least one 
GRI D command must be used. If more than one GRID 
command defines the same grid point, the coordinates 
of the last command are used. 

A condition label cann ot be used in conjunction with a GRID 
command . 
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Section 5.6 - Beam Member and Local Coordinate Definition 
5 . 6.1 - Beam Members 


Beam members are defined using the MEMBER command. A single 
member or a continuous string of members may be defined with a 
single MEMBER command. The user specifies the end point grid 
numbers' of the first member in the string and the number of 
members defined by the string. The increment between the user 
supplied grid point numbers for the first member of the string is 
assumed to be the same for all members in the string. Members are 
defined to run from the first grid point of the string to the 
second, from the second to the third, from the third to the 
fourth, etc. Thus members have a beginning grid point and a 
terrai nat.i ng ' gr i d point as shown in figure 2. 

Data for the MEMBER command takes the form, 

MEMBER, IMEMB,igridA,igridB,ne,matname 
where , 

IMEMB is the user defined member number for 
the first member in the string 

igridA is the grid point integer number at one 
end of the first member in the string 
that is being defined 

igridB is the grid point integer number at the 
other end of the first member in the 
string that is being defined 

ne is the number of members in the string 
being defined 

raatname is the number or alphanumeric name of 
a material property valid for all 
members in the string being defined; 
material properties are given by the 
MATPROP command 

Members run from end A to end B. 

1 ’ Member numbers in a string are assigned 
^ sequentially. If IMEMB is zero or blank, the first 

member in the string is automatically assigned a 
member number based on the order in which the MEMBER 
command appeared relative to other MEMBER commands. 
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Example 5.6-1: Defining a string of members 


MEMBER, , 1 0, 20, 1 0, 2$ FIRST STRING OF TEN MEMBERS 
ilO , 1 1 0 , 1 20, 2, 3$ SECOND STRING OF TWO MEMBERS 

In this example, two strings of members are defined. 
The first line of data is a MEMBER command defining 
a string of 10 beam members. The first member runs 
from grid point 10 to 20, the second from 20 to 
30,...,the tenth from 90 to 100. Material 
properties of all these members are defined by the 
MATPROP command number 2. Since IMEMB is blank and 
this is the first MEMBER command in the data set, the 
members are automatically numbered from 1 to 10. 

The second line is also a MEMBER command; It' 
defines a string of 2 members, the first running 
from grid point 110 to 120 and the second from grid 
point 120 to 130 with material properties defined by 
MATPROP command number 3. These members are 
numbered ^10 and 41. If the associated MATPROP 
commands were those given in Examples 5.4-1 and 5.4- 
2, then the first string of members are fully 
flexible and the second string of members are fully 
rigid. 

A condition label ca nnot be used in conjunction with a MEMBER 

command. 

N OT E _5 j_ 622 _^ There is no limit on the number of MEMBER commands 

that can be used, but the total number of members 
defined by all such commands cannot exceed the 
maximum value specified by the PORDER command. 

attempt is made to redefine an already defined 
MEMBER, the program will terminate at the pre- 
processor level and a fatal message will appear in 
the output file. 
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5.6.2 


Local Coordinate Systems 


In addition to the stationary global axis system, the user 
should be aware of the rotating local axis systems. Local axis 
systems (X',Y’),see figure 5.5-1, are not define d by the us er , but 
are automatically generated by the program. A local axis system 
is generated for each user defined beam member. This is 
accomplished internally by first establishing a direction for the 
member. In figure 5.5-1, this direction lies along a straight, 
line from end A to end B of the member and is denoted by a unit 
vector V. Thus V changes it's orientation as the member moves in 
space. 

At any instant in time V makes an angle, denoted as 0, with 
the positive global X axis. The program then automatically 
defines local cartesian axes whose origin coincides with that of 
the global axes, but rotate, (not translate), with the vector. 
These axes are denoted as "local axes" or "local coordinate 
system" since they are different for each member. The local X’ 
axis is always parallel to the local member vector and thus runs 
parallel to the member. The local Y' axis is normal to the member 
and makes the same angle with the global Y axis as the local X' 
axis makes with the global X axis. 
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Section 5.7 - Geometry Scaling 

The geometry of the entire structure may be scaled and/or 
rotated through the use of an opt ional SCALE command. This 
command is useful when parameter studies are being performed and 
the user desires an easy way to vary structural size or 
orientation. Data takes the form, 

SCALE , s cal ex , s cal ey rotate 
where , 

scalex is a scale factor for the global x- 
direction; default is 1 
scaley is a scale factor for the global y- 
direction; default is 1 

rotate is an angle in radians for rotating the 

entire structure (positive angle measures 
rotation from the positive global x axis 
towards the positive global y axis); default 
i s 0 

The order of operations is to scale the structure in the 
global x-direction, the global y-direction, and then rotate the 
structure. If "scalex" or "scaley" are zero or blank, a default 
value of 1 is used for the associated scaling in that direction. 

If "rotate" is blank, a default value of zero is used. If a SCALE 
command is not used, no scaling of geometry nor rotation occurs. 

Ex ampl e_5^72ij 

SCALE, 10. ,,0.8$ SCALE X-DIRECTION BY A FACTOR OF 10 
$ AND ROTATE STRUCTURE 0.8 RADIANS 

In this example, the geometry of the entire 
structure is expanded by a factor of ten in the 
global x-direction and then rotated by a positive 
angle of 0.8 radians. 

A condition label canno^ be used in conjunction with a SCALE 
command . 
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Section 5.8 - Spring Member Definition 

Optional linear springs are defined using the command SPRING. 
Springs may be either extensional or rotational and act in either 
a fixed direction or a variable direction. Data takes the form, 

SPRING , i gr i d1 , i gri d2 , idi r , sprngk , psspng ? Cj 

where , 

igridi is the grid point at one end of the 
spring 

igrid2 is the grid point at the other end of 
the spring 

idir is the direction in which the spring 
acts, that is, 

idir=0, spring acts in a direction 
parallel to a straight line 
joining its end grid points and 
thus its line of action varies 
idir=1, spring acts in fixed x 
direction 

idir=2, spring acts in fixed y 
direction 

idir=3, spring acts in rotational 
direction 

sprngk is the spring constant 

psspng relates to the initial stretch or 

rotation in the spring; see equations 
below 

Cj is the j th opti onal condition label which 
is defined with a Cj command; see section 
5.15 for further details. 
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N0TE_5^8-1: The spring force at igridl is given by, 


sprngk X [ d i spl acem en t ( i gi r d2) 

- di spl acement ( i gr i d 1 ) + psspng] 

The spring force at igrid2 is opposite in sign to that 
at igridl. 

The spring torque if dir=3 at igridl is given by, 

sprngk X [rotation ( i grld2) 

- r ot at i on ( i gr i d 1 ) + psspng] 

The spring torque at igrid2 will be opposite in sign 
to that at igridl. 

Ex ampl e_5 Defining a rotational spri ng 

SPRING, 1 0, 1 10, 3, 1 . E3, 0.1$ ROTATIONAL SPRING 

In this example a rotational spring with spring 
constant of 1000 and initial relative rotation of 
0.1 radians, acts between grids points 10 and 110. 

NOTE b.8-1: Up to 99 springs may be defined. 
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Section 5.9 - Damper Member Definition 

Opti onal linear viscous dampers are defined using the command 
DAMPER. The damper may either extensional or rotational and 
either act in a fixed direction or a variable direction. Up to 99 
dampers may be defined. Data takes the form, 

DAMPER , i grid 1 , i grid2 , i dir , cdamp ? Cj 

where, 

igridi is the grid point number at one end of 
the damper 

igridP is the grid point number at the other 
end of the damper 

idir is the direction in which the damper acts 
and is defined exactly as that for a 
spring 

member, (see Section 5.8) 

cdamp is the linear viscous damper constant 

Cj is the j th opti onal condition label which 
is defined with a Cj command; see section 
5.15 for further details 

Exampl e_5j;_9 “ij Defining a linear viscous da mp er 

DAMPER, 1 0, 1 1 0, 1 , 1 00. $ DAMPER NUMBER 10 

In this example, a damper with damping constant 100 
joins grid points 10 and 110 oriented in the x 
direction 


1 1 - 15-85 


5.9-1 


DAMPER 



Section 5.10 - Optional Additions to Inertial Coupling And Lumped 

Masses 

Inertial coupling of degrees-of -freedom is a natural 
consequence of the consistent mass approach used in the program, 
however the user may wish to specify additional inertial coupling. 
This is acccorapl i shed via optional ADMASS commands. ADMASS 
commands may also be used to define lumped masses at any grid 
point. Data takes the form, 

ADMASS , i gr id 1 , i dir 1 , i grid2 , i di r 2 , mass ? Cj 

where, igridi and idirl provide one degree-of-f reedom of 
the structure and igrid2 and idir2 provide a second 
degr eee-of- freedom of the structure; that is, 

igridi is a grid point for the first degree-of- 
freedom 

idirl is the direction for the first degree-of- 
freedom 

igrid2 is a grid point for the second degree-of- 
freedom; for defining lumped masses, nodel 
and node2 should be the same 

idir2 is the direction for the second degree-of- 
freedom; for defining a lumped mass, idirl 
and idir2 are generally the same 

mass is the magnitude of the additional inertial 
coupling term between the two degrees-of- 
freedom 

Cj is the j the optional condition label which is 
defined by a Cj command; see section 5.15 for 
further details 
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Example 5.10-1: Adding mass coupling and lumped mass 


ADMASS, 1 00, 2, 90, 2, 1 .-2$ ADD OFF DIAGONAL INERTIA 
90,2, 100,2, 1 .E-2$ ADD OFF DIAGONAL INERTIA 
1 1 0, 1 , 1 1 0, 1 , 5.E-2$ ADD LUMPED MASS 
1 1 0, 2, 1 1 0, 2, 5. E-2$ADD LUMPED MASS 

In this example, the first ADMASS command results 
in an additional inertial coupling with magnitude 
0.01 specified to act between the degree-of-freedom 
indicated by grid point 100 and direction 2, and 
grid point 90 and direction 2. This will result in 
a value of 0.01 added to the appropriate off- 
diagonal term of the structural mass matrix. 

Though not required, the second ADMASS command 
retains the symmetry of the mass matrix. The third 
and fourth ADMASS commands are used to add a lumped 
mass of 0.05 at grid point 110; that is, a diagonal 
term in the appropriate locations of the mass 
matrix . 

99 ADMASS commands may be used. 

NOTE^5^i022: This optional command permits the inclusion of 

incompressible fluid effects wherein the presence 
of the fluid gives rise to an additional full mass 
matrix. 
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Section 5.11 - Beam Member Preload 

An axial prestress may be applied to a fully flexible beam 
member by use of the op tional PRELOAD command. Data takes the 
form , 


PRELOAD ,m, preload 
where , 

m is a number of a fully flexible beam member; 

i.e., MATPROP command asociated with this member 
. has MTYPE=0 

preload is an initial member load, positive for 
tension and negative for compression. 

Ex ampl e_5^il2ij §P££i£yi.I 18 _fneraber _pr el oad 

MEMBER , 1 , 6 , 7 , 5, 2$ FIVE MEMBERS DEFINED 


PRELOAD , 1 , 1 000 . $ PRELOAD IN MEMBER 1 

2.2000. $ PRELOAD IN MEMBER 2 

3.2000. $ PRELOAD IN MEMBER 3 

n,-500.$ PRELOAD IN MEMBER 

5,-1000.$ .PRELOAD IN MEMBER 5 

In this example, member number 1 which joins grid 
points 6 and 7 has a tensile preload of 1000, 
member 2 has a preload of 2000, member 3 has a 
preload of 2000, member 4 has a preload of -500 and 
member 5 has a preload of -1000. 

A condition label canno^ be used in conjunction a PRELOAD command. 

N OT E _5 1_J_2l • Members for which no PRELOAD command is provided are 
assumed to have a zero pre-load. 
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Section 5.12 - Externally Applied Loads and Gravity 

Externally applied loads may be imposed at any grid point and 
defined to act in any direction in the global coordinate system. 

The time variation of the loads may be specified in a 
piecewise linear form through the use of a user supplied table of 
values. In addition, options are provided for loads which are 
constant or vary sinusoidally with time. 


Loads are defined in terms of load sets. A load set consists 
of a single time varying load with a particular magnitude and 
direction applied at a set of grid points. The following command 
is required for each set of applied loads. 

APPLOAD - to identify a load set and to indicate how, 

where, when and under what conditions, if any, 
a specified time varying load set is applied; 
an APPLOAD command is re quired when any lo ad 
set is defi ned 

The APPLOAD command takes the form, 

APPLOAD , nload , i memb , dir , s cal e , on , of f , i gri d 1 , i gri d2 ,... ? Cj 
where , 

nload is an integer load i dent i f cat i on number 

imemb is the identification number of the member 
whose local coordinate system is used to 
define the direction of the applied load; if 
imemb is zero or left blank, the global 
system is used 
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alpha= direction of load in radians as measured from 
either the 

global X axis (positive angles run from 
positive X to positive y axes) if imemb =0 
or the 

local x' axis (positive angles run from 
positive X' to positive y' axes) of member 
imemb when imemb>0; see figure 5.18-1 

alpha may be a user defined variable Qk, 
where Qk is the kth user defined variable, 
(see NOTE 5.12-1) 

alpha=T0RQUE for applied torques (positive torques 
act from positive global x to positive 
global y axis) 

scale is a factor used to scale the load magnitude 
values given by the LOAD, LOADSIN or LOADTIM 
commands, defined subsequently and if zero or 
left blank, default value of 1 is used; scale 
may be a user defined variable Qk, where Qk 
is the k th user defined variable (see NOTE 
5 . 12 - 2 ) 

on and off are the absolute values of time at which 
the load set is applied and removed, 
respectively 

igridi , igrid2, . . . are the grid points at which the 
load set is to be applied or Ln where n is 
the identifier of a defined list (see page 
5 . 2 - 1 ) 

Cj is the J th optional condition label which is 
defined with a Cj command; see section 5.15 for 
further details 
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NOTE; 5. 12-1 : 


"Alpha" may be specified via a user defined 
variable, Qk, where Qk is the k th user defined 
variable as described in section 5.16. Qk may be 
defined as a FORTRAN expression which depends on 
variables of the structural response such as 
velocity, acceleration, internal forces, etc. For 
example, follower forces can be specified by 
defining "alpha" as a fixed angle between a rotating 
beam member and the applied force. (See example 

5.16- 4.) Table 4.4-1 provides a complete list of 
admissible variables for Qk definitions. 

NOTE 5. 1 2-2; "Scale" may be specified via a user defined 

variable, Qk, where Qk is the k th user defined 
variable as described in section 5.16. Qk may be 
defined as a FORTRAN expression which depends on 
variables of the structural response such as 
velocity, acceleration, internal forces, etc. For 
example, this permits the user to define control 
forces with or without time delay. (See examples 

5.16- 2, 5.16-5 and 5.16-6.) Table 4.4-1 provides a 
complete list of admissible variables for Qk 
definitions. 

Example 5.12-1; APPLOAD Command for Applied Forces 
DEFINE;@6=1 .5707963 

APPLOAD,!, ,§6,5. 0,0. 0,40. 0,10, 15, 20$ APPLIED LOAD AT 
$ GRID POINTS 10, 15 AND 20 

In this example, an applied load set acts at grid 
points 10, 15 and 20. The load direction is 
prescribed to be 90 degrees from the global x axis 
(or parallel to the global y axis). The scale factor 
is 5 and the load set is applied at time 0 and 
removed at a time of 40 units. The first entry of 
the command identifies an additional command which 
identifies the time variation of the load set which 
is described subsequently. 
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Example 5.12-2: APPLOAD Command for Applied Torque 


APPLOAD, 2, .TORQUE, '4.5, 3.0, 12.0, 10, 15,20$ TORQUE AT 
$GRID POINTS 10, 15 AND 20 

In this example, an applied load set consists of 
torques at grid points 10, 15 and 20. The scale 
factor is 44.50 and the load set is applied at 3 
units of time and removed at 12 units of time. The 
first entry of the command identifies an additional 
command which identifies the time variation of the 
load set which is described subsequently. 

Any of the following commands may be used to prescribe the 
time variation of applied loads. 

LOAD - to specify a constant load in time (or a 

system response dependent load when a user 
defined Qk variable is employed as a scale 
factor in the APPLOAD command; see NOTE 5.12-2 

LOADSIN - to specify a sinusoidally time varying load 

LOADTIM - to specify a piecewise linearly time varying 
load 


NOTE 5. 1 2-3= Each of these commands may be referenced by many 
APPLOAD commands. The APPLOAD commands take the 
time variation specified in these commands and 
scale, orient and place them at prescribed grid 
points on the structure. If these commands are not 
referenced by APPLOAD commands, they are inactive. 
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The optional LOAD command takes the form, 


LOAD, nload, amp 
where, 

nload is a unique load identification number 
which is referenced by APPLOAD commands 

amp defines the load value 

Example 5.12-3: Specifing a constant load in time 
DEFINE: @6=1.5707963 $PI/2 DEFINED 

APPLOAD, 1 , ,@6,5.0, 0.0, ilO.0, 10, 15, 20$ APPLIED LOAD AT 
$ GRID POINTS 10, 15 AND 20 
LOAD: 1,20.$ LOAD OF 20 APPLIED 

In this example, a constant load is referenced by an 
APPLOAD command. Thus the load acts in a direction 
90 degrees from the positive global x axis; that is, 
along the positive y axis and acts at grid points 
10, 15 and 20 from time 0 to time of 40 units.. 

From the LOAD command, it has an amplitude of 20. 

The load amplitude is scaled by a factor of 5 to 
produce a total load amplitude of 100. 

A condition label cannot be used in conjunction with a LOAD 
command, but can be used in with an APPLOAD command which 
references it. 


11-15-85 


5. 12-5 


LOAD 



The optional LOADSIN command takes the form, 


LOADSIN , nload , freq .phase , amp 
where , 

nload is a unique load identification number which is 
referenced by APPLOAD commands 

freq is the frequency of the load in hertz 

phase is the phase angle of the sinusoidal load in 
radians 

amp is the amplitude of the sinusoidal load 

Use of LOADSIN will then result in a sinusoidal load from 
time=tloadon to tloadoff namely, 

(amp) sin[2ir(freq) (t-tloadon) + phase] 

where tloadon and tloadoff are prescribed in the APPLOAD 
command which references the LOADSIN command through the use 
of the LOADSIN identification number. 

Example 5.12-A; Sinuoidal Time Variation of Torque 

APPLOAD, 2, .TORQUE, A. 5, 0.0, 12.0, 10, 15,20$ TORQUE AT 
$ GRID POINTS 10, 15 AND 20 

LOADSIN, 2, 2. 0,0. 0,3.0$ SINUSOIDAL LOAD IN TIME 

In this example, a sinusoidal load is referenced by 
an APPLOAD command. Thus a torque acts at grid 
points 10, 15 and 20 from time 0 to time of 12 
units. From the LOADSIN command, it has a 2 Hertz 
i variation, no phase shift and an amplitude of 3- 
The amplitude will be scaled by A. 5 to yield an 
amplitude of 13-5. 
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Example 5.12-5: Cosine Variation of Load 


DEFINE,@1 =3.1 ^159265359$ PI DEFINED 
DEFINE, @2=0. 7853982$ PI/H DEFINED 

APPL0AD,9, ,@2,5.0,0.0,H0.0,10,15,20$ APPLIED LOAD AT 
$ GRID POINTS 10, 15 AND 20 


L0ADSIN,9,2.0,@1 ,3.0$CPSINE VARIATION IN TIME 

In this example, the value of n is defined somewhere 
in the data set. A sinusoidal load associated with 
the first APPLOAD command, having a frequency of 2 
Hertz, phase shift of n radians and amplitude of 3 
is defined. Thus a load with cosine type variation 
is indicated. From the APPLOAD command the time 
duration is from time 0 to a time of HO units and 
the load acts in a direction given by @6 or H5 
degrees from the positive global x axis and is 
applied at grid points 10, 15 and 20. 

The LOADSIN command cannot be used in conjunction with a 
condition label, but an APPLOAD command which references it can 
have a condition label. 
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The optional LOADTIM command takes the form, 

LOADTIM, nload, tt ,k , timel , loadi , time2,load2, . . . , timek, loadk 
where , 

nload is a unique identification integer which 
is referenced by APPLOAD commands 

tt is the time type, either ABS^ if time 

values used in this command are absolute or 
REL if time is measured relative to the 
time the load is applied 

k is the number of pairs of (time, load) values 
used to define a piecewise linear time 
variation 

timel, loadi is the first pair of (time, load) 
entries 

time2,load2 is the second pair of (time, load) 
entries 

timek, loadk is the k th pair of time load 
entries 

Time values, timel, time2, etc., must be given in ascending 
order. Linear interpolation is used between, the time, load pairs 
and linear extrapolation is used outside the range of the 
time, load pairs. 

Example 5.12-6: Piecewise Linear Load in Time 

APPLOAD, 4, , TORQUE, 11. 5, 0.0, 100.0, 10, 15,20$ TORQUE AT 
$GRID POINTS 10, 15 AND 20 

LOADTIM, 11, 6, ABS, 65.0,2.0, 10.0,2.0, 15.0,11.0, & 
20.0,4.0, 100.0,0.0, 1000.0,0.0 $ LINEAR PIECEWISE 
$VARIATION 

In this example, a piecewise linear load variation 
in time is to be applied in conjunction with 
instructions in an APPLOAD command. Thus the load 
set consists of torques acting at grid points 10, 15 
and 20. From the LOADTIM command, six pairs of 
time/load values are given to define the piecewise 
variation of the load. Linear interpolation is used 
between the pairs of values and linear extrapolation 
is used outside the range of the values. Thus, for 
instance, the load value at time 0 is 2 and that at 
time 2000 is 0. 
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A condition label cannot be used in conjunction with a 
LOADTIM command, but an APPLOAD command which references it can 
have a condition label. 
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The optional GRAVITY command uses the well known acceleration 
at earth’s sea level to set the gravity constant in Newton’s 
Gravitational Law. Gravity forces are applied to all grid points. 
To account for variations in gravity force throughout the 
structure and thus gravity gradient, the coordinates of the 
earth’s center in the global axis system is required as input. 

Data takes the form, 

GRAVITY,g ,xearth,yearth 
where , 

g is the acceleration due to gravity at sea level 
xearth is the global x axis location of the earth’ s 
center 

yearth is the global y axis location pf the earth's 
center 

Example 5.12-7: Specifying a Gravity Field 


GRAVITY,32.2,0.0,-2. Ile6$ fflAVITY FIELD DEFINED 

In this example, the units being used in the problem 
are feet and seconds. The acceleration due to 
gravity is 32.2 feet per second squared and the 
coordinates of the earth’s center in the global axis 
system are x=0 and y=-2.11e7 feet (or 3996 miles). 
The gravity load is applied to all grid points. 

A condition label cannot be used in conjunction with a 
GRAVITY command. 

NOTE 5. 1 2-^ ; Since the GRAVITY command places gravity forces at 
all grid points and orients them toward the earth’s 
center, an APPLOAD command is not used to reference 
a GRAVITY command. 

NOTE 5.12-5 : The GRAVITY command is not presently available. 

However, gravity can be approximated through the use 
of SPRING commands as done in the following example. 
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Example 5.12-8; Gravity Load Applied Through Uae of 
SPRING Commands 

SPRINGzlOO, 1 ,0,-e3.@4 $SPRING TO SIMULATE GRAVITY $ 
AT GRID POINT 1 

In this example, the earth's center is placed at 
grid point 100. The spring acts on a straight line 
joining grid point 100 and grid point 1 which is on 
the structure. The spring has an inital stretch 
which simulates the gravity force prior to 
structure deformation and motion, and a spring 
constant which simulates gravity gradient. The 
initial stretch and spring constant are specified 
with DEFINE commands using the formulas: 

psspng = (initial stretch)/ (spring constant) 

= [gm(Rg/R^)^]/[-2gm(Rg/R^)^/R^]= -R^/2 

2 

sprngk = spring constant = -2gm(Rg/R^) /R^ 

where g is the acceleration due to gravity at 

sea level, m is a quantity of mass lumped at the 

grid point to simulate the mass of all members 

and additional mass attached to that grid point, 

R is the radial distance from the earth's 
o 

center to the initial position of the grid point 

and R is the earth's radius, 
e 

These formulas represent a first approximation to 
gravity gradient phenomena. 
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Section 5.13 - Constraints 


The following types of constraints are currently available: 

(1) single degree-of-freedora constraints 

(2) linear multi-degree-of-freedora constraints 

(3) pin connection between two grid points 

(4) rigid member constraints 

(5) inequality constraints . ■ 

5.13.1 ~ SPFC (Single Degree-of-Freedom Constraint) Command 

Single degree-of-f reedom constraints enforce the motion at a 
specified grid point and direction to be a given fixed value 
through the use of a SDFC command. Data takes the form, 

SDFC,igrid,idir,rhs ? Cj 
wher e , 

igrid is the constrained grid point number 

idir is the contrained direction (idir=1,2 or 
3 constrains motion in the global x, y or 
rotational direction, respectively) 

rhs is the constant value to which the 
constrained degree-of-freedom is set 

Cj is the j th optional condition label which 
is defined with a CL command; see section 
5.15 for further details 

Example 5.13~1: Specifying a single degree-of-freedom 

constraint 


SDFC, 10, 1 ,0.0$ CONSTRAINT AT GRID POINT 10 

In this example, grid point 10 is constrained to 
zero in the direction parallel to the global x axis. 
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Note 5« 1 3~1 : The grid point number and direction specifications 

of the SDFC command are translated by the program to 
a particular problem degree-of-freedom. This 
degree-of-freedora , unless otherwise specified 
through the use of a INDDOF command, will be treated 
by the program as a dependent motion and will be ' 
eliminated from the problem. (Section 5.13.6 
describes the use of the the INDDOF command.) The 
reaction associated with the degree-of-freedom may 
be recovered for printing through the use of the 
REACT command described in Section 5.18. Reactions 
may also be used in condition statements or Q- 
variables as described in sections 5.15 or 5.16, 
respecti vley . 


11-15-85 


5.13-2 


SDFC 



5.13. 2 - Multl-Degree-Of-Freedom Constraints 

Multi-degree-of-freedom commands enforce a linear 
relationship between a collection of structural degrees-of- 
freedom. This is accomplished through the use of the optional 
MDFC command. Data takes the form, 

MDFC.k , igridi , idirl ,coef 1 , igrid2, idir2, coef 2, . . . ,rhs?Cj 
where , 

k is the total number of degrees of freedom 
appearing in the MDFC command 

igridi is the grid point number of the first 
degree-of-freedom in this command 

idirl is the direction of the first degree-of- 
freedom in this command; idirl =1,2 or 3 
for direction parallel to global x axis, 
parallel to global y axis or rotation, 
respectively 

Goefi is the coefficient for the first degree- 
of-freedom in this command 

igrid2 is the grid point number of the second 
degree-of-freedom in this command 

idir2 is the direction of the second degree- 

of-freedom in this command; idir2=1,2 or 
3 for direction parallel to global x 
axis, parallel to global y axis or 
rotation, respectively 

coef2 is the coefficient for the second 
degree-of-freedom in this command 

rhs is the right hand side of the linear 
multi-degree-of-freedom command 

Cj is the j th optional condition label which 
is defined with a Cj command; see section 
5.15 for further details 
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The linear multi-degree-of-freedora constraint takes the 
equation form, 

(coef 1 ) (dof 1 ) + (coef 2) (dof2) + ... + (coef n) (dof m) = rhs 

if Cj is TRUE 


where , 

dofi, dof 2, . . . ,dof m are the degrees-of-freedom 

associated with a grid point and direction 
appearing in the MDFC command 

Example 5.13~2; Multi-Pegree-of -Freedom Constraint 
Involving Three Degrees-of -Freedom 

MDFC, 3,30,2, 1 .0,35,2, 1 .0,i<0, 2,-1 .0,0.0 

In this example, the motion of grid point 30 in the 
global y direction plus the motion of grid point 35 
in the global y direction minus the motion of grid 
point 40 in the global y direction is constrained to 
be zero. 

Example 5.13~3: Single-Degree-of-Freedom Constraint 
MDFC, 1,10, 1 , 1 .0,0.0$CONSTRAINT AT GRID POINT 10 

In this example, grid point 10 is constrained to 
zero in the direction parallel to the global x axis. 
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NOTE 5.13-2: 


Grid point number, direction and axis system 
specifications are translated by the program to a 
particular problem degree-of- freedom. As the 
program processes the problem constraints, it 
automatically selects a dependent degree-of-freedom 
from each constraint equation when that constraint 
is processed. In turn that dependent degree-of- 
freedom is automatically eliminated from all 
currently unprocessed constraints and from the 
system equations of motion. When a constraint is 
processed, the dependent degree-of-freedom selected 
is that currently having the largest coefficient in 
the constraint equation, but not one which the user 
has declared to an independent degree-of-freedom 
through an INDDOF command. However, the 
coefficients of the constraints generally change 
during the constraint processing procedure and so 
the user cannot usually predict which degree-of- 
freedom will be treated as dependent in an MDFC 
command prior to program execution. The user will 
know though, that degrees-of -freedom specified in 
an INDDOF command are independent degrees-of -freedom 
and will not be eliminated from the problem. (See 
section 5 . 13*6 for further details on use of the 
INDDOF command.) 

Note 5.13~3: The reaction force associated with the dependent 

degree-of-freedom of an MDFC constraint can be 
recovered through the use of a REACT command as 
described in section 5.1i<, however, the user will 
usually not know a priori which degree-of-freedom 
has been treated as dependent. Of course when only 
one degree-of-freedom appears in the MDFC command it 
necessarily becomes a dependent degree-of-freedom, 
unless it is listed In an INDDOF command. See 
Section 5.13.6 for further details on use of INDDOF 
commands. 

NO TE 5. 1 3-^ : It is quite possible that some constraint equations 

may not be independent of other constraint 
equations. The program will automatically recognize 
and handle such a condition. The user will also be 
informed of its occurence. 
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5.13.3 - Pin Connected Grid Points 


«■ 


Pin connections between grid points can be accomplished by 
the use of two MDFC commands which would constrain the motions of 
the two points parallel to the X and Y global axes to be 
identical. If many pin connections exist in a structure, the 
composition of the required MDFC commands can become tedious. To 
expedite the input of pin connection specif cations , an optional 
PIN command may be used. This command results in the automatic 
generation of the appropriate MDFC commands. The PIN command 
takes the following form, 

PIN, igridi , igrid2 ? Cj 

where , 

igridi and igrid2 are the grid point numbers where 

the pin connection exists 

Cj is the j th optional condition label which is 
defined by a CJ command; see section 5.15 for 
further details 

Example 5.13~^: Pin Connection Using PIN Command 


PIN, 10, i|0$PIN CONNECTION BETWEEN GRID POINTS 10 AND 
$ 40 

In this example, a pin connection is specified 
between grid points 10 and 40. Consequently, the 
motions of the two points parallel to the global x 
and y axes are identical, but the grid points may 
rotate freely of one another. 
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-_R£g£d_Member_Cons traj^nts 

Rigid members impose constraints between the motions of their 
end grid points. These constraints are automatically accounted 
for when rigid members are specified in a MATPROP command; see 
Section 5.4. The user need not concern himself with defining 
these constraints. 
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5.13*5 - Inequality Constraints 


Inequality constraints may be applied by the user. Such 
constraints specify that a linear or nonlinear algebraic 
combination of variables must be less than or greater than some 
other linear or nonlinear combination of variables. Inequality 
constraints are applied through the use of ' cond it i on labels. See 
for instance examples of section 5.15. 


11-15-85 5 . 13-8 Inequality Constraints 



5.13.6 


User Selection of Independent Degrees-of -Freedom 


In general, the program processes constraints by eliminating 
certain degrees-of -freedom from the problem. Such eliminated 
degrees-of-f reedom become dependent degrees-of-f reedom. These are 
selected by the program during execution and the user has little 
control over which degrees-of-f reedom will be treated as 
dependent. However, the user can specify which degrees-of -freedom 
must be treated as independent by the program; that is, cannot be 
dependent degrees-of -freedom , through the use of the optio na l 
INDDOF command. 

The INDDOF command takes the form, 

INDD0F:k,igrid1,idir1,igrid2,idir2,... 

where, 

k is the number of degrees-of-freedom specified 
independent by this command 

i gr i d 1 , i gr i d2 , . . . are the grid point numbers associated 

with independent degrees-of-freedom 

idirl , idir2, . . . are the directions in global axis system 

associated with independent degrees-of- 
freedom at grid points node 1 , node 2 , . . . , 
respectively 

E X a m p 1 e _5 1_3 2,5 j §P®£i£i£§.tion_of_^n dependent degrees -of - 

freedom 


INDD0F:3, 10, 1 ,4,2,25,3 $ INDEPENDENT DOF'S 

In this example, three degrees-of-freedom are 
specified which the program must treat as 
independent. They are as follows; 

grid point 10 in direction parallel to global x 
axis 

grid point 4 in direction parallel to global y 
axis 

grid point 25 in rotational direction 

A condition label cannot be used in conjunction with a INDDOF 
command . 
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Not e 5.13-6: INDDOF commands override any other specifications so 

that an attempt to declare a constraint on degrees- 
of- freedom specif ied in an INDDOF command are 
ignored by the program. A warning message will be 
provided in the output file to indicate this 
occurrence. 
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Section 5.1^ 


Initialization 


Initial displacements and or velocities in the global 
coordinate system may be defined by the user. In itial conditions 
on all degrees-of -freedom are assumed to be zero unless they are 
defined through the use of PIS a nd VEL commands. Data takes the 
form, 


INIT , scaledis , scalevel 
DIS,igrid,idir, coord, dis 
VEL, igrid,idir, coord, vel 

where , 

scaledis and scalevel are scale factors for 
all initial displacements and 
velocities respectively; if the INIT 
command is not used, the default 
values are unity 

igrid is the grid point of the degree-of- 
freedom to be initialized 

idir is the direction of motion that is to be 
initialized; idir=1, 2 or 3 for motions 
parallel to the global x axis, parallel 
to the global y axis or positive 
rotations 

coord is the identification number of a member 
in whose local coordinate system the 
initial condition(s) are described; if 
coord =0 or is left blank, the global 
system is implied 

dis and vel are the initialized values of the 
displacement and velocity 
r especti vely 

Ex ampl e_5^i^-]_j Specif i^cat^on_of_^n£t^al_cond it ions 

DIS , 20 , 1 , , 5 . 2 
VEL, 20, 2, 0,-1 .0 

In this example, the displacement at grid point 20 
in the global x direction is initialized to 5.2 
while the velocity of the same grid point in the 
global y direction is initialized to -1. Since no 
INIT command is used the scale factor default of one 
isassumed. 
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general philosophy used in the program is to 
allow the program to execute if at all reasonably 
possible. Especially is this philosophy important 
in the processing of initi.al conditions and 
constraints. (See also section .) In 
speci.fiying initial conditions it is quite likely 
that the user may impose conditions which are- 
contradictory to user specified constraints, 
including rigid member constraints. In such a case, 
the problem constraints will override the user 
supplied initial conditions. The program will 
employ some of the user supplied initial conditions 
and calculate the others so that the constraints are 
satisfied. The program institutes a priority for 
deciding which user supplied initial conditions are 
employed and which are calculated from the 
constraints. 

Priority is on a first appearance basis, so. that 
when a contradiction between initial conditions and 
constraints occurs the program discards the last 
user supplied initial condition and checks to see if 
a con t r adi c i t i on still exists. If there is no 
longer a contradiction, program execution proceeds. 
If a contradiction still exists the next to last 
supplied initial condition is discarded and a 
contradiction check made again. The process 
continues until the contradiction is eliminated. 

(In some cases it may not be possible to eliminate 
the contradiction. For example, when only one non- 
zero initial condition exists and the degree of 
freedom it is associated with is constrained to 
zero.) Warning messages will appear in the output 
file indicating which initial conditions have. been 
discarded and the results at the starting time will 
reveal the effect of the constraints on the initial 
conditions. 


0 
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The user will find this procedure to have certain 
advantages. For example, if degrees-of -freedom are 
slaved to one another through constraints, the user 
need not unduly be concerned about what could be a 
very complicated kinematics problem in determining 
the appropriate complete and consistent set of 
initial conditions. The user may decide not to 
specify the initial conditions at those locations 
which are believed to be slaved, via the 
constraints, to other locations where initial 
conditions are specified. The program will first 
discard the defaulted initial conditions as it 
attempts to produce a well defined problem. The 
effect is that the kinematics of the problem are 
automatically taken care of by the program. 


1.1 -1 5-85 


5. 1 M-3 


INIT, DIS & VEL 



Section 5.15 - Condition Statements and Labels 

Condition statements are opti onal logic expressions written 
by the user in FORTRAN-like syntax, which at any time step may be 
TRUE or FALSE. A condition statement is identified by a condition 

label. The condition statement is specified through the use 

of a CLj (or Cj ) command. The command takes the form, 

CLJ , FORTRAN expression 

where j is an integer identification number. 

Associated with the condition statement is a condition 

label, denoted Cj . LATDYN commands are made conditional by 
attaching a question mark, " ? " , followed by a condition label Cj 
to the end of the command. If the condition statement associated 
with the condition label is TRUE then the data commands to which 
the condition label is attached are active. If the associated 
condition statement is FALSE then the data command to which the 
condition label is attached are inactive. 

I 

Condition labels may be attached to the following commands: 

ADMASS 

APPLOAD 

DAMPER 

FREQ 

INCMASS 

MATPROP 

MDFC 

PLOTINC 

PRINT 

SDFC 

SPRING 

STIFF 

STOP 

TIMSTEP 

SET 

RESET 

KEEP 

ROTLOCK 

In general, the use of such conditions gives the user 
considerable flexibility in executing the program. For example, 

Cj commands may be used to construct inequality constraints by 
defining a unique condition label to be associated with an 
inequality condition statement written in FORTRAN and attaching 
the condition label to a constraint command. 
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Whenever condition labels activate or deactivate commands 
which can affect the mass and/or stiffness matrix, these matrices 
are automatically updated.. A message in the output file will 
indicate this occurrence. Condition labels oh the following 
commands can lead to updating of mass and/or stiffness matrices: 

ADMASS (mass matix updated) 

DAMPER (mass matix updated) 

SPRING (stiffness matrix updated if STIFF command is 
used) 

INCMASS (mass matrix updated) 

MATPROP (mass matrix updated and stiffness matrix also 
updated if STIFF command is used) 

MDFC and SDFC (mass matrix updated and stiffness matrix 

also updated if STIFF command is used) 
STIFF (stiffness matrix updated) 

TIMSTEP (mass matrix updated if DAMPER command is used 
or INTGTYP is i 5 

The logic expressions in the CJ command are written in 
FORTRAN-like syntax using a special set of variables and 
parameters related to the variables of the problem. The variables 
may be displacements, velocities, accelerations, internal and 
external forces, axial strains, member rotations, structural, 
frequency and time. 


The admissible variables in a condition statement are listed 
in Tables 4.4-1 and take the form, 


XLOC ( n) 
YLOC(n) 
ROT ( m ) 


X(p ,n) 


Y(p ,n) 


for the current location of grid point n 
parallel to the global x axis 

for the current location of grid point n 
parallel to the global y axis 

for the current rotated position of member m 
measured positive in a direction from the 
positive global x to the positive global y 
axis 

for component of motion quantity 
(displacement, if p=D; velocity, if p=V; 
acceleration, if p=A) at grid point n parallel 
to the global x axis 

for component of motion quantity 
(displacement, if p=D; velocity, if p=V; 
acceleration, if p=A) at grid point n parallel 
to the global y axis 
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XP(p,n,l) 


YP(p,n,l) 


PHI ( n) 


MAG ( p , n ) 


ANGLE(p ,n) 


AX(tn) 

SHA(m) , 
SHB(tn) 

BM ( m , i end ) 

FEXTX( j ) , 
FEXTY ( j ) 

FEXTM( j ) 

FEXTA( j ) 

STRAIN ( m ) 


for component of motion quantity 
(displacement, if p=D; velocity, if p=V; 
acceleration, if p=A) at grid point n parallel 
to a local x'axis (parallel to a designated 
member , 1 ) 

for component of motion quantity 
(displacement, if p=D; velocity, if p=V; 
acceleration, if p=A) at grid point n parallel 
to a local y'axis (transverse to designated 
member , 1 ) 

for rotation at grid point n in the positive 
angular direction; from positive global x axis 
to positive global y axis 

for magnitude of motion quantity 
(displacement, if p=D; velocity, if p=V ; 
acceleration, if p=A) atgrid point n 

for angular orientation of motion quantity 
(displacement, if p=D; velocity, if p=V; 
acceleration, if p=A) at grid point n; 
positive angles are measured from positive 
global X axis towards positive global y axis 

for axial force in designated member m 
(positive for tension) 

for shear force at end A or B, respectively, 
of beam member m 

for bending moment at A, iend=0, or end B, 
iend =1 , of beam member m 

t h 

for component of j external force parallel 
to the global x or y axis, respectively 

for the magnitude of the j external force 

t hi 

for angular orientation of the j external 
force; positive angles are measured from the 
positive global x to positive global y axis 

for axial strain in member m 
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REACTX(n) for reaction force at grid point n in the 

global x-direction 

REACTY(n) for reaction force at grid point n in the 

global x-direction 

REACTT(n) for reaction torque at grid point n in the 

global x-direction 

REACTM(n) for the magnitude of the reaction force at 

grid point n 

REACTA(n) for the angular orientation of the angular 

force at grid point n 

REACTC(n,l) for the component of the reaction force at 

grid point n parallel to the local coordinate 
system of member 1 

T f or t ime 

NTSTEP for time step number 

OMEGA(k) for frequency- in radians per unit time of 

designated mode k, where the absolute value of 
k designates the k th lowest mode if k is 
negative and the k th highest mode if k is 
positive 

t h 

CLj for the j logical condition; provides 

logical TRUE or FALSE for condition label Cj 

where , 

p is a parameter indicating the type of motion desired 
as , 

D, for displacement 
V, for velocity 
A, for acceleration 

1 designates the member number associated with a local 
axis system 

m is a member number 

n is a grid point number 


CLj 
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The defined variables may be combined as logical expressions 
writteii in FORTRAN using the following logical operators: 


. EQ . 
.GT . 

. LT . 

. GE . 

. LE . 

, AND . 
.OR . 


equal to 
greater than 
less than 

greater than or equal to 
less tha or equal to 
logical "and" 
logical .inclusive "or" 


In addition to the variables defined above, the following 
constants may also be used: 


PI for TT, PIH for ir/2, PIQ for tt/4 and PI2 for 2tt. 


CLj 
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The program uses the condition statement to automatically 
create a FORTRAN subroutine. The subroutine employs the user 
given algebraic expressions of the condition label exactl y as 

legal FORTRAN operation is available to the user 
in writing the algebraic expressions of the condition labels. For 
example, the algebraic expressions may make use of the intrinsic 
FORTRAN in-line functions ABS, TABS, SQRT, SIN, COS, TAN, etc. If 
FORTRAN errors exist in the condition label expressions, compiler 
errors will result. Hence if compilation errors arise it is 
suggested that the user first examine all FORTRAN input 
expressions . 

Ex ample_5 Condition st atement used to constru ct a 

nonlinear spring 


X(D,iO $ DEFINES Cl 
A 

end^of FORTRAN expression 

In this example, the condition statement to be 
associated with a condition label Cl is defined. 

The condition is true when the displacement in the 
global X direction at grid point 6 is greater than 
the displacement in the global x direction at grid 
point This condition statement may be applied 

to any of the data commands listed above by simply 
attaching a question mark, "?", followed by the 
condition label Cl to the right hand end of the 
data command such as, . 

SPRING , 4 , 6 , 1 , 1 0 . 0 , 0 . 0 ? Cl ^SPRING WITH CONDITION 
$ LABEL 

This command defines a spring between grid points 4 
and 6 which acts in the global x direction only 
when condition Cl is true. Effectively a nonlinear 
spring which can only take tensile strains is 
defined. 


CL1 , X(D , 6) .GT . 
A 

start of FORTRAN expression 
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N0TE_5^15-l£_ In many problems it is useful to have multiple 

commands of a single type some of which may have 
attached condition labels. The use of multiple 
occurrences of TIMSTEP, INCMASS, STIFF, FREQ, , 
MATPROP, PRINT AND PLOTINC, allows the user to vary 
the data supplied by these commands during the 
execution according to rules the user specifies 
through condition labels attached to these 
commands. 

Which of the multiple commands is active at any 
time during execution depends upon the sequence in 
which the commands appear in the data set and 
whether their attached labels are TRUE or FALSE. 
Throughout the manual this is often referred to as 
condition-order priority. 

C ondi t i on- or der priority implies that the program 
examines the last command of such a multiple set. 

If it has no attached condition label, then it is 
always active, whereas if it has a condition label, 
it is active only if the condition is TRUE. If the 
condition is FALSE then the program examines the 
preceeding command of the multiple set in the same 
fashion. The process continues until an active 
command is found or all commands in the multiple 
set are exhausted. Thus, the last command of a 
multiple set with a condition label has priority if 
the condition is TRUE. (See subsequent examples 
5.15-2, 5.15-7, 5.15-8 and 5.15-9.) For 
convenience in using this document, condi tion-order 
priority is defined wherever commands which depend 
on this process appear in the text. 


CLj- 
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Example_5^i5-2£__Condi t £on_s t at ernent_u3ed_to_change_t i^me 

st ep d u r I n g J o b e X e.c u 1 1 6 n 

CL2, T .GE. 5.$DEFINES C2 

TIMSTEP,.01 $ INITIAL VALUE OF TIME STEP 

TIMSTEP , /OOl? C2$ DECREASED TIME STEP WHEN C2 TRUE 

In this example, multiple TIMSTEP commands are 
used. Because the order of their appearance in the 
data set is important, no contradictory 
instructions to the program result. A condition 
statement with label C2 is defined and then applied 
to the second TIMSTEP command. The first TIMSTEP 
command sets the time step to 0.01, however, the 
subsequent TIMSTEP command when active, will 
supersede it. Thus, when 5 time units are reached, 
condition C2 is TRUE, the second TIMSTEP command is 
activated and the time step is changed to 0.001. 
Effectively, the variation of the solution time 
step is thereby controlled by the user. 

ne qual_i^y_cons t r a£nt__cons t r uct ed_w l^th 
co n diti on st at ement and single-degree- 
of-freedom constraint 


CL3, Y(D,10) GE. 2.$DEFINES C3 

SDFC , TO , 2 , 2 . 0 ? C3$ CONSTRAINT ACTIVE WHEN C3 TRUE 

In this example, condition statement with label C3 
is defined and then applied to an SDFC command. 

The displacement of grid point 10 parallel to the 
global y axis is constrained to be less than 2. If 
grid point 10 attempts to displace from its initial 
position more than 2 units parallel to the global y 
axis, a single degr ee-of-f reedom constraint is 
applied which fixes the global y displacement at 2 
units. Effectively, the user has constructed a 
capture mechanism at y=2 which captures grid point 
1 0 


CL j 
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E X am p 1 e _5^152i : LH®, 

be hav io r constructed from condition 

£££££!2£££_£££_£l.£ii£l^£l£££l££l££££^£!!! 

constraint 


CL4, Y(D & $ LINE ENDS WITH & SYMBOL THUS IS 
$ CONTINUED ON NEXT LINE 

,10) .GE. 2.0 .AND. REACTY(IO) .LE. 0.0$ DEFINES C4 
SDFC , 1 0 , 2 , 2 . 0 ? C4$ CONDITIONAL CONSTRAINT AT GRID 
$ POINT 10 IN Y DIRECTION 

This example is similar to that of the previous 
example, 5.15-3. However, whereas in the previous 
example, grid point 10 was captured permanently by 
some mechanism at y=2, in this example, grid point 
10 is only constrained at y=2 if the reaction force 
at grid point 10 is negative. Effectively, the 
user has constructed a wall at y=2 and allowed grid 
point 10 to rebound from that wall when the force 
between the grid point and the wall is pulling 
rather than pushing on the structure at that point. 
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To accomplish this rebound behavior it was 
necessary to use the admissible variable 
REACTY which recovers the reaction force at the 
constraint;, that is, in the y direction at grid 
point 10. Prior to impact, C^l is FALSE since the 
displacement of grid point 10 parallel to the 
global y axis is less than 2 and the SDFC command 
is inactive. Consequently, the reaction force is 
zero bef or e impact . When grid point 10 impacts the 
wall, the reaction force is still zero so that 
will now be TRUE and the SDFC command will be 
activated. Thereafter, the reaction force will 
continue to be monitored. ' If the reaction force 
changes sign, C^ becomes false, the SDFC command 
will be deactivated and rebound will occur. 

Exampl e_5^15 ^5 j t at ement__used_to_def i ne_early 

jo b termination 

CL5 ,MAG( V, 20) . GT . I.E^I .AND. T . GT . 2.0$ DEFINES 

$ C5 

STOP? C5$ JOB TERMINATES WHEN Cb IS TRUE 

In this example, a fifth condition statement is 
defined which is true when the magnitude of the 
velocity vector at grid point 20 exceeds 10,000 and 
time exceeds 2 units. The condition label is 
attached to a STOP command which will terminate 
execution if the fifth condition statement is true. 


CLj 
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Ex ampl e 5 j^ 15“6 ] used to control ma ss 

INCMASS.lOO? C6$ MASS MATRIX UPDATE WHEN C6 IS TRUE 

CL6,(PHI(V,1)+PHI(V,2)+PHI(V,3)+PHI(V,4)+PHI(V,5)& 

+ PHI( V,6)+PHI(V,7)+PHI(V,8)+PHI(V,9) & 

+ PHI( V, 1 0) )/1 0 . .GT. 1 000.$ C6 DEFINED 

In this example, the mass matrix is not usually 
updated since an INCMASS command is not usually 
active. Only when the condition statement 
associated with label C6 is true, is the INCMASS 
command active. The condition statement depends on 
the average rotational velocity of the first ten 
grid points of the structure. For additional 
information on mass matrix updating see section 
6 . 2 . 

Ex ampl e_5^i52X2 £oud^t i^on_s t at em en t used to desig na te 

application of applied for ce 

APPLOAD,! , ,0.75, 0.0, 1.0, 30$ SCALE FACTOR IS 

$ USUALLY ONE 

APPLOAD , r ,, 0 . 75 , 0.0, 2.0, 30? C7 $SCA LE, F A CTOR TWO 

$ WHEN C7 TRUE 

L0ADSIN,1, 10.0, 0.0, 15.2$ TYPE OF LOAD IS 

$ SINUSOIDAL IN TIME 

CL7,MAG(A, 30) . LT . 1 . E.3 $ DEFINES C7 

In this example, a sinusoidal load in time is 
applied to grid point 30 with amplitude 15.2 in a 
direction 0.75 radians from the positive global x 
axis. However, when the acceleration at grid point 
30 drops below a magnitude of 1000, the applied 
load scale factor is doubled; that is, the 
sinusoidal load amplitude becomes 30. M. 
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Exam£i e_5_^i_5 “8j 

MEMBER , 1 2 , 1 2 , 1 3 , 1 , 1 , 0$ THIS IS THE 12 TH DEFINED 
$ MEMBER 

MATPROP, 1 , 10.5E6, .1 , .O'l ,2.588E-4$ NOMINAL MATERIAL 
$ PROPERTY 

MATPROP, 1; 0.0, .1, .01, 0.0? C8$ FAILED MEMBER 

$ PROPERTY 

CL8, (ABS(AX(12)) .GT. 1.E4)$ DEFINES C8‘ 

In this example, the 12 th member is associated 
with the first material, property. Two MATPROP 
commands are designated as material property number 
1 . The second appearing- MATPROP command is 
associated with the removal of a member. The 
material property elastic modulus and density are 
made to vanish when the internal axial force of 
member 12 exceeds a failure value of 10,000 as 
defined by the eighth condition label. Thus the 
axial force determines the failure of the member 
and the member is effectively removed. 

Example 5.15~9: Multiple condition statments used to 

remove failed member 

MEMBER , 1 2 , 1 2, 1 3 , 1 , 1 , 0$ THIS IS THE 12 TH DEFINED 
$ MEMBER 

MATPROP , 1 , 1 0 . 5E6 ,. 1 ,. 01 , 2 . 588E-4$ NOMINAL MATERIAL 
$ PROPERTY 

MATPROP, 1,' 0.0, .1, .01, 0.0? C8$ FAILED MEMBER 

$ PROPERTY 

CL8, ( ABS(AX( 1 2) ) .GT. 1.E4)$ DEFINES C8 
MATPROP , 1 , 0 . 0 ,. 1 ,. 0 1 , 0 . 0? C9 $ FAILED MEMBER 
$ PROPERTY WHEN TIME EXCEEDS 2 
CL9,T .GT. 2.0$ DEFINES C9 

This example is similar to 5.15-8 except that here 
an additional condition statement with label C9 is 
also used. Since the MATPROP command with the C9 
label attached appears after that with the C8, label 
attached, the statement associated with C9 will 
take precedence. Thus the member stiffness will be 
reduced to 5.25E6 or half its nominal value after 
time exceeds two units whether or not the axial 
load has exceeded 10,000. Thus the failure 
condition with label C8 attached will have no 
relevance after a time of two units. 
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The above examples display some of the versatility in using 
condition labels such as in creating nonlinear members which act 
only in tension (or compression), controlling time steps, setting 
lock-up and rebound conditions, changing mass matrix update 
increment, instructing possible early job termination and removing 
failed members. They may also be used to apply control forces 
which depend on structural deformation quantities, contol printed 
and plotted output, to add mass to the structure due to impact 
with a foreign object and a host of other functions too numerous 
to 1 i St . 
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Section 5.16 - User Defined Variables (also denoted as Q- 

variables) 

5.16.1 - General Discussion 


In section 5.4 the command DEFINE was introduced. This 
allowed the user to define in a shorthand fashion, constants which 
often appear in the data set. In addition, the program allows the 
user to define variables which depend on the time varying 
quantities XLOC, YLOC, ROT, X, Y, XP, YP, PHI. MAG, ANGLE, AX, 

SHA, SHB, BM, FEXTX, FEXTY, FEXTM,FEXTA, OMEGA, STRAIN, NSTEP and 
T as defined in section 5.15. This enables the user to vary 
certain quantities appearing in data commands subject to the time 
varying quantities of the analysis. User defined variables may be 
used in any condition label, in a TIMSTEP command, in a MATPROP 
command and as a scale factor and/or load direction in APPLOAD 
commands. The use of user defined variables in TIMSTEP o A T P R 0 P 
commands is not available in the current version. 
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5^1 6 . 2 9 . C .® 

User defined variables are prescribed with the opt i on al SET 
and RESET commands. 

To define a constant parameter Qn (which once calculated 
remains fixed), the SET command is used. Data takes the form, 

SET , Qn=FORTRAN arithmetic expression ? Cj 
where , 

n is an identification integer 

Cj is the j th op tional condition label which is 
defined with a CLj command as discussed in 
section 5.15; C 2 _£|}oul d_not_cont ai n_a_us er 


The SET command creates a constant denoted as Qn where n is 
the identification integer of the SET command. Qn is initally set 
to zero. If the SET command has an attached condition, Qn is set 
to a fixed value equal to the arithmetic expression of the SET 
command when the condition first becomes TRUE. If no condition 
label is attached to the SET command, the arithmetic expression is 
evaluated at the first time step. Once the variable is set, its 
value remains fixed even if the condition label later becomes 
FALSE . 

2£22822_2Z _£*22_222_2£_222£*22£_§.££_222!5222 2£_2_££§£l_22!!!!n£2£^ 
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Exarapi e_ 5 ^l_ 62 lj 

SET comma nd 

APPLOAD, 12, ,0.0, Q1 ,0.0, 10.0,5 $ VARIABLE SCALE 
$ FACTOR 

LOAD, 12, 1.0$ SCALED LOAD 

SET,Q1 =(X(A, 6) + XU, 7) + , X ( A , 8 ) ) / 3 . 0 ? C2 $DEFINES 

$ Q1 

C2,X(D ,7) ' .GT. 0.5$ DEFINES C2 

In this example, a variable Q1 is set to the average 
value of the accelerations parallel to the global x 
axis at three grid points when the condition label 
C2 first becomes true; that is, when the 
displacement in the global x direction of grid point 
7 first exceeds one-half. Q1 is then used as a 
scale factor of an APPLOAD command which defines a 
load at grid point 5 over a maximum time interval of 
0 to 10 time units. Since all user defined 
variables are automatically initialized to zero, the 
scale factor is initially zero, and thus the applied 
load, does not exist at all before C2 becomes TRUE. 
Since Q1 is defined with a SET command, its value is 
determined when the condition statement associated 
with label C2 becomes TRUE. The value of Q1 will 
not change thereafter evenif the condition 
statement subsequently becomes FALSE. 
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To define a quantity whose value varies, the RESET command is 
used. It takes the form, 

RESET, Qn=FORTRAN arithmetic expression ? Cj 

where , 

n is an identification integer 

Cj is the j th optional condition label which is 
defined with a CLj command as discussed in 
section 5.15; C 2 _^canno^_cont a^n_a_us er defined 

condi t ion • 


The RESET command defines a variable Qn, where n is the 
identification integer. At every time step for which the 
condition statement attached to the RESET command is true, the 
value of the arithmetic expression of the RESET command is 
calculated and the variable reset to that value. If no condition 
statement is attached, the variable is reset to the value of the 
arithmetic expression at every time step. The variable may depend 
on current or past values of the structural response. In this way 
RESET may be used to define control forces on the structure 
including time delay. 


6 
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Exampl e_52_l62§.i User defined variable as load scale 

f ac^or_to_y^el d_veloc i t y_dependej 2 t 
co ntrol 

DEFINE;§6=1 .5707563 $ PI/2 

RESET, Q2= YP(V,8)$ DEFINES Q2 PROPORTIONAL TO 
$ VELOCITY ■ 

APPLOAD.3, ,@6,Q2,0.0, 10.0,8$ VARIABLE SCALE FACTOR 
LOAD , 3, 206 . 5$ LOAD VALUE OF 206.5 BEFORE SCALING 

In this example, Q2 is used as a scale factor for an 
applied load at grid point 8. From the RESET 
command, the scale factor, Q2, is equal to the 
velocity of grid point 8 parallel to the global y 
axis. The load amplitude of 206.5 in the third LOAD 
command, which is referenced by the APPLOAD command, 
will be multiplied by the variable scale factor. A 
velocity dependent control force is thus 
constructed. Unlike the previous example where the 
SET command was used, the RESET command permits the 
variable scale factor to continuously vary as the y- 
direction velocity of grid point 8 varies. To 
construct the same example with a time delay in the 
control force, see example 5.16-7. 

Ex ampl e_5^l623j User _vaMable_£n_a_condi^t^on_3 tat ement 

MEMBER , 7 , 7 , 8 , 2 , 2$ DEFINE MEMBERS SEVEN AND EIGHT 
RESET, Q3= YP(A,8,7) - YP(A,7,7) $DEFINES Q3 AS 
$RELATIVE TRANSVERSE ACCELERATION OF TWO GRID 
$P0INTS COMPRISING MEMBER 7 

RESET, Q4= YP(A,9,8) - YP(A,8,8) $DEFINES AS 
$RELATIVE TRANSVERSE ACCELERATION OF TWO GRID POINTS 
$ COMPRISING MEMBER 8 

CL22, Q3 .GE. Q4$ DEFINES C22 IN TERMS OF Q3 AND Qi( 

' INCMASS, 1 000 $MASS MATRIX UPDATE AT LEAST EVERY 1000 
$TIME STEPS 

INCMASS, 250? C22$ MASS MATRIX UPDATE MORE OFTEN IF 
$C22 TRUE 

In this example, beam members 7 and 8 are defined 
and the relative acceleration between their end 
points is set to user variables, Q3 and Q^. These 
variables then constitute a condition statement, 
CL22, which governs mass matrix updating. 
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Example 5.16-^: Combining user defined variables 


RESET, Q3= YP(A,8,7) - YP(A,7,7) $DEFINES Q3 AS 
$RELATIVE TRANSVERSE ACCELERATION OF TWO GRID 
$POINTS COMPRISING MEMBER 7 

RESET, Q4= YP(A,9,8) - YP(A,8,8) $DEFINES Q4 AS 
$RELATIVE 

$ TRANSVERSE ACCELERATION OF TWO GRID POINTS 
$COMPRISING MEMBER 8 

RESET, Q5=( Q3+Q4) /2 . $AVERAGE OF Q3 AND Q4 
APPLOAD, 1 8, ,Q5, , , 1 6$APPLIED LOAD WITH VARIABLE 
$MAGNITUDE AT GRID POINT 16 

LOAD, 1 8, 1 .0$THIS LOAD HAS A SCALE FACTOR OF Q5 

In this example, variables Q3 and Q4 of example 
5 . 16-3 are again used. They are averaged to produce 
Q5 which is employed as a scale factor for a load at 
grid point 1 6 . 
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E X ampl e_5 ^1_625 ^ E. i 

MEMBER, 1 i|, 9, 1 0, 1 , ALUMINUM$ MEMBER DEFINED 
RESET, Q9 = PHI ( 1 0) +PI/2$ DEFINE Q-VARIABLE NORMAL 
$T0 ROTATION OF MEMBER 1A AT GRID POINT 10 
APPL0AD,3,Q9, 3.0, , , 10$ LOAD IN DIRECTION Q9 AT GRID 
$ POINT 10 WITH MAGNITUDE 3 
LOAD , 3, 1 . $CONSTANT LOAD 

In this example, a Q-variable is reset at every 
time step to an angular value given by the rotation 
of grid point 10 plus ir/2. A constant in magnitude 
load is thereby oriented always normal to the beam 
members rigidly connected to grid point 10. 

NOTE 5.16-7: Up to 99 user defined variables may be employed. 
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5.16.3 


Time Delayed Variables. Via KEEP Command a nd Q-Histor.y 
Function ~ 


As already discussed, user defined variables, Qk , (k=1 
through 99), are defined in SET or RESET commands. They depend 
upon the analysis variables listed in Section 5.15 and 
Table i4.4-1. During the program execution, Qk is assumed to be 
evaluated at the current time step. The value of Qk at previous 
values of time are generally destroyed. However, the program can 
be instructed to retain values of certain user defined variables 
at a number of previous values of time. This is accomplished with 
the use of the KEEP command. Data for the KEEP command takes the 
f orm , 


KEEP : Qi , Q j , Qk , ns 

where, 

Qi , Qj , Qk are the ith, j th and kth user defined 

variables which may appear as arguments of QH 
functions 

ns is the number of previous time values, before the 
current time, for which the values of Qi , QJ or Qk 
are to be retained in memory; the values of the Q 
variables will be saved at all ns previous time steps 

Exampl e_5 . j_6-6 j T o_r e^a£n_£r e v^ous_val ues__of _a_us er 

~ defined varia ble 

RESET, Q12 = AX(22)$ Q-VARIABLE EQUAL TO AXIAL FORCE 
$ IN MEMBER 9 

KEEP, Q9,1$ Q9 RETAINED FOR ONE PREVIOUS TIME STEP 

In this example, a Q variable is defined to be the 
axial force in member 9. The value of Q9 is 
retained at one previous time step. Thus at any 
time during execution, the current value of Q9 and a 
retained value of Q9 one time step in the past is 
available to the user for employment in other . 
commands . 
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If the user wishes to employ a retained value of a user 
defined variable, rather than the current value, in a command 
permitting user defined variables as data entries, then a Q- 
History ( QH ) function is used, as 

QH(Qk,m) 

where, 

Qk is the k th user defined variable and appears in a 
KEEP command 

m is an integer specifying the number of previous time 
steps at which Qk is to evaluated; only the magnitude 
of m is used by the program and the magnitude must lie 
between 1 and the value of ns in the KEEP command 


QH functions may be used wherever a Qk variable is permitted. 
Exampl e_5j^i627] 

f act or_to_y^el d_velo city _depen dent 
co ntrol with time delay 


KEEP,Q2,5 $ PROGRAM TO RETAIN VALUES OF Q2 FOR 
$ PREVIOUS FIVE TIME STEPS 

RESET, Q2= YP(V,8)$ DEFINES Q2 PROPORTIONAL TO 
$ VELOCITY • 

APPL0AD.3, ,@6,QH(Q2,5) ,0.0, 10.0,8$ VARIABLE SCALE 
$ FACTOR 

$ EVALUATED FIVE TIME STEPS PREVIOUS TO CURRENT TIME 
LOAD , 3. 206 . 5$ LOAD VALUE OF 206.5 BEFORE SCALING 

This example is similar to that of 5.16-2 except 
that a time delay equal to five time steps is used 
in calculating the control force. 
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NOTE_5^2^"ll Only one KEEP command is permitted in a data setj;^ 

NOTE The absolute value of m in the argument of the QH 

function must lie between 1 and t'he value of ns in 
the KEEP, command. Since only the magnitude of m is 
us'ed by the' program i t is recommended; as a memory 
aid to, the user, that m be inserted as a negative 
integer in the' argument of the QH function. 

NOTE_5^1_623j_ Any or all of the defined Q-variables in the KEEP 
~ command may appear as arguments of QH functions. 
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Section 5.17 - Joint Lock-Up 
5.17.1 - General Discussion 


The user may desire to allow pin connected members to become 
rigidly connected when certain conditions are met. The grid 
points which make up the pin connection are denoted herein as the 
"joint" and the transition of the joint from pinned to rigid is 
denoted herein as "lock-up". Lock-up often occurs in analyzing 
the deployment of str uctures . and in some mechanisms. A lock-up 
occurence involves an impact and consequently involves transfer of 
momentum across a joint. Such momentum transfer depends upon the 
properties of the locking joint and the members which come into 
the joint. The LOCKUP command is used to provide the properties 
of the joint when it locks-up. If the user does not wish to model 
the dynamics of the locking joint, the ROTLOCK command may be used 
instead. This command conserves the momentum of the entire 
structure without detailing the structural behavior in the 
neighborhood of the locking joint. 


5 . 17.2 - LOCKUP Command 


The joint model prescribed by the opt io n al LOCKUP command 
consists of a nonlinear spring and viscous damper acting in 
parallel at the joint. The user provides the nonlinear spring and 
damper properties or uses certain default values. In addition, 
the user provides the condition under which lock-up is to occur in 
the form of a condition label.. 

When this condition is met, the nonlinear spring and damper 
are automatically activated. The spring and damper will quickly 
force the relative velocity between the grid points comprising the 
joint to zero. Once the relative velocity reaches zero, the 
spring and damper are automatically removed and replaced by a 
mul ti- degr ee- of - f r ee dom constraint which enforces rigidity of the 
joint. 

If no LOCKUP command is used, then the mul t i - degr ee- of - 
freedom constraint to enforce rigidity of the joint is applied as 
soon as the condition label defining lock-up is TRUE. 
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The nonlinear spring properties and the lock-up condition are 
provided by the user through the use of the op tional LOCKUP 
command. Data takes the form, ' 

LOG KU P,sl ,s2,c1 ,c2,n1 ,n2,n3i . ?C j 
where , 

si and s2 are coefficients describing the 

nonlinear lock-up spring constant, k, 
such that , 

k = si + ( s2) ( Qj ) ^ 

in which Qj is the current value of 
a user defined Variable; see below 
and section 5.15 

cl and c2 are coefficients describing a 

nonlinear joint damper whose damper 
value, c, is'given by, 

c = c1,+ (c2)(QJ)^ 

n 1 , n2 , n3 , n4 , . . . i 3 a string of grid point 

numbers comprising one or more joints; n1 
and n2 comprise the first joint, n3 and 
n4 comprise the second joint, etc. 

Cj is the j th condition label. As described 
in detail in section 5.15, the command to 
which the condition label is attached, in 
this case the LOCKUP command, is activated 
when the condition label is true. Thus, in 
this case, Cj expresses the condition which 
signifies joint lock-up. The condition 
label, Cj , is defined with the CL command as 
a logical expression written in FORTRAN. 

For use with the LOCKUP command, Cj mu st 
take the form, 

Qj .GT. 0.0 
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where Qj is an arithmetic FORTRAN expression 
defined by the j th RESET command as 
described in section 5.16. The 
value of Qj may depend on current 
values of problem variables, such as, 
displacements, velocities or 
accelerations . 


Exampl e_5^1_72l 

LOCKUP, 100. ,0.0, 10.0,0.0, 1 1/12 ? C5$ LOCK-UP AT GRID 
$ POINTS 11 AND 12 

C5,Q5 .GT. 0.0$ DEFINES C5 IN TERMS OF Q5 

DEFINE, §7 = 3. 1 'll 59265359$ PI DEFINED 

RESET, Q5=ABS(R0T( 1 1 )-R0T( 1 2) )-§7$ DEFINES Q5 

In this example, grid points 11 and 12 constitute a 
lockable Joint which locks-up when the fifth logical 
condition, C5, is true. C5 is true when Q5 exceeds 
zero; that is, when the magnitude of the relative 
rotation between grid points 11 and 12 exceeds ir . 
(Note that a DEFINE command has been used to define 
the value of ir . ) When C5' becomes true the LOCKUP 
command is activated and in this example, a linear 
rotational spring and viscous damper are 
automatically placed between grid points 11 and 12. 
The spring and damper are in addition to any that 
might be defined with the SPRING or DAMPER commands. 
From the LOCKUP command data, the spring stiffness 
is 100 and the damper coefficient is 10. When the 
relative rotational velocity at grid points 11 and 
12 becomes equal, the spring and damper will be 
automatically removed and replaced with a rigid 
mul t i - de gr ee- of - f r ee dom constraint which will make 
the joint rigid and no further relative rotation 
between the two nodes of the joint will occur. 
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N 0TE_5 ^1_7 - Momentum is always, conserved during- rock-up , but- 

energy is not generally conserved*. The f-.ollowing 
rules^ govern the use of the LOG KUP • command as they 
relate to energy conservation. 

If only flexible members (MTYP-E = 0.) come, into the joint 
and the user des->ires that energy be conserved, during lock-up, 
then the spring and damper properties of the LOCKUP command 
aresettozero. 

If any member involved in the lock-up of a joint is 
rigid (MTYPE=3) then energy- cannot be conserved during lock-up 
and nonzero spring and/or damper properties should be used to 
remove energy . 

sl=i000 
s2=0 
cl =1 00 
c2=0 

NOT E_bj_1722 At the present time the LOCKUP command is not 
available. 
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5j_i_7^3_2_5,2Ik2.9,5_9.2!!l!B®E5i l2_k£2!Sl2£_?.2k2ki2i_?lf!!lk®22_k£i28 

Co ns ervation of Momentum, Buk_^ik!l22k 

Deta il Analysis In the Neighborhood of 
the Locking~Joint 


The ROTLOCK command takes the form, 

ROTLOCK: igri d1 , igrid2 ? Cj 
where, 

igridl and igrid2are the grid point numbers involved in 

the lockup which prior to lock-up have 
no rotational constraint defined 
between them 

t h. 

Cj is the J condition label; £_£ 0 ! 2 d^t^on_i abel_j_s 
mandatory for this command 
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Section 5.18 - Output Control 
5 ^ 1_8 . 1 "General Discu ssion 

The user is' given considerable versatility in the form and 
quantity of output desired. The following results from the 
program may be output in printed or plotted form: 

displacements, velocities, accelerations 

internal beam member forces; i.e., axial forces, shear 
forces, bending moments at both ends of the member. 

reactions 

mass matrix 

tangent stiffness matrix (if STIFF command has been used) 
frequencies and mode shapes (if FREQ command has been used) 
user defined variables ( Q- var i abl es ) 
status of condition labels 

The user may print displacements, velocities and 
accelerations or internal member forces at selected time intervals 
or plot displacements, velocities, accelerations or internal 
member forces in the time or frequency domain. The user may also 
choose to print the mass or tangent stiffness matrix or to plot a 
trace (non-zero locations of terms) of these matrices at selected 
time intervals. Frequencies and mode shapes which have been 
selected for calculation in the FREQ command are automatically 
printed. The user may plot the translating, rotating and 
deforming total geometry of the structure at selected time 
intervals. In addition, user defined variables and the status of 
condition labels are automatically printed and may be plotted if 
desired. 
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5.18. 2 - Sign Convention.s For Interpretation of- Grid; Point Motions 

As shown in figure. 5... 6-V, t;ranslation mot.ipns, and are 

positive parallel to the global X and Y; a.xe.s respectively. 
Rotational grid point motion, ib positive if rot.atio.h is in a 
direction from the positive glb,b.al’X axis to the. positive global Y 
axis. 


Also, as shown in figure 5.18-2, member internal axial 
forces, and Pg are positive parallel to the local V vector. 

The local V vector runs from end A t,o end B of the member. 
(MEMBER commands determine ends A and B. Members run from end A 
to end B.) Shear forces, S. and S„ are positive parallel to the 

A ‘ D ' 

local transverse Y' axis and bending moments M. and M^ are 

A D 

positive in the positive rotational direction. 
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Figure 5il8-1. Sign Conventions For Internal Member Forces 
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Sect i j Prin t 1 ng 

Printing is accomplished by using the opt i o,nal PRINT command 
in conjunction with the' GRDPRT -and/or MEMPRT ‘commands . Data for 
the PRINT command takes the form', 

PRINT:k.1,i1,k2,LG, k3,k4 ? Cj 

where , 

k1 is an integer specifying the printing of 

displacements, velocities and accelerations 
every k1 time steps; if k1 is zero or left 
blank, there is no request for printing of 
d i spl a.cements vel oc i t i es and accelerations 

i1 is an integer specifying that printing of 

displacements, velocities and accelerations is 
to be in the local axis system which rotates 
with member i1; if il is zero or left blank, the 
global system is used 

k2 is an integer specifying the printing of 

internal forces (axial force, shear force and 
bending moment at each end of the member) every 
k2 time steps; if k2 is zero or left blank, no 
request is made for printing of internal forces. 

LG is a letter (either L for Local or G for Global) 
specifying that printing of internal member 
forces are to be in the local or global axis 
system which rotates with each member; if LG is 
left blank, the Global system is used 

k3 is an integer whose absolute value specifies the 
printing of the mass matrix every lk3| time 
steps; if k3 is zero or left blank, no request 
is made for printing the mass matrix; if k3 is 
positive, the mass matrix is printed; and if k3 
is negative traces of the mass matrix are - 
printed 
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k is an integer whose absolute value specifies the 
printing of the stiffness matrix (including 
incremental stiffness if requested in STIFF 
command) every jk4| time steps; if k4 is zero or 
left blank, no request is made for printing the 
stiffness matrix; If k^t is positive the 
stiffness matrix is printed; and if k4 is 
negative traces of the stiffness martrix are 
printed 

Cj is the j th o pt i onal condition label defined 

with a Cj command; see section 5.15 for further 
details 

NOTE 5.18- 1; In addition to the even increment printing specified 
in the PRINT command, printing will always occur at 
the start time in the TIMSPAN command and when 
certain condition statements become TRUE. 
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Example 5 . 1 87 1 


TIMSPAN:: 1 0 . , 2000. 

P R I NT : Y, 0 0 , 5 , 10 0 , 5 , 1 0 0 , 0 

In this example, the mass matrix, displacements, 
velocities, accelerations and all internal member 
forces (axial forces, shear forces and bending 
moments a.t both ends of designated members) are to 
be printed at the starting time of 10 as given in 
the TIMSPAN command and every 100 time steps. 

R es'ul ts wi 1 1 'be given in the time varying local 
coordinate system which rotates with member number 
5 . The tangent stiffness matrix is hot printed. 

Note that the first entry following’ the command name 
instructs printing of displacements, velocities and 
accelerations simultaneously at designated grid 
points, while the. fourth entry follow'ing the command 
name instructs printing, of all internal forces 
simultaneously in the designated members. 

Designated grid point numbers and member numbers for. 
printing are provided using the commands GRDPRT and 
MEMPRT described next. 
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Unless otherwise specified, printed output of displacement, 
velocity and accelerations will occur for all grid points and 
printed output of forces will occur for all members. The commands 
GRDPRT and MEMPRT limit printed output to subsets of grid points 
and members . 

The GRDPRT command takes the form, 

GRDPRT, igridi ,igrid2 

where, 

igridi, igrid2 constitutes a list of grid points by 

specifying actual grid point numbers 
and/or predefined list(s) of grid 
point numbers (see page 5.2-2) 

Example 5.18-2: Use of GRDPRT Command 

LIST,To, 1 , 2, 3.5T1 565"“$ DEFINE A LIST OF 
$ FOR ■ LATER USE 

GRDPRT, LI 0, 11 $ SUBSET OF GRID POINTS 
$ FOR PRINTING 

In this example, a list of grid points was defined 
somewhere in the data set. This list was denoted 
as "10" and contains grid points 1, 2, 3, and 5, 10 

and 15. The GRDPRT. command creates a subset of 
grid points for printing by referencing the 
previously defined list as "L10" and adding to it 
grid point H . 
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The MEMPRT command takes the form, 


M E M P R T , i m e m 1 , i m e m 2 , 

constitutes a list of m'eraber numbers 
by speci'fylng. actual member numbers 
and/or predefined listCs) of member 
numbers (see page 5.2-2) 

Example 5.18-3': Use of MEMPRT commah'd 

LIST , 8 , 20TH0B5 $ DEFINE A LIST OF 
$ FOR LATER USE 

MEMPRT, L8 $ SUBSET OF MEMBER NUMBERS 
$ FOR PRINTING 

In this example, a list of member numbers was 
defined somewhere in the data set. This list was 
denoted as "8" and contains 20, 25-, 30, 35 and 40. 

The MEMPRT command creates a subset of member 
numbers for printing by referencing- the previously 
defined list as "L8" . 


where , 

i mem 1 , i mem2', 
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Reactions are printed by using the REACT command. Data takes 
the form, 


REACT, igrid,idir,m 
where , 

igrid is a grid point number , 

idir is a direction; idir=1 for x-direction, idir=2 
for y-direction and i.dir = 3 for rotation 

m is the member number of the local coordinate 

system in which the reaction is to be calculated; 
ra equal to zero or left blank indicates the 
global axis system 

Reactions associated with grid points and load 
directions specified in a REACT command which are 
not currently dependent degrees-of-freedom due to an 
active SDFC or MDFC command, or rigid member 
constraint, or for some other reason, are calculated 
as zero. See also Notes 5 . 13~1 through 5 . 13 ~ 3 . 

A condition label cannot be used in conjunction with a REACT 
command . 
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5.18.M; Plotting 


If resultsare to be. plot ted, the program must be instructed 
to create an plot output file. Any or all of the .data on the 
plot output file can be accessed through the use of plotting 
commands. The plotting commands provide instructions to the post- 
processor and may either, be, placed in the data set and thereby 
executed immediately f ollowiri-g the completion of. all LATDYN 
calculations or used at a later time during a poat-'procesSi ng 
session. Thus any data that has been placed on the plot output 
file and appropriately saved on disk or magnetic tape, can later 
be accessed and the results plotted. (For further discussion see 
pos t- pr ocessi ng in Appendix A.) 
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5.18.M.1: Creating a plot output file 


Two commands are required to instruct creation of the plot 
output file. They are, 

PLOTINC - specifies the density of data points used to 
generate plots 

PLOTLIMIT - specifies the time duration during which data is 
placed on the plot output file 

These commands take the form, 

PLOTINC , r?Cj 

'Where , 

r is an integer specifying that all plots are to be 
generated from data calculated every r th time step 

Cj is the j th optional condition label; see section 
5.15 for further details 

N0TE_5 j. Multiple PLOTINC commands may be used, but they have 
a condition-order priority as discussed in NOTE 
5 . 15-1 of Section 5.15. Thus, the first appearance 
in the data set of a’ PLOTINC command cannot have an 
attached condition label. See also Section . 5 . 

PLOTLIMIT, t1 , t2 

where. 


1 1 

is the 
output 

1 ower 
file 

time 

limit 

for 

all 

data on 

the 

plot 

t2 

is the 
output 

upper 

file 

time 

limit 

for 

all 

data on 

the 

plot 


NOTE _5 ^ PLOTLIMIT command is not used, the time limits 
of the TIMSPAN command are employed. 
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Once the LATDYN e xe c u t i o n. : Qf results is complete • the post- 
processor can be used- to ^piot ^ariy quantities placed on the plot 
output file. The us er i ns t r.uc.ts . t he post-processor using commands 
which fall into one of six categories, 

(1) generation of time history plots 

(2) generation of user defined variables ( Q- var i abl es ) 
plotted aganist one another 

(3) generation of pictoral view of structure during motion 

and deformation , 

(^) generation of frequency domain plots 

(5) zooming for macro and micro viewing of time and frequency 
domain plots 

(6) user assistance during pos t- pr ocessor interactive session 

The pos t- process or commands may appear in the data set and be 
executed directly after LATDYN computations are finished, or may 
be used at a later time in an interactive session with the post- 
processor . 

NOTE 5 . 18 -5 : The use of a command in category (1) generates a set 
of time histories of response on one plot; that is, 
the time' variation of all designated quantities in 
the command will appear on one set of axes. 
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Post-processing Commands - Category 1 


The following commands generate time history plots of data 
that reside on the output file. 

PLOT XLOC , n 1 , n2 , n3 , . . . to plot time histories of global X 

position at prescribed grid points 

PLOT YLOC , n1 , n2 , n3 , . . . to plot time histories of global Y 

position at prescribed grid points 

PLOTROT , ml , m2 , m3 I . . . to plot time histories of angular 

position of prescribed members; 
positive angle measures in direction 
from positive global x to positive 
global Y axis 


where , 

n1,n2,n3)... represents a set of grid points numbers 
m1,m2,m3»... represents a set of member numbers 
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Po st-processing Commands - Category -1 (c ontin ue d) 

The following commands g,en,erate time his.tory pLots. of data 
that resides on the oatpu,t file. 


PLOT MA G , p , n 1,, n2 , n3 ». ...... 


PL0TX,p,n1 ,n2, n3i . . . 


P LOT Y,p,n1 ,n2,n3. . . . 


PLOTPHI , p , n 1 , n2, n3 , . . . 


to plot time, hi stor i es. 0 ;f magnitude 
of mot.ion quantities (displacement, 
p=:D; velocity, p=V; or acceleration, 
p = A) at a set of grid, points 

to. plot time histories of motion 
q.uantities (displacement, p = D; 
velocity, p=V; or acceleration, p=A) 
parallel to global X axis at a set of 
grid points 

to plot time histories of motion 
quantities (displacement, p=D; 
velocity, p=V; or acceleration, p=A) 
parallel to global Y axis at a set of 
grid points 

to plot time histories of rotational 
quantities (displacement, P=D; 
velocity, P=V; or acceleration, p=A) 
at a set of grid points; positive 
angular rotation is in a direction 
from from positive global X axis to 
positive global Y axis 


where , 

p is a parameter designated as , 

D - for displacement 
V - for velocity 
A - for acceleration 

1 represents the member number of a local 

coordinate system; if zero or blank the 
global axis system is assumed 

nl,n2,n3,... represents a set of grid points numbers 
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Post-processing Commands - Category 1 (continued) 


The following commands generate time history plots of data 
that resides on the output file. 

PLOTXP , 1 , p , n1 , n2 , n3 . . . . to plot time histories of motion 

quantities (displacement, p=Dj 
velocity, p=V; or acceleration, p=A) 
parallel to a local X' axis of member 
1 at a set of grid points 

PLOT YP , 1 , p , n 1 , n2 , n3 , . . . to plot time histories of motion 

quantities (displacement, p=D; 
velocity, p=Vj or acceleration, p=A) 
parallel to a local Y' axis of member 
latasetofgridpoints 

P LOT AN GLE , p , n 1 , n2 , n3 » • • . to plot time histories of angular 

orientationof a motion quantities 
(displacement, p=D; velocity, p=V; or 
acceleration, p=A) at a set of grid 
points; positive angular orientation 
is measured in a direction from 
positive global x axis to positive 
global y axis 

PLOTSTRAIN , ml , m2 , m3 I • . • to plot time histories of axial 

strain in a set of members 

PLOTFEXT , ma , Idl , ld2 , to plot external force magnitude if 

ma is M or orientation angle if ma is 
A for APPLOAD sets designated by 
load 1 , load2 , . . . 

where , 

p is a parameter designated as , 

p - for displacement 
V - for velocity 
A - for acceleration 

n1,n2,n3,... represents a set of grid points numbers 

ml, m2, m3,... represents a set of member numbers 

1 represents the member number of a local 

coordinate system; if zero or blank the 
global axis system is assumed 
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P OS t- pr 0 c essi ng Commands - Category 1 (continued) 


The following cpmmands generate time history plots of data 
that reside on the output file. 

PLOTC , c 1 , c 2 , c 3 > • . . to pl'ot TRUE or FALSE time histories 

6f a set of condition labels 


P LOT Q , q 1 , q2 , q 3 . . . . to plot time histories of a set of 

user defined variables 


PLOT AX, ml ,m2,m3. • • . 


PLOTSHA.ml ,m2,m3, . . . 


to plot time histories. of axial 
forces in a set of members (tension 
is positive) 

to plot time histories of shear force 
at end A in a set of members 


PL0TSHB,m1 ,m2,m3i . . . 


to plot time histories of shear force 
at end B in a set of members 


PLOTMOMA.ml ,m2,m3, . . • 


to plot time histories of bending 
moments at end A in a set of members 


PLOTMOMB , m 1 , m2 , m3 I . . . to plot time histories of 

moments at end B in a set 


where , 

p is a parameter designated as , 


bending 
of members 


D - for displacement 
y - for velocity 
A - for acceleration 


m1,m2,m3.*-- represents a set of member numbers 

c1,c2,c3.**» represents a list of condition label 
numbers 


q1.q2,q3i*.. represents a list of user defined variable 
numbers 
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The following command creates a plot of user defined 
variables ( Q- var i abl es ) where one user defined variable is plotted 
aganist another. Data takes the form, 

PL0TQVQ.Q1 ,Q2,Q3, . • . 

where user defined variables Q2, Q3, etc. are plotted aganist Q1. 
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Po3 t- proc ess£ng_C ominands_ 2 _C at e gbr y _3 

The following command generates pictoral view's' of the structure 

SNAPSHOT , ts , te , k ,F, M , N to plot the di spl aceraentV def or mat i on 

of the structure every k th time step 
from time "ts" to "te" with drawn 
forces if F. appears; member numbers 
if M appears'; grid point numbers if N 
appears; drawn axes if A appears; and 
overlaid snapshots if 0 appears 

When the SNAPSHOT, command i's used without the presence of ts 
and te, d i spl acem en t / def or mat i on plots of the structure are 
generated at the following times: 

starting time as prescribed in the TIMSPAN command, (thus 
producing an initial deformation plot), 

at time tl as prescribed in the PLOTLIMIT command 

and every k th time step after time tl. 
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Po3t-proce33 i ng Commands - Catego ry ^ 

Results on the output file may be reduced to produce 
frequency domain plots. . Power spectral density and power spectral 
level plots may be generated. Data takes the form, 

PSD - for power spectral density (units are power/Hertz) 

PSL - for power spectral level (units are power/ spectral 
resolution bandwidth) 

The PSD and PSL commands apply to the time history command which 
immediately proceeded them. 
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Post-processing Commands - Cate go ry 5 

The TIMEWINDOW command is used to set the axes pf the. time history 
plots. Effectively the user can zoom-in on results pf ' i nt eres t . 

The command , ta k es the form, 

TIMEWINDOW, TIME1 ,TIME2 

where TIME1 and TIME2 provide upper and lower bounds for the 
time scale on time history plots and plots of Q-variables versus 
one another. The time window must be within thie bounds of the 
PLOTLIMIT command. If TIME1 and TIME2 are left blank, the bounds 
of the PLOTLIMIT command are usedi' ' ' 

The position in the data sequence determines which time 
history plots are affected by a particular TIMEWINDOW command as 
each command remains active. until anpther appears inthe data set. 


The FREQWINDOW command is used to set the axes for plots in the 
frequency domain. Effectively the user can zoom-in on results of 
interest. 


The command takes the form, 

FREQWINDOW, FREQ 1 , FREQ 2 

where FREQ1 and FREQ2 provide upper and lower bounds for the 
frequency scale on plots in the frequency domain. 

NOTE 5.18 -6 ; The position in the data sequence determines which 
frequency domain plots are affected by a particular 
FREQWINDOW command as each command remains active 
until another appears in the data set. PSD and PSL 
commands must follow a FREQWINDOW command. 
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During an interactive post-processing session the user can 
get assistance from the post-processor program via the commands 
HELP and INFO. 

The HELP command provides on-line quick reference to a 
summary of pos t- processor commands and their function. 

The INFO command provides on-line quick reference to 
important characteristics of the data on the output file and 
current status of the post-processor. This includes the job title 
and subtitle, time limits, current time window, number of time 
steps, grid point and member numbers, Q-variable identification 
numbers and condition label numbers. The user will find this 
command very useful, since what is on a particular data file may 
be forgotten . 
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P lot t£ng_Exampl es 

iot_D l^£P^acem^n t s_^n_G^opal_X 2 
Directi on 

PL0TINC,100$ SET PLOTTING DENSITX 
PLdTLIMIT/p;0,20.d$ SET UPPEN AND LOWER TIME 
$’ LIMITS f6r plots 

PLOTX.D ,2, r, 6, 8$ PLOT DISPLACEMENTS AT GRID 
$PTS. 2; 4,6,8 ■ 

In this example, the plotting specifications, 
as supplied via the PLOTINC and PLOTLIMIT 
commands, are the generation of time history 
plots from time 0 to 20 using results 
calculated every 100 time steps. The PLOTX 
command indicates that quantities to be 
plotted are components of vectors parallel to 
the global x axis. The plot parameter "p" 
being set to D arid the string of integers 
indicate that displacement time histories are 
desired at grid points 2, 4, 6 and 8. Since 

only one PLOTX command is used, a single plot 
is created containing the four time histores 
of the four designated grid points. Scaling 
of the plot will be performed automatically on 
the basis of the results for all four 
histories . 
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Example 5.18-4; Plot Acceleration in Local System 
PLOTINC.IOOf SET PLOTTING DENSITY 
PLOTLIMIT, 0.0, 20.0$ SET UPPER AND LOWER TIME 
$LIMITS FOR PLOTS 

PLOTYP, 2, A, 2, 4$ ACCEL PLOTS TRANVERSE TO 
$MEMBER 2 AT $ POINTS 2,4 

In this example, the PLOTYP command indicates 
quantities to be plotted are components of 
vectors parallel to a local y' axis which is 
transverse to a structural member. The first 
integer value indicates that the local axis is 
associated with structural member 2. The 
parameter "p" being set to "A" and the string 
of integers which follow, indicate that time 
histories of acceleration are desired at grid 
points 2 and 4. Plot specifications are the 
same as in the previous example. 
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Example 5. l8-b;. P-'lot Bending Mom ent 

RLO.TINC ,.1 00$. SET: PLOTTING DENSITY; 
RL0TLIML.T.,0.b,20.. 0$ SET UPPER A-NO LOWER TIME 
LIMITS FOR $ PLOTS 

RLOT'MOM.1 ,.V2,.20$ INTERNAL BENDIN-G. MOMENT 

In this- example, the. PL0TM0M1 command 
indicates quanfcites to be plotted are. bending 
moments- at end 1 of designated members. The 
designated members- are prescribed by the 
string of values, namely, 12 and 20. The plot 
specifications- are the same as those of 
Example 6.18-3. 

E X a m p 1 e _5 ^18 - 6 j Use of Condition label to ou tp ut 

additional plots during certa in 
phases of the calculati ons 
PL0TINC,100$ SET PLOTTING DENSITY 
PL0TINC,25? C10$ SET ALTERNATE PLOTTING 
$DENSITY 

C10,& $ CONTINUATION LINE FOLLOWS 
( ABS(SH1 ( 1 2) ). .GT. 1.E4 .OR. ABS ( SH 1 ( 20 )') & 

.GT. I.Ell) $ PREVIOUS LINE SETS CONDITION 
$LABEL CIO 

PLOTLIMIT,0.0, 20;..0$ SET UPPER AND LOWER TIME 
$LIMITS. FOR PLOTS 

PLOTSHA, 1 2, 20$ SHEAR FORCES IN MEMBERS 12 AND 
$20 

In this example, the PLOTSHA command indicates 
that time histories of shear forces at end 1 
of designated members are to be plotted. 
Designated members are 12 and 20. Plotting 
normally uses results generated every 100 time 
steps, but when condition CIO is TRUE, 
(magnitude of shear force in either member 
exceeds 10,000 units), plotting uses results 
generated every 25 time steps. 
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Example 5.18-7: Plot Magnitude of Velocity Vector 
PLOTINC.IOO^ SET PLOTTING DENSITY 
PLOTLIMIT, 0.0, 20.0$ SET UPPER AND LOWER TIME 
$LIMITS FOR PLOTS 

PLOTMAG, V, 5, 1 0, 1 5, 20$ MAGNITUDE OF VELOCITY 
$VECTOR AT GRID POINTS 5, 10, 15 AND 20 

In this example, the PLOTMAG command indicates 
that time histories associated with the. 
magnitude of a designated quantity is plotted. 
The "V" entry indicates the designated 
quantity is velocity. The following string of 
integers are the grid points at which the 
magnitudes of the velocity vectors are to be 
plotted . 

Ex a m p 1 e_5 ^i 8 2 §.J 

Plots 

TIMSPAN, 0. 0, 50 . 0$ STARTING AND STOPPING TIME 
PLOTLIMIT , 6 . 0 , 1 5 .0$ PLOT FROM TIME OF 6 TO 15 
SNAPSHOT ,, 1 000$ DEFORMATION PLOTS EVERY 1000 
$TIME STEPS 

In this example, structural deformation plots 
are indicated at time equal to 0, (the initial 
time in the TIMSPAN command), time equal to 6, 
(the initial plot time in the PLOTLIMIT 
command), and every 1000 time steps thereafter 
up to time equal to 15, (the final time in the 
PLOTLIMIT command). 
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SECTION 6 - USER GUIDELINES 

User guidelines for selection of the following data values 
are provided in this section: 

Time Step In TIMSTEP Command 

Increment For Mass Matrix Update In INCMASS Command 


Section 6.1 - Time Step In TIMSTEP Command 

Selection of the time step is very important to successful 
program execution. When using any explicit time integration 
operator (INTGTYP set to 1 through 4), the time step must be 
sufficiently small so as to avoid numerical instability. In 
linear problems it can be shown that the time step must be less 
than twice the reciprocal of the highest system . frequency . 

However, most of the problems that LATDYN would normally be used 
for are nonlinear. Nevertheless, the linear criterion for 
numerical stability usually provides a good guideline when the 
highest instantaneous system frequency is used in the criterion. 
This frequency can be derived from the current tangent stiffness 
and mass matrices. However, this information is not generally 
available to the user prior to program execution. It is suggested 
then that the following procedure be used when a satisfactory 
selection of the initial time step cannot be made prior to program 
execution. 
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The suggested procedure for determining a satisfactory 
initial time step i.s a .trial and error approach based on the 
general behavior of explicit, time integration al g^'r ithms . In 
general, for these algorithms, the larg'est time step^ for which the 
solution is numerically- siable is also, the most' efficient and 
provides results, which -are, nearly as accurate a's any smaller time 
step. Thus, the pr o;c;ed ur'e is to start' with a numerically stable 
time step and increase that time step until numerical Instability 
is found and then, to the degree desired, iterate', to the largest 
stable time step. The suggested steps in this procedure follow. 


(1) Find the highest system frequency which usually will be 
associated with a localized vibration mode. To do this, 
locate the shortest flexible finite elementof the 
model. Calculate its bending and extensional vibration 
frequencies using the formulas, 


0) 


( k L ) c/ sL 


and m = c/ ( 2L ) 


and select the higher value where, 

kL is the char act er i s t i c beam value and it is 

suggested to use the value for a fully clamped 
beam, namely, 4.730 

c i s the elastic wave speed in the rnembtjr; that is, 
/(E/p) 

which is about 200,000 inches /second for aluminum 
L is t h o e 1 e m e n t length 

s is the slenderness ratio; that is, L/[/( I/A ) ] 
RLsYoung’smodulus 

P is the mass density per unit volume 
I is the element cross sectional moment of inertia 


A is the element cross sectional area 
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(2) Perform a series of small trial runs to establish a more 
efficient time step. This may be done by executing the 
program, (preferably interactively), using 2/o) as a 

m 3 X 

first choice for the time step and using a terminating 
time in the TIMSPAN command equivalent to about 100 time 
steps. If results grow unbounded, reduce the time step. 
When a stable time step has been found, double its value 
and rerun. Continue this process until results grow 
unbounded . 

(3) Use the last stable time step or if desired iterate to 
find a better selection. However, it is not usually 
advantageous to expend much effort in iteration since 
the problem being solved is usually nonlinear and could 
possibly become unstable later in the analysis. 

(^) Execute program with selected time step. If INTGTYP is 
3 or 4, time step will automatically be adjusted and 
numerical solution will likely remain stable. In any 
case, the time step may be changed during execution via 
user options in the input data. The following partial 
runstream serves as an example in accomplishing this: 
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Examp^e_ 62 j_j R.uns tream For -User.. Contr^ol of ~T Ime . S t ep : Bas e 'd o n 

Pr og ram' Calcu-lati on'~ of Hig he st Sy'stem^Frequeh cy 

STIFF.,ap0,0 , UP.pA.TE TANGe'nt' STIFFNESS mATRTX 
INCMA SS . 1 b b 0$ ,U P DAT E-' ;M A S S M A T R I X . 

FREQ, 1 0.0’0,.1 , 0$' calculation' AND PRINTOUT OF HIGHEST 
$FREQUENCY‘. , . .. 

RESET , Q.1 =0 .95*( 2/, OMEGA) , 

TIMSTEF’ ,( seiecte'd"’ value from steps 1 to' 3)$' INITIAL TIME 
$ STEP 

C5, T ,GE. 1 0 00 * ( s el ect ed time step from steps 1 to 3) 
TIMSTEP , Qi ?C5$' new' TIME STEP EQUALS Q1 

In this sample par t i al . runs tr earn', the highest system 
frequency is calculated by the program every 1000 time 
steps using updated stiffness and mass matrices. A user 
defined variable, in this case Q1 , is defined as twice 
the reciprocal of the highest frequency and an 
additional 5$ reduction is arbitrarily added. The 
variable Q1 is used as the new time step after the first 
1000 time steps. A condition label, in this case C5, is 
defined to activate the new time step after the initial 
time step has been used 1000 times. 

NOTE 6.1 -1 : It is not advisable to .calculate the highest 

frequency too often during a run as this' can increase 
computer solution time considerably. 


Time Step Selection Example 
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Selection of member type, MTYPE in the MATPROP command, 
influences time step selection. If a member is very stiff and is 
treated as a flexible member, MTYPE=0 in the MATPROP command, with 
a high stiffness or cross sectional property, the highest 
frequency of the system will be greatly increased, thus decreasing 
the necessary time step for numerical stability when INTGTYP is 1 
through It would then be desirable to use a rigid member, 

MTYPE=3 in the MATPROP command, as this would eliminate the high 
frequency from the problem. However, the use of rigid members 
requires updating of the mass matrix as discussed in the next 
section. Consequently, the user must trade-off the small time 
step of stiff members against the mass matrix updating. Each 
results in incresed computer solution time. In general, though, 
reduced time step has a greater adverse affect on solution time, 
than mass matrix updating. 
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Section 6.2 - Incremental Updating of Mass Matrix 

LATDYN uses a consistent mass formulation in its solution 
procedure.' Because of the way in which the finite elements are 
modeled, the consistent ma'ss m at rix will change very litt le during 
program execution even if large angular changes between members 
occur, pro vided all members are fully fle xi ble, MTYPE=0 in the 
MATPROP command. Thus when all members are flexible, few if any 
updates to the mass matrix are required for large angular motions 
of members. Other considerations, such as updating of the mass 
matrix for frequency calculation as used in Section 5.1 will then 
govern the use of INCMASS. 

If rigid members, MTYPE=3 in the MATPROP command, are used, 
the mass matrix will need to be updated. The updating increment 
depends upon the angular rotation rate of the rigid members. That 
is why the INCMASS command permits an incremental angle entry 
which instructs mass matrix updates when any rigid member rotates 
through the specified incremental angle since the previous update, 
independent of what caused the previous update. The default value 
of five degrees for this angle is based on studies performed on 
mass matrix updating. The sample problem of Section 6.1 
ill ustrates , the effect of these updates on solution accuracy. As 
described in Section 5.3i mass matrix updating will also occur at 
even increments of time steps in accordance with the first 
numerical entry in the INCMASS command. 

Mass matrix updating must also occur when constraints become 
activated or deactivated. This is handled automatically by the 
program. The user has no control over such updates; they are in 
addition to to the periodic update specified in the INCMASS 
command. 

Independent of what caused a mass matrix update, a message 
about its occurrence will appear in the output file for user 
information. 
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SECTION 7 - SAMPLE INPUT AND OUTPUT 


Presented in this section of the User's Manual are four 
representative sample problems which illustrate the use of LATDYN. 
Included with each problem are the actual input data cases along 
with a brief description of the commands used. In addition, 
sample output is presented in the form of printed output from 
LATDYN and various plots generated through the pos tpr ocessor . 

Problem 7.1 - Rotating Rigid Beam 


This sample problem consists of a single rigid beam having 
concentrated tip masses and given an initial rotational velocity 
as shown in Figure 7.1-1. The main purpose of this example is to 
illustrate the use of all the required commands, some of the 
optional commands, comments, continuation lines and the free field 
input data format. 

Input Data 


The input data for this problem is shown in Figure 7.1-2. 

This case uses in addition to the required commands, 6 optional 
commands, each of which is discussed below (Note that although the 
SUBTITLE command is optional, it is considered essential and 
discussed along with the other required commands). 

Listed below are the required commands necessary for 
successful execution of LATDYN, 

TITLE 

SUBTITLE (See note above) 

PORDER 

TIMSPAN 

TIMSTEP 

MATPROP 

GRID 

MEMBER 

Note that although the program will execute without the PRINT or 
PLOTINC commands, it is essential to use at least one of them in 
order to get any meaningful output from the program. 

The TITLE and SUBTITLE commands each consist of one line of 
alphanumeric data which is printed out as part of the output. The 
user is encouraged to keep accurate descriptions through these 
commands for ease in identifying individual data cases. 

The PORDER command specifies the maximum number of finite 
elements and grid points in the model. To maximize program 
efficiency it is recommended that the values specified be the same 
as the actual values in the problem. In this example, one beam 
element and two nodes are specified. 
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The TIMSPAN command gives the starting and terminating times 
for the analysis. For this problem, the start and stop times are 
0.0 and 1.0 respectively.’ Note that through the 'use of a STOP 
command (used in Problem 7.2) the user’ can terminate the analysis 
prior to the specified terminating time if a previously specified 
condition becomes true. 
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Y 

4 ^ 



Vj = 18.85 in/sec 
RIGID ELEMENT 


^ODE 1 
• rad/ sec 


NODE rad/sec 





Vg = -18.85 in/sec 


Figure 7.1-1 Rotating Rigid Beam With Initial Velocities 
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:$ 

STHIS sample: PRLDE^tLEM IS- A RIGID BEAM WITH- ATTACHED. TIP MASSES WHICH 
■THAS AN INITIAL ROTATIONAL VELOCITY SUCH THAT THE. BEAM WILL ROTATE 
•^EXACTLY PI RADIANS IN ONE SECOND. 
t- 

TTITLE: DATA CASE TITLE, U UP TO-. 8.0 CHARACTERS) 

TITLE: T9S - FIRST TEiST: CASE FOR LATDYN- USER'S GUIDE 

t- 

^SUBTITLE: DATA CASE- SUBTITLE (UP TO SO CHARACTERS) 

SUBTITLE: ONE RIGID.' ELEMENT WITH; INITIAL VELOCITY 
T 

^DEFINE: <aN=CONSTANT 

DEFINE: @1=. 0002588 ^DEFINE MATERIAL DENSITY 

DEFINE: @2=3.141592654 T DEFINE PI 

@3= l.OE-3 ^DEFINE VALUE OF CONCETTRATED MASS 
•T 

REORDER: NELM, NODES 

PORDER: 1, 2 tONE BEAM ELEMENT AND TWO NODES ARE DEFINED 

1- 

$T I MSP AN: START, STOP 

TIMSPAN: 0.0, 1.0 f ANALYSIS TO RUN FOR 1 SECOND 

t- 

$TIMSTEP:DELTA T 7C0NDITI0N LABEL 

TIMSTEP: l.OE-3 INTEGRATION TIME STEP IS .001 SECONDS 

t- 

•TMATPROP: MATNAME,MTYPE, E , A , I , RHO ? CONDITION LABEL 
MATPROP: ALUMINUM, 3, 10.5E06, @2, .78540, @1 T -USING RIGID MEMBER 

•T WITH ALUMINUM MATERIAL PROPERTIES 
■t 

TGRID:NSr,# NODES, XI , Y1 ,X2 , Y2 ,INC 

GRID: 1, 2, 0.0, 0.0, 12.0, 

0.0, 1 TTHIS IS A CONTINUATION OF THE LINE ABOVE — TWO GRID POINTS 

T HAVE BEEN DEFINED 
T 

TMEMBER: IMEMB ,N1 ,N2 , NE ,MAT 

MEMBER: ,1, 2, 1, ALUMINUM TONE MEMBER FROM NODES 1 TO 2 USING 

TTHE PF^OPERTIES OF "ALUMINUM" FROM THE MATPROP COMMAND 
$ 

T 

■TPRINT:DVA INC, COORD.INTF INC, COORD, MM INC,SM INC ‘^’CONDITION LABEL 
PRINT: 100, O, 100, L, 0, O ^^PRINT DISPLACEMENTS/ 

^VELOCITIES/ ACCELERATIONS EVERY 100 TIME STEPS AND FORCES EVERY lOO 
■T . 

■TPLOl INC: PLOT INC ? CONDITION LABEL 

PLOT INC: 10 TSTORE RESULTS EVERY 10 TIME STEPS FOR PLOTTING 

T 

TINCMASS: NSTEPS , ANGLE7C0NDITI0N LABEL 

INCMASS: 1, .0875 TUPDATE MASS MATRIX EACH TIME FOR RIGID ELEMENT 

•f 


Figure 7.1-2 Input Data For Problem 7.1 
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«:ADMHSS: 

ADMASS: 


NDEl , 
1, 
1 , 
2, 

O 


$VEL,6RID PT, 
VEL ; 1 , 

^ If 
1 , 

JL. n 

— EDI /TOP — 


IDIRl , 

NDE2 , 

IDIR2, 

MASS 

7C0NDITI0N LABEL 



1 , 

1, 

1, 

<kZ 

•i^LUMPED MASSES ARE 

DEF INED 

AT EACH 

/-> 
X. , 

1, 


@3 

fNODE ACTING IN THE 

X AND Y 

DIRECTIONS 

1 , 

2, 

1 , 

@3 

i^WITH MAGNITUDE AS 

GIVEN IN 

THE DEFINE 



1 

©3 

^COMMAND 




DIR, COORD ,INIT VEL 

3, 0, ~&2 INITIAL ROTATIONAL AND NORMAL 

3, O, ~&2 $ VELOCITIES ARE DEFINED AT 

2, 0, +18.849556 ^GRID POINTS 1 AND 2 

2, 0, -18.849556 


Figure 7.1-2 continued 
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The TIMSTEP command indicates the integration,:-step ,'size for 
the transient response. Note t’hat both the TTMSP^N and TI'MSTEP 
commands must use a consistent set of time units. The s.pecified 
time step size for thi.s example is .001. The TIMSTEP com.mand can 
be modified via a condition label to provide variable step sizes 
during the analysis. 

Material and structural properties of the beam members are 
specified by the .MATPR-GJP command. -Each beam member has a unique 
identification name (.cprap.osed of any a 1 p.hanuraer i c characters) 
which is specified by the MEMBER command to designate the 
properties of th,e elements being defined. The Young's modulus (E) 
and the mass density/unit volume (RHO) are 

specified material properties whereas the cross sectional area (A) 
and the area moment of inertia (I) are specified structural 
properties. Note that all of the above quantities must be 
specified in a consistent set of units. In this sample problem, 
the material and structural properties of an aluminum bar with 
circular cross section and radius of 1 have been specified using 
the f oot / i nch/ pound/ second system of units. A unique feature of 
LATDYN is the ability to vary these properties during execution of 
the program by using a condition label associated with the MATPROP 
command. Further discussion on this feature can be found in Note 
5 .^- 5 . 

Grid points of the finite element model are defined using the 
GRID command (Note that in this manual the terms g£id po^jit and 
node are used interchangeably). Using a single GRID command the 
user can generate a series of equally spaced grid points between 
two specified points or define a single grid point. In this, 
example, the beginning and ending coordinates of two grid points 
are specified as XI, Y1, X2, Y2 and the starting grid point number 
is defined as grid point 1. Then since the number of grid points 
and the increment in numbering are specified as 2 and 1 
respectively, two grid points are defined. Note that although the 
GRID command can be used to define a string of grid points, the 
user must use a MEMBER command to ensure that the grid points are 
connected properly. 

Finite element beam members are defined through the MEMBER 
command. As in the previous command, this command can be used to 
define a string of members or a single member. The first 
parameter of the MEMBER command is an optional user defined member 
number which allows the user to specify the numbering of the. 
elements. If the member number is not specified then the program 
automatically assigns unique member numbers. In addition, the 
user provides the beginning and ending grid point numbers of, the 
first member and the number of members being defined with that 
particular command. If the number of members defined is equal to 
1, then a single beam member is defined between the two grid point 
numbers specified, however, if the number of members is greater 
than 1, then a string of members is defined as starting from the 
ending grid point of the first member and continuing until the 
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total number of members in that string is equal to that specified 
on the command. The increment between the beginning and ending 
grid point numbers of each of the remaining members in the string 
is assumed to be the same as that of the first member. In 
addition, the properties of each member are labeled with an 
identification name which must correspond with the identifier of 
a valid MATPROP command. In this example, one member is defined 
between grid points 1 and 2 and the properties of the member are 
given on the MATPROP command with the name "Aluminum". 

The optional commands used in this example are, 

DEFINE 

PRINT 

PLOTINC 

INCMASS 

ADMASS 

VEL 

Each of these are discussed below. 

The DEFINE command is a convenient way for the user to define 
constants which are used often in the data set. The form of this 
command is 


@n=Constant 

where the "n" is a unique integer from 1 to 99 (thus up to 99 
constants can be defined). Once a constant is defined, the 
shorthand definition can be used anywhere in the data case and the 
preprocessor will substitute the appropriate constant at that 
point. In this case, three shorthand definitions (@1, @2, and @3) 
appear in the MATPROP, ADMASS, and VEL commands. 

The PRINT command is the means for achieving printed output. 
The user can request the printing of displacements, velocities, 
and accelerations as well as internal forces (axial force, shear 
force, and bending moment) in either a local or global coordinate 
system. In addition, p r ov i s i ons • exi s t for the printing of the 
mass and stiffness matrices. With the PRINT command, the user 
specifies the increment at which each of the above quantities will 
be output. Thus in this example, the displacements, velocities, 
and accelerations of all the grid points will be printed every 100 
time steps in the global coordinate system, and the internal 
forces of all the members will be output every 100 time steps in 
the local coordinate system of each respective member. Also, no 
request has been made for the periodic printing of the mass and 
stiffness matrices. The PRINT command can be controlled via a 
condition label such that the user can vary the frequency and 
content of the output during different portions of the analysis. 

The ability to create plots is obtained by using the PLOTINC 
command. This command consists of a user specified integer, "r", 
which indicates that at every "r" time steps all of the results 
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will be stored on a seperate file which can be accessed through 
the postprocessor. The user is encouraged to read Section 5.19.^ 
for the details of the various plotting options available through 
the postprocessor. An important point to note is that th'e PRINT 
and PLOTINC commands are completely independent since they operate 
on different user specified increments. In this example, the 
results are sampled every 10 time steps so that the plots of this 
analysis will contain 10 times more data than that printed. The 
plotting increment on the PLOTINC command can also be changed by a 
condition label if desired. If the user saves the plot files as 
described in Section 5.19 and the Appendices, the postprocessor 
can be executed separately from LATDYN after the LATDYN execution 
iscomplete. 

Periodic updating of the system mass matrix is accomplished 
via the INCMASS command. The user specifies an integer "n" and 
the mass matrix is updated every "n" time steps (Note that the 
mass matrix will not be printed out unless the user specifies the 
same increment on the PRINT command). In addition, the user can 
specify the maximum angular rotation that any Mgi^d member can 
undergo without the mass matrix being updated. The default value 
of this rotation is .0875 radians (5 degrees). See Section 6.2 
for a discussion of the necessity of updating the mass matrix and 
its affect on system momentum. For this problem, the mass matrix 
is updated at every time step and the value of the rotation angle 
is specified as .0875. INCMASS can also be altered during 
execution by a condition label. 

Inertial coupling of degrees-df -freedom or lumped masses can 
be added to the model via the ADMASS command. In this example, 
lumped masses have been added to grid points 1 and 2 such that the 
masses affect the x and y directions but not the rotational 
direction. Note that the §3 which was defined by an earlier 
DEFINE command is used to specify the magnitude of the lumped 
mass. It is possible, through the use of a condition label, to 
add or remove the ADMASS command or to change the data on an 
existing command. See Section 6.2 for a discussion of the affect 
on system momentum due to mass additions during execution. 

Initial velocity conditions are imposed on the structure 
through the use of a VEL command. Velocities can be prescribed in 
the X, y or rotational directions either in the global coordinate 
system or in the local system of a designated member. In this 
case, both rotational and normal initial velocities have been 
prescribed for both grid points such that the beam will rotate 
exactly ir radians in one second. 

Several points can be made from an examination of the data in 
Figure 7.1-2. First, note the use of comments (which are preceded 
by a "$" sign) can be used anywhere in the data. Second, the free 
field format allows any number of spaces between data items. 

Thus, the number 1.0 E-3 would be interpreted as two pieces of 
data and should be input as 1.0E-3- Third, if a command is used 
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several times in succession, it is not necessary to repeat the 
command name since the program will assume all lines of data apply 
to that name until it encounters a new name (as demonstrated in 
the DEFINE, ADMASS and VEL commands of the sample problem). 

Fourth, a line of data can be continued to the next line by using 
a and completing the remaining data on the following line (as 

demonstrated with the GRID command). 

Sample Printed Outp ut 

Shown in Figures 7.1~3a through 7.1-3f is the printed output 
resulting from the execution of LATDYN using the data set of 
Figure 7.1-2. Figures 7.1-,3b through 7.1~3d show how the input 
data is printed out after it has been read and assimilated by 
LATDYN (Note that the order that the commands are output is not 
necessarily the same as they were input). It is often useful to 
compare this information with the data case to ensure that all of 
the data provided by the user is correct for the purpose it was 
intended. 

Following the input data printout, the program outputs any 
preliminary computations and presets which were carried out 
( Figur e 7 . 1 “3d ) . First the mass properties of the structure are 
printed. These entail the location of the centroid, the total 
mass, and the mass moments of inertia about, the x, y and z axes. 
Next, a summary of the constraints imposed on the structure is 
listed. For this problem, the maximum number of constraints is 
given as three since for a two noded rigid beam element only three 
of the six degr ees - of - freedom are independent. Then any preset 
displacements and velocities imposed by a DIS or VEL command are 
displayed. Also, information concerning the dynamic direct access 
memory usage is output. 

Next, as also shown in Figure 7.1-3d, an informative message 
is printed dealing with the integration procedure and the manner 
in which the displacements, velocities and accelerations are 
computed at each time step. It should be noted that most of the 
printed output mentioned above is typical of that which each run 
of LATDYN will generate, however, the printing of the assimilated 
input data and some of the preliminary computations and presets 
will vary from case to case depending on the input. 

Given in Figures 7.1~3e and 7.1-3f ane the results for this 
case at selected time increments. For every case in which a PRINT 
command is used, the output will always begin with the information 
known at time zero before the first integration step is taken. 

Then results will be printed according to the increment which the 
user specified on the PRINT Command. In this example, the 
displacements, velocities and accelerations as well as the forces 
were specified to be output at every 100 time steps, thus with a 
time step of 0.001 the next printout is at 0.100 seconds (Figure 
7.1~3e). For the case of a rigid beam, the values of the internal 
forces and moments are zero. The first two columns of results 
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listed under the heading "SOLUTION VECTORS" consist of the 
sequence number (the order in which the grid poi'n.ts were entered 
in the data) and the grid point number (the actual number which 
the user wishes assigned to the grid point). Listed under the 
subheading "GRID POINT TRANSLATION" and "GRID POINT ROTATION" are 
the X, y and phi displacements, velocities and accelerations of 
each of the grid points in its own local coordinate system. 
Directly underneath the local coordinate information are the 
values of "XLOC" and "YLOC" which are the coordinates of each node 
in the global system. Finally, the results listed under the 
heading "MEMBER DATA" are the rotations and strains of each beam 
member. 

Figure 7.1-3f shows the printout given for the last time step 
in the analysis which also includes a summary of the cpu time used 
in the major subroutines and ass erabli ng/ i nver t i ng the mass matrix. 
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Problem 7.2 


Flexible Boom With Rigid Hub And Applied Torque 


This sample problem involves the response of a flexible boom 
with a rigid central hub subject to an applied torque as shown in 
Figure 7.2-1. An antisymmetric boundary condition is used such 
that only one-half the structure is modelled. This sample problem 
demonstrates the application of external loads, boundary 
constraints, and condition labels. 

I n put Data 

Figure 7.2-2 contains the input data for Problem 7.2. 
Condition labels are used throughout the problem and the user is 
referred to Section 5.16 for a complete description of condition 
labels. The condition labels are written as, 

CL1 :T . GT . 1 . 0 
CL2:T.GT.5.0 

The first condition label (CL1) states that whenever the analysis 
time, T, exceeds the value of 1 second, the condition is "TRUE" 
and any commands which are dependent on this label are activated. 
Similarly for the second condition label (CL2) is "TRUE" whenever 
the time becomes greater than or equal to 5 seconds. The above 
implies that whenever a condition label is not "TRUE" it is 
"FALSE", thus any commands dependent on it are not active. In 
this example, condition label 1 is used in the TIMSTEP, PRINT and 
PLOTINC commands and condition label 2 is used in the STOP 
command . 

In addition to the required and optional commands used and 
discussed in Problem 7.1, the following optional commands are used 
in the current problem,' 

PLOTLIMIT 
INTGTYP 
STOP • 

APPLOAD 

LOAD 

SDFC 

SCALE 

The PLOTLIMIT command indicates a specific time duration for 
which the results are stored for plotting with the postprocessor . 
This command is used in conjunction with the PLOTINC command such 
that the user can obtain a large amount of data for plotting, but 
only during a specified portion of the analysis. 

The INTGTYP command specifies the algorithm used in the 
transient analysis. In this example, the integration procedure 
used is central differences (INTGTYP=1). 
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Figure 7.2-1 Flexible Boom With Rigid Hub And Applied Torque 
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s 

tTITLE; DATA CASE TITLE (UP TO 80 CHARACTERS) 

TITLE: T2S - SECOND TEST CASE FOR LATDYN USER'S GUIDE 
■$ 

^SUBTITLE: DATA CASE SUBTITLE (UP TO 80 CHARACTERS) 

SUBTITLE: FLEXIBLE BOON WITH RIGID CENTRAL HUB AND APPLIED TORQUE 
^CONDITION LABELS DEFINITION 

CLl: T.GT. 1 .0 f FIRST CONDITION TRUE WHEN TINE EXCEEDS 1 SECOND 
CL2:T.GE.5.0 -^CONDITION TRUE WHEN TINE REACHES 5 SECONDS 
■$ 

^DEFINE: @N=CONSTANT 

DEFINE: @1 = 11.2E06 ^^DEFINE VALUE OF YOUNG'S NODULUS 
DEFINE: @99=3 . 3376E09«:UP TO 99 CONSTANTS CAN BE DEFINED 
# 

BORDER: NELN, NODES 

BORDER: 6, 7 i^DEFINE UP TO 6 ELENENTS AND 7 NODES 

t- 

■fT INSPAN: ST ART, STOP 

TINSPAN; 0.0, 10.0 ^ANALYSIS INTENDED TO RUN FOR 10 SECONDS 

tTINSTEP: DELTA T 7C0NDITI0N LABEL 
T INSTEP: 3.0E-5 ^ORIGINAL TINE STEP 

3.0E-4 ?1 «:N0DIFY TINE STEP IF CONDITION ONE IS TRUE 

f:NATPROP: NATNANE,NTYPE, E , A , I , RHO ? CONDIT ION LABEL 
NATPROP: 1, O, @1, 48.440, 193.60, 2.5B8E-4 tDEFINITION OF NATERIAL 

tAND STRUCTURAL PROPERTIES WITH NATERIAL NANE 1 

tGRID:NST,# NODES, XI , Y1 , X2 , Y2 , INC 

GRID: 1,. 7, 0.0, 0.0, 100.0, 0.0, 1 -^SPECIFY 7 NODES EQUALLY SPACED 

tFRON ZERO TO 100 ALONG THE X-AXIS 
■$ 

i^NENBER; INENB ,N1 ,N2 , NE ,NAT 

NENBER: ,1, 2, 6, 1 ^SPECIFY,6 ELENENTS CONNECTING THE 7 NODES 

f- 

■S 

•SPRINT; DV A INC, COORD,INTF INC, COORD, NN INC,SN INC ?CONDITiaN LABEL 
PRINT: 16666, 0, 16666, L, 0, 0 :fPRINT RESULTS EVERY 16666 

1666, 0, 1666, L, 0, 0 '?! f DECREASE PRINT FREQUENCY 

«:WHEN THE TINE STEP IS INCREASED -IF CONDITION 1 IS TRUE 
■s 

#PLOTINC! PLOT INC ? CONDITION LABEL 

PLOT INC: 1666 tSTORE PLOTTING INFO. 10 TINES AS OFTEN AS PRINTOUT 

166 71 ^CHANGE PLOT INC CONSISTENT WITH PRINT INCRENENT 

tPLOTLINn -.LOWER TINE LIN IT, UPPER TINE LINIT 

PLOTLINIT: 0.0, 5.0 tONLY DESIRE PLOTS FOR FIRST 5 SECONDS 

Figure 7.2-2 Input Data For Problem 7.2 
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f BATCH PLOTTING COMMANDS 
SNAPSHOT , 0.0 ,5.0, 400^0 
PL0TY,D,7 
PLOTPH 1 , D , 1 

«:INCMASS: NSTEPS , ANGLE7.C0ND I T I ON LABEL 

INCMASS: O, 0.0875 #N0 UPDATING OF THE MASS MATRIX IS REQUIRED 

S 

TINTGTYP: INTG,(BETA) 

INTGTYP: 1 

■I- 

#STOP: 7 CONDITION LABEL 

STOP: 72 .rSTOP ANALYSIS WHEN CONDITION 2 IS TRUE 
t 

#ADMASS: NDEl, IDIRl , NDE2, IDIR2, MASS 7C0NDITI0N LABEL 

ADMASS: 7, 1, 7, 1, 248.0 SLUMPED MASS AT NODE 7 IN THE 

7, 2, 7, 2, 248.0 SX AND Y DIRECTIONS 

1, 3, 1, 3, @99 SADD ROTATIOMAL INERTIA TO NODE 1 

S 

SAPPLOAD: LOAD ID, IMEMB ,DIR , SCALE , tSTART , TSTOP , GR I D PNTS 7 CONDITION 
APPLOAD: 1, , TORQUE, 10.0, 0.0, 10.0, 1 SAPPLY TORQUE TO NODE 1 

S 

SLOAD:LOAD ID, AMPLITUDE 

LOAD: 1, 1.0904E07 SSPECIFIES VALUE OF TORQUE USED BY APPLOAD 

S 

s 

SSDFC:GRID NO , CON DIR, CONSTANT DOF 7 CONDITION LABEL 
SDFC: 1, 1, 0.0 SCONSTRAIN MOTION OF NODE 1 IN X DIRECTION 

1, 2, 0.0 tCONSTRAIN MOTION OF NODE 1 IN Y DIRECTION 

S 
S 

SSCALE: SCALEX , SCALEY , ROT 

SCALE: 20.0, 1.0, 0.0 SSCALE GEOMETRY IN X DIRECTION BY 20 

SHAS THE NET EFFECT OF MAKING THE TOTAL STRUCTURE 2000 UNITS LONG 


Figure 1.2-2 continued 
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As mentioned earlier, the STOP command can be used to 
terminate program execution prior to the stop time specified on 
the TIMSPAN command. The program can be stopped anytime during 
execution once the condition label associated with the STOP 
command becomes true. Only one STOP command can be included in 
the data set. In this example, the program is stopped once 
condition label 1 becomes true. 

Using the APPLOAD command the user specifies the direction, 
time duration and location of an externally applied load. The 
user has the option of specifying the direction of an applied load 
in terms of the local coordinate system of any member, if not 
specified then the direction is assumed to be relative to the 
global coordinate system. This command is required when a load is 
to be applied to the structure (with the exception of GRAVITY) and 
must be used in conjunction with any of the LOAD, LOADSIN, or 
LOADTIM commands. In this example, APPLOAD is used along with the 
LOAD command to apply a torque at node 1 which is to act for 10 
seconds. The scale factor of 10 is applied to the magnitude of 
the load specified on the LOAD command. The application of loads 
through the APPLOAD command can be controlled with condition 
labels to turn loads on and off during the analysis. 

The LOAD command specifies a constant load with a given 
magnitude. Each LOAD command has a load identification number 
which is also used by the APPLOAD command to designate the 
location, direction and scaling of the load. For this case, the 
load identified as "1" is a torque at node 1 with a magnitude 
which will be scaled by 10 when it is applied via the APPLOAD 
command. 

Single degree-of -freedom and boundary constraints are 
specified with the SDFC command. The user can specify that the 
motion of a grid point be fixed at a given value. In this 
example, the x and y motions of grid point 1 are constrained to 
zero although it is allowed to rotate. Thus, grid point 1 has 
effectively been pinned to simulate the antisymmetric boundary 
condition desired. It is possible to make the SDFC constraints 
conditional such that they can be imposed or released throughout 
the analysis. 

By using the SCALE command the user can scale the dimensions 
of the structure in the global x or y directions as well as rotate 
the entire structure. This command affects al l of the coordinates 
of the grid points defined by the GRID command. This problem uses 
SCALE to increase the length of the structure by 20 in the x 
direction while maintaining the original dimensions in the other 
directions. 


1 1 -15-85 


7.2-5 


Sample Problem 2 



Sample Plots 


Figures 7.2-3 through 7.2-5 show the plots generated (via the 
postprocessor ) by the plotting commands SNAPSHOT, PLOTY, and 
PLOTPHI. These commands are included in the input data set for 
this sample problem a‘s described in Section 5.19.^. 

The command u’s'ed to produce Figure 7.2-3 was, 

SNAPSHOT, 6.0,5. 0, 400,0 

and the result is a series of structural deformation plots 
overlaid . 

Figure 7.2-4 depicts the time history of the displacement of 
node 7 parallel to the global y axis and was produced by the 
command , 

PLOTY, D, 7 

Finally, a time history of the rotational motion of node 1 is 
given in Figure 7.2-5. The command used to obtain this plot was, 

PLOTPHI, D,1 


11-1 5-85 


7.2-6 


Sample Problem 2 



FIN TINE=5. 00000 
INI TIMErO. 00000 


I 



Figure 7.2-3 Plot Created by SNAPSHOT Command 
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Figure 7.2-4 Plot Created by PLOTY Command 
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Figure 7.2-5 Plot Created by PLOTPHI Command 
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Problem 7.3 - Unfolding Deployment of One Bay 


The problem modelled in this example is the unfolding 
deployment and lock-up of two flexible beams as depicted in Figure 
7 . 3 ^ 1 . The application of a rotational spring at a pin joint is 
demonstrated as well as the method for locking a joint. In 
addition, user defined variables are introduced. 

I nput P at a 

The data set for this problem is shown in Figure 7.3~2. All 
of the required commands and most of the optional commands which 
were discussed in the previous problems have been used in this 
example. The following optional commands are also used, 

SPRING 

PIN 

ROTLOCK 

SET 

RESET 

The SPRING command is used to define a spring connected 
between two grid points. Springs can be rotational or extensional 
and can act in a fixed direction or follow the line of action 
between two grid points. In addition to the grid point numbers 
and the direction, the user specifies the spring constant and any 
initial tension/ compression (positive for tension and negative for 
compression). In this example, a rotational spring ( di r ect i on=3 ) 
is defined between points 2 and 3 with a spring constant of 1000.0 
and an initial compression of t: radians. Thus, the torque 
provided by the spring is, a linearly decreasing function of the 
rotation of the two nodes such .that the torque reaches zero when 
the total rotation of the nodes is ir (see sketch below). 


Torque 

(N-m) 



Total Rotation 


Springs can be added or removed at anytime during the analysis by 
using a condition label. Thus, it is possible to synchronize the 
deployment of a multi-bay truss by activating the deployment 
springs at prerarranged times. 
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7.3-1 Unfolding Deployment Of One Bay 
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■$ 

f TITLE: DATA CASE TITLE (UP TO 80 CHARACTERS) 

TITLE: T IS - THIRD TEST CASE FOf< LATDYN USER'S GUIDE 
S- 

SSUBTI TLE; DATA CASE SUBTITLE (UP TO 80 CHARACTERS) 

SUBTITLE: UNFOLDING DEPLOYMENT OF ONE BAY 
■$ 

f DEFINE: @N=CQNSTANT 

DEFINE: @1=3.14159265 tDEFINE PI 

«:P0RDER: NELM, NODES 

P0RDER:2, 4 ^MAXIMUM OF 2 ELEMENTS AND FOUR NODES 
STIMSPAN: START, STOP 

TIMSF'AN: 0.0, 5.0 fFIVE SECOND ANALYSIS 
■$ 

fTIMESTEP: DELTA T 7C0NDITI0N LABEL 

TIMSTEP: . 25E-3 .fTIME STEP OF 0.00025 SECONDS 

J:MATPR0P:MATNAME,MTYPE, E , A , I , RHO ? CONDITION LABEL 
MATPROP-.TYPEl , O, 70. E9, .0003688, 5.7625E-6, '2712.0 t-DEFINIT ION 

:SOF MATERIAL AND STRUCTURAL PROPERTIES WITH NAME "TYPEl" 


J^GRID:NST,# NODES, 

XI 

, VI 

Y 

, A , 

Y2 ,INC 

GRID: 1, 2, ( 


0, 

0, 

25, 1 ^DEFINE TWO NODES 1?< 

3, 2, ( 


25. , 

0, 

0, 1 ^DEFINE TWO NODES 3^!/. 

4^- 

2TMEMBER: IMEMB ,N1 


N2 , NE 

,MAT 


MEMBER: ,1, 2, 


TYPEl 

t DEFINE 

MEMBER FROM NODE 1 TO NODE 2 

,3, 4, 

^ 9 

TYPEl 

$DEF INE 

MEMBER FROM NODE 3 TO NODE 4 




SPRINT: DVA INC, COORD.INTF INC, COORD, MM INC.SM INC 7C0NDITI0N LABEL. 
PRINT: 1000, 0, 1000, L, 0, 0 *:PR1NT EVERY 1000 

^PLOTINC: PLOT INC 7C0NDII ION LABEL 
PLOT INC: 100 #PLOT EVERY 100 STEPS 
■I 

SBATCH PLOTTING COMMANDS 

SNAPSHOT , O . O , 3 . O , 1 2 , 0 , A 

PLOT ROT, 2 

PL0TPH1,D,3 

PLOT 0,3, 5 

PLOT 0 , 5 

•TINCMASS: NS TEPS, ANGLE 7C0NDITI0N LABEL 

I NCMASS: 200000, .0875 4: INCREMENT MASS MATRIX EVERY 200K STEPS 

■$ 


Figure 7.3~2 Input Data For Problem 7.3 
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«:SPR I NG : NDE 1 , NDE2 , I D I R , SPNGK , PSSPNG ? COND I T I ON LABEL 

SPRING: 2, 3, 3, 10pp., -@1 iOEFINE ROTATIONAL SPRING BETWEEN 

TNODES 2 AND 3 WITH MAGNITUDE OF 1000 

t 

■T ... 

«:p I N : GR I DPT 1 , GR I DPT2 ?CONb I T I ON LABEL 

PIN: 2, 3 :^DEFINE PIN JOINT BETWEEN NODES 2 AND 3 

# 

•fCLI: CONDITION LABEL DEFINITION 

CLl: ABS(PHI (D,2)-PHI (D,3) ) .GE.Q2 '^CONDITION TRUE WHEN NODES 
$: 2 AND 3 HAVE ROTATED THROUGH A SUM OF PI RADIANS 
CL2: NTSTEP.GE.l -^CONDITION 2 TRUE AT FIRST TIME STEP 
■S- ^ 

«:R0TL0CK:GRD PTl , GRD PT2 7C0NDITI0N LABEL 

ROTLOCK: 2, 3 ?1 ^ROTATIONAL LOCKUP OF NODES 2 AND 3 WHEN 

J:CONDITION 1 IS TRUE 
•f 

SSET: QN=FORTRAN 7C0NDITI0N LABEL 

SET: Q1=ABS<YL0C<1)-YL0C(2) ) 72 :TDEFINE ELEMENT LENGTH AT TIME ZERO 
SET: Q2=4. *ATAN ( 1 . O) «:ALTERNATE DEFINITION OF PI 
t- 

PRESET: QN=FORTRAN 7C0NDITI0N LABEL 

RESET: Q3=ABS (PHI (D , 2) -PHI (D , 3) ) ^TOTAL ROTATION BETWEEN NODES 2 AND 3 
RESET: Q4=ABS ( YLOC ( 1 ) -YLOC (2 ) ) 72 «:RELATIVE DISPLACEMENT BETWEEN NODES 2 AND 
RESET: Q5=Q4/Q1 72 ^:CHANGE IN DISPLACEMENT NORMALIZED BY ORIGINAL LENGTH 


Figure 7.3"2 continued 
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The PIN command is a convenient method of creating a pin 
joint between two grid points. This command is used in lieu of 
two MDFC commands to constrain the translational degrees-of- 
freedora of the points together while allowing them to rotate 
independently. In this example points 2 and 3 are pinned 
together. It is also possible to create or remove pins during the 
analysis with condition labels. 

Using the ROTLOCK command, the user can simulate the 
rotational lock-up of a joint. This option is used when it is 
desired to lock a joint without modelling the details of the joint 
as discussed in Section 5.17.2. In this example, the pin joint 
between grid points 2 and 3 is locked up when condition label 1 is 
"TRUE". From the CL1 command, condition label 1 is "TRUE" when 
the total rotation of both nodes exceeds the user defined variable 
Q2, which is defined as ir (see discussion of SET command below). 

The SET and RESET commands are user defined commands which 
allow the specification of variables in the form of FORTRAN 
expressions which can involve any of the physical quantities 
measured by the program at any time during the execution (the user 
is referred to Section 5.17 for a more complete description of 
these commands). In this ' example, the SET command is used to 
compute the length of one element as the absolute value of the 
difference in the y locations of nodes 1 and 2 in the global 
coordinate system. This command is activated when condition label 
2 is "TRUE", i.e,, at the first time step (see CL2 command in : 

Figure 7.3“2). Note that once activated, the value of the Q- 
variable specified by a SET command remains constant, and. can then 
be used elsewhere in the data set. The second SET command is used 
to define the value of ir which is used in the CL1 command. 

Also used in this example are three RESET commands. Once 
activated, the variables specified by each of the RESET commands 
are recomputed at every time step, in contrast to the SET command 
whose Q-variable iS’ evaluated only once. The first RESET command ■ 
computes the total rotation between points 2 and 3 at every time 
step (note that the variable Q3 could replace the lefthand side of 
the expression in the CL1 command). The second RESET command 
specifies the computation of the relative displacements of points 
1 and 2 as was done in the first SET command. Finally, the third , 
RESET command allows the calculation of the change in the relative 
displacements . of points 1 and 2 normalized by the original length 
by dividing by Q1. 
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Sample Plots 


Plotting is requested in this problem due to the presence of 
the SNAPSHOT, PLOTPHI , and PLOTQ commands in the input data 'set. 
Figures 7.3"3 through 7.3"5 contain the plots generated .by these 
commands . 

The deployment sequence of the two beams is depicted in 
Figure 7.3~3 by't'he s'e’r i es of ' o ver 1 ai d structural deformation 
plots at various times. The command which generated this plot 
was , • ' 


SNAPSHOT , 0 . 0, 3 . 0, 1 2, 0 , A 

Figure 7.3"^ is a time history of the rotation of node 3 
which was produced by the command, 

PLOTPHI, D, 3 

The magnitudes of the two user defined variables Q3 and Q5 
are depicted in Figure 7.3~5. The command used to plot the Q- 
variables was, 

PLOTQ, 3, 5 

Note that the joint lock-up time is the time at which the user 
defined variable Q3 reaches -n radians. 
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Problem 7.^ - Rotating Rigid Beam With Control Forces 


Problem 7.4 involves a rigid beam element with an initial 
rotational velocity of -ir rad/sec as in Problem 7.1. Q-variables 
are used to construct control forces normal to the' beam and 
proportional to its velocity which act to slow it down first to - 
it/ 2 rad/sec and then to -ir/4 rad/sec. Of interest in this problem 
is the use of condition labels to activate and deactivate the 
control forces as well as the user defined variables to provide 
the appropriate direction and magnitude of the forces. 

The basic model for this problem is shown in Figure 7.4-1 and 
the input data is given in Figure 7.4-2. There are five condition 
labels specified in the data. These are used to specify the 
instances at which the control forces are activated. Condition 
label 1 is "TRUE" when grid point 1 has rotated through an angle 
equal to the value specified by the Q-variable Q5, namely 1 0 
degrees. Condition label 2 becomes "TRUE" when CL1 is "TRUE" and 
the rotational velocity of point 1 is greater than or equal to tr/2 
(note the use of the logical .AND. in the condition label). 
Effectively, CL2 is used to turn on the control force when the 
rotational velocity of the beam is greater than that desired. 

Once the new velocity state is reached, CL2 becomes "FALSE" and 
the control force is removed. 

Condition labels 3, 4 and 5 are used in a similar manner to 

allow the beam to reach the second rotational velocity state 
specified by Q6 and Q7, whereQ6 and Q7 are defined through the 
SET command . . 

The SET command is used to define the desired rotational 
velocity states specified in the condition labels as Q6, Q7 and 
Q8. These values are computed once at the first time step and 
their values remain constant throughout the analysis. 

Using the RESET command, the angle of the velocity vector at 
grid points 1 and 2 is computed at each time step as the variable 
names Q1 and Q3. These angles are used in the APPLOAD command to 
specify the direction of the control forces normal to the beam. 

In addition, Q2 and Q4 are computed as the magnitude of the 
velocities of points 1 and 2 multiplied by a scale factor of 0.01, 
and are the scale factors for the control forces specified in the 
APPLOAD command. 

The control forces are specified via the APPLOAD and LOAD 
commands. The APPLOAD commands specify the direction, scale 
factor, time duration and location of the control forces. The 
first two APPLOAD commands are activated when condition label 2 is 
"TRUE" and deactivated when the desired velocity state specified 
on the CL2 command is reached. Likewise, the next two pairs of 
APPLOAD commands specify the application of the control forces 
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when conditions ^ and 5 respectively are ''TRUE,'' . N.d'te' that the 
direction of the control forces is given by the angle of the 
velocity vector normal to the beam (Q1 and Q3) . ah’tf-- the scale 
factor is proportional to the magnitude of the velocity of the 
beam (Q2 and Q4). Each APPLOAD' command is' ass i gned; a- load 
identification number;, whi ch' is. ref erenced by a particular type of 
load. In this exampl'^^V the lda,d‘ type is a constant Toad of 
amplitude 1.0 given’ b'y-, the L^OA'd;, command . Thus, the scale factor 
used in the APPLOAD command (which varies' wi'th the magnitude of 
the beams velocity') is used to scale the amplitude' of the LOAD 
command thereby creating an external load whose magnitude can vary 
ateachtimestep. 

Sampl e_P ]L ot s 

In this example problem, 5 plots are' requested through the 
plotting commands SNAPSHOT, PLOTC, PLOTMAG, PLOTQ, and PLOTROT 
which are in the input data set. The plots are shown' in Figures 
7.A-3 through 7.^-7. 

A series of overlaid plots of the beam at various times are 
shown in Figure 7.A-3. These results were produced by the 
command , 

SNAPSHOT , 0 . 0 , A . 0 , 55, 0 

The logical status of the 5 condition' labels during the 
analysis is shown in Figure 7.^-4, an'd the command which generated 
this plot was , 

PLOTC, 1 , 2, 3, 5 

Figure 7.A-5 shows a time history of the magnitude of the 
velocity vector at node 1 which was created by the command, 

PLOTMAG, V,1 

The magnitude of the user defined Q-variable Q2 during the 
analysis is presented in Figure 7.^-6. The command used to 
produce this plot was, 

PLOTQ, 2 

Finally, the rotation of the beam element is depicted in 
Figure 7.A-7 and this corresponds to the command, 

PLOTROT,! 
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V 2 = -18.85 in/sec 


(a) INITIAL CONDITIONS 



(b) APPLICATION OF CONTROL FORCES 


Figure 7.^-1 Rotating Rigid Beam With Control Forces 
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■f 

tlHlS PROBLEM IS THE SAME AS THE FIRST PROBLEM BUT IT HAS CONTROL 
^FORCES WHICH ARE PROPORTIONAL TO THE VELOCITY 

■f- 

^TITLE: DATA CASE TITLE ( UP TO 80 CHARACTERS) 

TITLE: T9SA - FOURTH TEST CASE FOR LATDYN USER'S GUIDE 
# 

iJSUBTITLE: DATA CASE SUBTITLE (UP TO 80 CHARACTERS) 

SUBTITLE; ONE RIGID ELEMENT WITH INITIAL VELOCITY AND CONTROL FORCES 

■t- 

^DEFINE: @N=CONSTANT 

DEFINE; @1=. 0002588 ^DEFINE MATERIAL DENSITY 

DEFINE: @2=3.141592654 * DEFINE'pI 

(§3= l.OE-3 ^DEFINE VALUE OF CONCENTRATED MASS 
•T 

«:CLI; CONDITION LABEL DEFINITION 

CLl: ABS(PHI (D, 1) ) .GE.Q5 «TRUE WHEN ROTATION OF NODE 1 = Q5 
CL2; CLl . AND. ABS (PHI (V, 1 )). GE. PI/2. 0 tTRUE WHEN CLl IS TRUE AND 
■TTHE ROTATIONAL VELOCITY OF NODE 1 IS GREATER THAN PI/2 
CL3: ABS(PHI (D, 1) ) .GE.Q6 «^TRUE WHEN NODE 1 ROTATES GREATER THAN Q6 
CL4; CL3. AND. ABS(PHI (V, 1) ) .GE. PI/4.0 #TRUE WHEN CL3 IS TRUE AND' 

^THE ROTATIONAL VELOCITY OF NODE 1 IS GREATER THAN PI/4 

CL5: ABS(PHI (D, 1) ) .GE.Q7 $TRUE WHEN NODE 1 ROTATES GREATER THAN Q7 

«:P0RDER; NELM, NODES 

PORDER; 1, 2 TONE BEAM ELEMENT AND TWO NODES ARE DEFINED 

T 

TTIMSPAN: START, STOP 

TIMSPAN: 0.0, 4.0 T ANALYSIS TO RUN FOR 4 SECOND 

T 

TTIMSTEP: DELTA T 7C0NDITI0N LABEL 

TIMSTEP: l.OE-3 T INTEGRATION TIME STEP IS .001 SECONDS 

T 

TMATPROP:MATNAME,MTYPE, E , A , I , RHO ? CONDITION LABEL 
MATPROP: ALUMINUM, 3, 10.5E06, @2, .78540, @1 T USING RIGID MEMBER 

T WITH ALUMINUM MATERIAL PROPERTIES 
T 

TGRID:NST,# NODES, XI , Y1 , X2 , Y2 , INC 
GRID: 1 , 2 , 0.0, 0 . 0 , 1 2 . 0 , S< 

0.0, 1 TTHIS IS A CONTINUATION OF THE LINE ABOVE — TWO GRID POINTS 

T HAVE BEEN DEFINED 
T 

TMEMBER: IMEMB ,N1 , N2 , NE , MAT 

MEMBER: ,1, 2, 1, ALUMINUM TONE MEMBER FROM NODES 1 TO 2 USING 

TIHE PROPERTIES OF “ALUMINUM" FROM THE MATPROP COMMAND 
T 
T 

TPRINT:DVA INC, COORD , INTF INC , COORD , MM INC,SM INC 7C0NDITI0N LABEL. 
PRINT: 100, 0, 200, L, O, O TPRINT DISPLACEMENTS/ 

TVELOCITIES/ACCELERATIONS EVERY 100 TIME STEPS AND FORCES EVERY 200 
25, 0, 25, L, 0, O 72 TPRINT MORE OFTEN 

TWHEN CONTROL FORCES ARE ACTIVATED 


Figure 7.H-2 Input Data For Problem 7.4 


1 1 - 15-85 


7.4-4 


Sample Problem 4 



tPLOTINC: PLOT INC ? CONDITION LABEL 

PLOT INC: 10 «STORE RESULTS EVERY 10 TIME STEPS FOR PLOTTING 

3 ?2 SSTORE RESULTS MORE OFTEN WHEN CONTROL FORCES ON 
t 

»BATCH PLOTTING COMMANDS 
SNAPSHOT ,0.0,4.0,55,0 
PLOTC, 1 ,2,3,4,5 
PLOTMAG,V, 1 
PL0TQ,2 
PLOTROT , 1 

tINCMASS: NSTEPS , ANGLE7C0NDI T ION LABEL 

INCMASS: 1, .0875 ^UPDATE MASS MATRIX EACH TIME FOR RIGID ELEMENT 

$ 


^rtDMASS: 

NDEl , 

IDIRl , 

NDE2, 

IDIR2, 

MASS 

7C0NDITI0IM LABEL 



ADMASS: 

1, 

1, 

1, 

1, 

@3 

:fLUMPED MASSES ARE 

DEFINED 

AT EACH 


1, 

2, 

1, 

2, 

@3 

$NODE ACTING IN THE 

X AND Y 

DIRECTIONS 


2, 

1, 

2, 

1, 

@3 

MWITH MAGNITUDE AS 

GIVEN IN 

THE DEFINE 




2, 

2, 

@3 

^COMMAND 




S 

*APPLOAD:LOAD ID,IMEMB , DIR , SCALE , TSTART, TSTOP , GRID PNTS 7C0NDITI0N LABEL 
APPLOAD: 1, , Ql, 02, 0.0, 4.0, 1 72 tAPPLY CONTROL FORCES 

1, ,03, Q4, 0.0, 4.0, 2 72 tAT NODES 1 AND 2 WHEN 

*CL1 .IS TRUE, WITH DIRECTIONS NORMAL TO THE BEAM AND SCALE FACTORS 
^PROPORTIONAL TO THE MAGNITUDE OF THE VELOCITY 



1, 

,01, 

02, 

0 . 0 , 

4.0, 

1 

74 


1, 

,03, 

04, 

0.0, 

4.0, 

2 

74 


1, 

,01 , 

02, 

0.0, 

4.0, 

1 

75 

^LOAD: LOAD 

1, 

ID, 

,03, 

AMPLITUDE 

04, 

0.0, 

4.0, 

2 

75 

LOAD: 1, 


1.0 i^DEFINES 

CONSTANT 

LOAD 

WITH 

AMPLITUDE 


* 

■$ 

tSET: QN=FORTRAN 7C0NDITI0N LABEL 


SET: 05=10. 0*PI/ 180.0 
06=90. 0*PI / 180. 0 
Q7=P1 
t- 


PRESET: QN=FORTRAN 7C0NDITI0N LABEL 

RESET: Ql= ANGLE(V,1) ifCOMPUTE DIRECTION OF CONTROL FORCE AT NODE 1 
02= -MAG(V,l)*1.0E-2 t-DEFINE SCALE FACTOR FOR FORCE AT NODE 1 
03= ANGLE<V,2) ^COMPUTE DIRECTION OF CONTROL FORCE AT NODE 2 
Q4= -MAG(V,2)*1.0E-2 $DEFINE SCALE FACTOR FOR FORCE AT NODE 2 
■t- 


*VEL,GRID PT, 
VEL: I, 



'"J 


DIR, 

3, 



COORD , INI T VEL 

0, -@2 $ INITIAL ROTATIONAL AND NOf<MAL 

0, -@2 $ VELOCITIES ARE DEFINED AT 

0, +18.849556 ^GRID POINTS 1 AND 2 
O, -18.849556 


Figure 7.4-2 continued 
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Figure 7.M-3 Plot Created by SNAPSHOT Command 
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Figure 7.^-^ Plot Created by PLOTC Coromand 
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Figure 7.^-5 Plot Created by PLOTMAG Command 
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Figure 7.4-7 Plot Created by PLOTROT Command 
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APPENDIX A; LATDYN STRUCTURE 
Components of LATDYN 

LATDYN consists of three components which operate in 
combination with a FORTRAN compiler and are linked through data 
files and job control language commands. The components are the 
pre-processor, the LATDYN calculation routines, and the post- 
processor. The pre-processor acts as the command language 
interpreter for the program. The post-processor provides the 
graphics portion of the program. 

When the LATDYN program is activated, using the instructions 
provided in Appendix B for CDC/NOS execution or Appendix C for 
VAX/VMS execution, all or some of these components are utilized 
depending on user selected options. 

When LATDYN is activated the user selects the components of LATDYN 
to be executed. For CDC/NOS equipment, the options are selected as 
arguments of a system control command as shown in Appendix B. For 
VAX/VMS equipment, the user is queried on the screen concerning 
option selection. 

In addition to the three components, LATDYN also makes use of 
a FORTRAN compiler. The use of a FORTRAN compiler gives the user 
access to the power of a high level language, albeit in a highly 
structured environment. This requires that the FORTRAN compiler be 
incorporated as part of the execution process of LATDYN and that 
there be a means of linking the compiled code with the main 
program. For CDC/NOS and VAX/VMS equipment these requirements are 
met and utilized automatically. 

Pre-processor Component 


The principal function of the pre- processor portion of the 
program- is to act as a command language interpreter. It allows the 
user to modularize the input data by permitting multiple input 
files. It also serves to check for many syntactical errors which 
may exist in the input data commands. As described later under 
this heading, it cannot detect all types of errors, but will detect 
and flag all command language syntactical errors. 

The control instructions which activate the pre-processor 
depend upon the computer system being used. Appendices B and C 
provide these instructions for the CDC/NOS and VAX/VMS computers. 
For the CDC/NOS machines, the pre-processor can be used in either 
an interactive or batch mode. For the VAX/VMS machine, the pre- 
processor is used in an interactive mode. Usually the interactive 
mode will prove most useful. 
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“ After" latdyn has been 

activated and the pre-processing' option in the i’nteracti've mode 
selected, the pr e-processor * s' first action is- to' ask for the name 
of the input data file. This query appears on the- screen as: 

"Where are the input directives?" 

The response;' to j^t'h is question is a file name or the word "ME". 
If the response to tlhTs" que’st i on i's- "ME" then input data commands 
are to be supplied rtfanual^ry by' the user directly from the terminal. 
The user will be inStructe'd when to supply this' data. 

The pre-prdcessoh creates' certain output which may be useful 
to the user. Pre-pr ocessor output is divided into three 
categories, 

* An echo of all input command data 

* Warning messages 

* Fatal errors 

Hence, the next action of the pre-processor is to inquire about the 
output it creates. 

In the interactive mode the program asks for the name of a 
file for echoing input data: 

"Where should the directives be echoed?" 

If the response to this question is the word "ME" (or a 
carriage return), then all pre-processor output (not to be confused 
with calculated output) is directed to the screen. 

If the response is a file name, then input directives will be 
echoed to that file. When run interactively, the pre-processor 
always outputs warning and fatal messages to the screen, but in 
addition the user can place these on files. The program will ask, 

" Where should any warning messages be written?" 

" Where should any fatal error messages be written ?" 

The user should respond to these questions with file names or if 
only screen messages are desired, with carriage returns. Saving 
files of these messages can prove useful for large problems wherein 
many errors may exist. 
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If the answer to the question "Where are the input 
directives?" was "ME", a question mark will appear on the screen 
after the pre-processor output questions have been answered. The 
user should then enter the individual LATDYN commands following 
each with a carriage return. The program will respond with an echo 
of the command if the answer to the question "Where should the 
directives be echoed?" was "ME". A question mark on the screen 
indicates that the next command should be entered. When all input 
data commands have been entered, the user should respond with a 
single carriage return to terminate manual entry of input data 
commands . 

Whether input data is supplied from a file or from the screen, 
the user will be queried for additional input data with the 
question, 

"Where are any additional directives?" 

The user responds with a file name, "ME" or "NO". This 
additional request allows input data to be supplied via multiple 
files or a combination of files and screen input. 

Pre-processing in the Batch Mo de - This option is only 
available on the CDC/NOS machines. The instructions for activating 
the pre-processor in the batch mode are provi ded on page B-1. 
Arguments in the instruction command provide the name of an input 
data file which contains the LATDYN commands. 

Files Created by the Pre-pro cessor - When all input data has 
been supplied, the program will process the data and create the 
files LATIN, UDAC and POSTIN. The file LATIN contains the pre- 
processed instructions for the calculation portion of LATDYN. The 
file UDAC contains any FORTRAN coding supplied by the user in Q- 
variable definitions and condition statements. UDAC is compiled 
and linked to LATDYN in the form of subroutine UDAC. In addition, 
if the input file called for graphics output, a file named POSTIN 
is created. POSTIN contains all the plotting commands of the input 
data. However, if a fatal error is detected in the input data, 
none of these files will not be created. 

When the program has completed pr e-processi ng all of the input 
data and no fatal errors have been detected the program will inform 
the user that the pr e-processor has completed its task. 

Note that the existence of these files does not mean that the 
data is free from all errors. The pre-processor is able to detect 
certain syntactical errors in the input data. However, there are 
errors that the pr e-processor is unable to detect. The pre- 
processor cannot detect FORTRAN errors in the Q-variables and the 
condition statements. Also it cannot check to see that the number 
of members and grid points used in the model does not exceed that 
given in the PORDER command. These types of errors will be 
detected in the FORTAN compiler or in the LATDYN pre-processor . 

The pr e-processor will inform the user of those errors which it can 
detect . 
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If a fatal FORTRAN error exists in the condition statements or 
in the Q-variables then the FORTRAN compiler will detect it only if 
it is a syntactical error. However, errors such as incorrect or 
non-existent names of in-line FORTRAN functions are not detectable 
by the FORTRAN compiler. These will cause execution errors. If 
such occurs they are caused by errors in the input data. They can 
be detected by checking the cross reference listing in the 
compilation listing of Subroutine UDAC and by checking the load map 
for unsatisfied externals. There should be no undefined variables 
in UDAC and the loaded LATDYN program should contain’ no unsatisfied 
externals for either functions or subroutines. If these errors 
exist, then the user should correct the spelling or other error in 
the input data and rerun the pre-processor. 


LATDYN Calculation Routines 


After successful execution of the pre-processor portion of the 
program, the calculation routines of LATDYN may be executed. The 
instructions for doing this depend on the computer being used. 
Appendices B and C provide these instructions for the CDC/NOS and 
VAX/VMS computers. 

Po st- processor Component 

The purpose of the post-processor component of LATDYN is to 
produce graphical output for nearly all quantities evaluated in the 
calculations portion of the program. When the input data contains 
a PLOTINC command, LATDYN creates a plot file which contains 
calculated results ati even time intervals as specified in the 
PLOTINC command. Section 5.18 describes the PLOTINC command. 

These results may be plotted using the post-processor portion of 
the program. Commands for the post-processor are described in 
Section 5.18. These commands determine the actual plots to be 
created. The pos t- processor may function in an interactive or 
batch mode on the CDC/NOS machines, but only interactively on the 
VAX/VMS machines. 

Post-processor in the Ba tch Mode - On CDC/NOS, the batch post- 
processing mode is activated using the instructions of Appendix B. 
Batch post-processing requires a file POSTIN. When post-processor 
commands are placed in the input data the pre-processor creates 
POSTIN which contains the post-processor commands of the input data 
preceeded by the word "batch' and followed by the word "exit". The 
file may be edited using the system editor and commands added or 
deleted. Similarly, a POSTIN file can be completely created with 
the system editor. 


A-H 


LATDYN Structure 
Pos t- processor 



P o3 t- proc es sor in the Interactive Mode - Instructions for 
activating the post-processor in an interactive mode are supplied 
in Appendices B and C. In the interactive mode, each post- 
processing command is entered directly from the terminal. During 
an interactive session, the user has access to two helpful online 
commands, namely, HELP and INFO. The HELP command allows the user 
to review the available post-processing commands of section 5. 1 8 
with explanation of their functions. The INFO command allows the 
user to review the data associated with the output plot file which 
contains the results of the LATDYN execution. This command is 
useful for it is easy to forget what data is on the output file. 

The INFO command provides the the data case title, the initial 
time, the final time, the number of time steps, the current time 
window, the number of model grid points, the number of windows and 
lists of grid point identification numbers, Q-variable numbers and 
condition numbers. 
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APPENDIX B: EXECUTION ON CDC CYBER 170 SERIES 

COMPUTERS UNDER THE NOS ‘ OPE RAT IN G SYSTEM 


This Appendix presents the procedures for executing LATDYN on 
the- CDC/NOS computer system operating at the NASA Langley Research 
Center. Procedures for other computer installations will be 
similar . 

LATDYN consists of three components: a pre- processor , LATDYN 

calculation routines and a post-processor. The user has a number 
of options available to execute the various components of LATDYN 
using procedure files existing on the LATDYN library (UN=LATDYN) at 
LaRC. These include: 

1) execution of pre-processor and LATDYN calculation 
routines 

2) execution of post-processor 

3) execution of all three components 

For any of these options, execution is carried out via 
procedure files. At LaRC these files are obtained from UN=LATDYN 
using the following, 

GET,"PR0C FILE N A ME " / UN = LATD YN 

Also note that most of the procedure files require the user to have 
a file named "ZLATDYN" on his user number which contains the 
following information, 

USER, charge number , password . 

CHARGE, xxxxxx,LRC. 

DELIVER delivery Information 


1_) E x ecution of Pre-Processor and LATDYN Calculation Routinesj^ 

a) The user can execute LATDYN in a batch mode, an 
interactive mode or a combination of each. The pre- processor may 
be executed interactively to enable the user to carry out a 
syntactical check on the data set. To execute the pre-processor 
interactively and the calculation routines in either batch or 
interactive modes the procedure files "LATRUNB" or "LATRUNI" are 
used. The commands take the form, 

LATRUNB, 

or DATA= , PL0T= ,MACH= ,0UT= 

LATRUNI, 

where. 
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DATA=name of the input data file 

PLOT = name of the output data file for post-'processing 

( default = TAPE1 2) 

MACH = c ompu t er ■ 3 pec i f icat i on for batch execution ( def ault = RHC ) 

OUT=narae of output file for printed results ( def aul t=LATOUT ) 

Following t he , s uc c es sf ul exec ut i on of the pre-processor in an 
interactive mode, the LATDYN calculation routines are executed in 
the batch mode if "LATRUNB" was used or the interactive mode if 
"LATRUNI" was used. The output data file for post-processing 
produced by either a "LATRUNB" or "LATRUNI" procedure file will be 
in indirect access form. 

b) A batch job can also be submitted from the terminal 
without executing the pr e-processbr interactively. Thus all 
syntactical checks are performed in the batch mode, and the user 
will not immediately know if the Job is able to be executed. 

In this case, the procedure files "LATBTCH" or "LATDIR" are 
used. "LATBTCH" is used if the output data file for post- 
processing is to be indirect access and "LATDIR" if the output data 
file for pos t- process! ng is to be direct access. The command takes 
the form, 

LATBTCH, 

or DATA=^ , PL0T = ,MACH = ,OUT= 

LATDIR, 

2) Execution of Post-processor : 


a) Interactive pos t - process i ng can be done at any time 
following the execution of LATDYN calculation routines provided the 
output data file for post-processing has been saved. Following the 
interactive session the user has the option to get hard copies of 
all the plots made by following the instructions which appear on 
the screen. 

For interactive pos t- process! ng , use procedure file "LATPLOT" 
when output data has been saved in indirect access form and use 
"LATPLTD" when output data has been saved in direct access form. 
Commands take the form, 

LATPLOT , PLOT= 

or 

LATPLTD , PLOT= 

Following the interactive plotting session the user will 
automatically receive instructipns for obtaining hard copies of the 
plpts via the procedure file LATPOSt. ilATPOST contains the 
commands for converting the plot file, created by LATPLOT or 
LATPLTD, into a SAVPLT file necessary for the creation of hard 
copies. The command for obtaining hard copies is. 
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LATPOST, HPLOT 


, DEVICE= 


where, 

HPLOT=name of the savplt file to store the hard copy plotting 
information ( def aul t=SAVPLT ) 

DEVICE=type of hard copy device to create hard copies 
(currently limited to VARIAN) 

b) The user can also perform pos t- process i ng in the batch 
mode by using the procedure file "PLTBTCH". This file requires, in 
addition to the output data file for post-processing, that the user 
have a command file which contains the plotting instructions for 
executing the post-processor. The form of the command is, 

PLTBTCH, PLOT DAT = , DATA = , DEV = ,TYPE = , , 

MACH= ,PLTSV= 


where , 

PLOTDAT=name of the LATDYN output data file for post- 
processing ( def aul t=»TAPE1 2 ) 

DATA=name of the file containing the plotting instructions 
(default=POSTIN) 

DEV=type of hard copy device to be used and can take the 

form, 

VAR (Varian plotter), (default) 

VARF (Varian fanfold) 

CALI 1 (Calcomp) 

CALCM (Calcomp 33 inch flatbed) 

MOVIE (Movies) 

TYPE=output plot file specification 

TYPE=IND if the file "PLOTDAT" is indirect access 
(default) 

TYPE=DIR if the file "PLOTDAT" is direct access 
MACH=computer specification ( def ault=RHD ) 

PLTSV=name of the savplt file which the hard copy plotting 
information is to be stored under (default=SVPLT) . 

3) Execution of All Three Components: 


The user can also submit a single batch job which will 
execute the pr e-processor , the LATDYN calculation routines and the 
pos t- processor . The plots generated by the execution of the post- 
processor are in accordance with the plotting instructions included 
in the LATDYN data file as described in Section 5.19. 

j To execute in this fashion use procedure file LATDYNA in the 
following manner, 

LATDYNA ,DATA = , PLOT = , POST DAT = ,TYPE = , 

DEV= ~,MACH= ,0UT= ,PLTSV= 
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where , 


DATA = name of input, data file 

PLOT=narae of output data file for post-processing 
( def aul t=TAPE1 2) 

POSTDAT = name of file c.orita i hi ng- the plotting instructions 
(defaults POSTIN) 

TYPE, DEV, mACH; OUT and' PLTSV are as defined above 
Note: This execution option is hot currehtly avai-lable. 
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APPENDIX C: EXECUTION ON DEC VAX COMPUTER UNDER THE 

VAX/VMS OPERATING SYSTEM 

The VAX computer is a 32 bit virtual memory machine and is 
well able to handle the computational requirements for LATDYN. 
Though LATDYN was originally written for the 60-bit CYBER 6000 
series of computers, it can produce good results using the 32-bit 
arithmetic of VAX's single precision mode. However, only the 
double precision version of the program should be used if 
independent verification of results are unavailable. 

There are two files for executing LATDYN on the VAX system. 

The files are 'LATINT.COM' for the interactive execution of the 
program and 'LATSUB.COM' for the batch execution of the program. 
Because of the long execution time required for the double 
precision version of LATDYN, it is recommended that the user use 
the batch mode unless the 'run is very short. Both files are used 
interactively and require user response to queries. 

The VAX VMS operating system uses the § symbol to initiate the 
execution of a command file. File LATINT.COM and LATSUB.COM are 
executed by typing '@LATINT' or '§LATSUB'. The '.com' is not 
included in the command. 

Interactive Execution 

To execute LATDYN interactively on the VAX/VMS system, the 
instruction "§LATINT" is used. When entered, the program will 
request the names of the data files to be used and created. The 
files to be named are the pre- processor echo file, restart input 
data file, restart output data file, plot output. data file and 
print output data file. The pre-processor echo file is not 
supplied by the user, but is used only to echo the input data file 
onto the print output data file. The restart input file is 
required for a restart and is generated by a previous LATDYN 
execution. The default name is F0R010.DAT. The restart output 
data file must be saved if the user wants the to perform a restart 
in the future. The default name is FOR010.DAT. The plot output 
file is created if the PLOTINC command is present in the input 
data. The default name is F0R012.DAT. The print output file name 
must be given if the user does not wish to use the default name of 
LAT0UT.DAT. 

Batch Execution 

To execute the calculation routines of LATDYN in the batch 
mode with the pre-processor component of LATDYN executed in the 
interactive mode, the instruction "SLATSUB" is used. After 
successful pre- processing , a submittal job file is automatically 
created and submitted for batch processing. The use of LATSUB is 
recommended for all but very small jobs. 

After entering the instruction "^LATSUB" the program will 
request a job root number. The root number will allow the user to 
submit simultaneous cases. Next, requests will be made for the. 
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names of files to be used or created just as is don§ in the case of 
interactive execution. 

If the pre-processor executes successfully, the submitted job 
will be placed in the system queue for execution. The user is able 
to check the status of the job by using the VAX DCL command "SHOW". 
The submittal command uses the VAX DCL command "NOTIFY" to inform 
the user when the program’ has" completed execution'. The user is 
able to logoff the system or’ perform another task if so desired. 
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REFERENCE: LATDYN THEORY 


Following is the reference document which describes the 
theoretical basis for the LATDYN computer program. 
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CONVECTEO TRANSIENT ANALYSIS FOR LARGE 
SPACE STRUCTURES MANEUVER AND DEPLOYMENT 


Oerrold Housner.* 

NASA Langley Research Center 
Hampton Virginia 2366S 


Abstract 


An application of convected finite element 
transient analysis to maneuver and deployment of 
flexible multi-member trusses and frames for appli- 
cation to large space structures is presented. The 
use of convected analysis within the finite element 
framework permits accurate treatment of both large 
elastic deformations and large rigid body rotations, 
while permiting the utilization of a vast storehouse 
of existing computional finite element technology. 

A consistent mass (rather than a lumped mass) formu- 
lation of the equations of motion is developed so 
that exact rotational inertial properties are re- 
tained in the analysis. This is critical for 
space structures which undergo unlimited rotations. 

A shortcoming of consistent mass formulations for 
rotating elastic systems is that they lead to iner- 
tial tenns which depend upon elastic motion and thus 
require up-dating. It is shown that this shortcom- 
ing may be avoided by making appropriate approxima- 
tions with little loss in accuracy. The analysis 
is used to study maneuver of slender booms, deploy- 
ment of axisymmetric hoops composed of flexible 
members and deployment of planar unfolding multi- 
flexible-meinber structures. It is found that 
multi-member structures can possess natural deploy- 
ment patterns such as a sequential (member-by-mem- 
ber) pattern. Such patterns may be highly desir- 
able features of deployment design and their 
achievement without the use of complex control 
mechanisms could represent considerable cost 
savings. 


Nomenclature 


(F) 

(G) 






[Kc] 




rs 


vector of external forces 

vector of nonlinear inertial terms 
for the i^^ finite element. 

total vector of nonlinear inertial 
terms 

components of (G^); (g!|)= 

j=l 

entries of (6^^^)) as given in 
Appendix B; m£ 6; n^ 6 

bending displacement of i^^ finite 
element beyond the i^^ convected 
coordinate, eq, (13) 

integrals of the flexural shape 
function; U i£ N; 5, see 

Appendix B 

rotational spring constant at joints 
of deployable structures 

spring constant matrix, see eq. (30) 

length of i^^ finite element 

mass of i^^ finite element 

entries of l£ r£ N ; 

s£ N, see Appendix A 




coefficients of i^^ beam element 
flexural shape function 


^i 11 • ^i,12 integrals of flexural shape func- 
^i*21 * ^1*22 given in Appendix B 


[B] 


connectivity matrix, eq. (20) 


[M] 

[m;] 


[C] 

[D] 


(di) 


(d) 



(f) 


constraint matrix, eq. (22) 
damping matrix 

vector of displacement degrees- 
of-freedom for i^^ finite elenient 

full vector of displacement degrees- 
of-freedom, see eq. (20) 

vector of generalized forces for i^^ 
finite element, see Appendix A 

full vector of generalized forces 


N 


r, R 


Rh 

T 


assembled mass matrix 

mass matrix for i'-^ finite element 

components of mass matrix; 

3 

[M^] = £ 

j = l 

total number of members 

total number of independent degrees 
of freedom 

local radius of hoop member dur- 
ing deployment at its center and 
either end, respectively 

fully deployed hoop radius 

local dimensionless coordinate 
along convected beam finite 
element; 0 < s < 1 


* Aerospace Engineer, Member AIAA 
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[Ti] 

convected coordinate transformation 
matrix for i^*^ finite element 

U, V 

local displacement of beam element 
parallel to x and y axis respec- 
tively 

X, y 

global cartesian coordinate system 

“i 

angle formed by 1^^ undeformed 
finite element and x axis 

'’i.ll • '‘i.12 . 
•’1,21 » •’1,22 

integrals of flexural shape 
functions, see Appendix A 

6 

variational operator 

%,i 

mean axial strain of straight 
line joining end points of 
ll^ finite element 

^o,i 

mean axial strain of i^^ deformed 
finite element 

■•i ,1 •''i ,2 

end point displacements of i^*^ 
finite element parallel to x axis 

6i 

angle formed by straight line 
joining the end points of i^*^ 
finite element and the x axis 

^1.1 ’ ^1.2 

end point displacements of i^^ 
finite element parallel to y axis 

*1,1 • *1,2. 

rotations of i^*^ finite element 
end points 

'C 

angle formed by hoop member and 
vertical axis, see figure 7 


Introduction 

Recently, interest has grown in establishing 
permanently manned space stations. Space stations 
may be composed of multi -body interconnected struc- 
tural components some of which are thin shells 
(habitation and work modules), panels (solar arrays 
and radiators) or beams (space booms, towers and 
robotic manipulator arms). These cor^ponents or 
the entire space station may, from time to time, 
be maneuvered or, as in the case of robotic arms, 
manipulated through large rotations. Such maneuvers 
result in dynamic loads on the various space station 
components. It is likely that to obtain structural 
component loads for preliminary design, maneuvers 
of large beam-frame models will be analyzed. 

In addition, many station components may be 
sent aloft packaged and then deployed in orbit. 
Deployment concepts often involve large rotations 
of flexibl e truss-type members. Since deployment 
experiments are both difficult and expensive to per- 
form either on the ground or in orbit, dynamic 
analyses will be relied upon to validate performance 
of proposed deployment concepts. 


In any of the three application areas, maneuver, 
manipulation or deployment, dynamics of intercon- 
nected flexible components undergoing large rota- 
tions is involved. . In some cases involving flexi- 


bility and large angle rotation, flexible and 
rotational responses are uncoupled and a linear 
solution is possible. 

This situation can occur when the following condi- 
tions are sati sfied: 

(i) the structural modes (including rigid body 
modes) do not change with* time 

(ii) the contributions of the externaj forces in 
each modal shape (including rigid shape) do not 
depend on the structure's rotated position 

(iii) the structural deformations are small 

(iv) the nonlinear centrifugal and Coriolis forces 
are negl igibl e. 

When only some of these conditions are met, it 
may still be accurate to perform a quasi -linear 
analysis in which the modes and/or external force 
contributions in the modal shapes are updated peri- 
odically in the course of the analysis. However, 
there exists a large class of problems in which 
the nonlinear coupling between flexibility and 
large angle rotation must be accounted for. 

In space structures problens requiring a non- 
linear analysis, the strains in the structural 
components are small even though rotations and 
hence displacements are large. To avoid the compu- 
tations associated with accurate strain measures 
for large motions, it is desirable to separate the 
rigid body rotations from the structural deforma- 
tions. This separation may be accomplished by 
using an updated or convected coordinate system. 

Over the past decade, convected coordinates 
have been succesful 1y used in conjunction with 
finite element models in performing transient res- 
ponse analyses using explicit time integration 
algorithms. (For example, see references 1-4). 

Thus, convected transient analysis techniques offer 
a viable approach for analyzing maneuver, manipula- 
tion and deployment. Two additional benefits to 
this approach are the utilization of a vast store- 
house of finite element computational technology and 
the capability of the convected analysis within the 
finite element framework to account for both large 
elastic deformations and large rigid body rotations. 

Even though convected transient analysis has 
seen considerable use, its application to fundament- 
al configurations which are likely to arise in 
space stations and their components has not receiv- 
ed much attention. Previous developments in con- 
vected frame analysis have been limited to struc- 
tures containing large concentrated masses, such 
as automobiles where the engine, transmission, 
differential, etc, may be treated as lumped masses. 

A lumped mass formulation has considerable numeri- 
cal advantages over a consistent mass formulation. 
However, many flexible components of large space 
systems do not lend themselves to an accurate 
lumped mass treatment. (See reference 5.) For 
example, large space booms and antenna structures 
have uhifonp or slowly varying distributed mass. 
Lumping*^ the d1 stributed mass leads to overestimat- 
ing the rotational inertia and thus underestimating 
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its rigid body rotational rate. In addition, the 
coupling between rigid body motion and structural 
deformation due to nonlinear centrifugal and 
Coriolis forces or large structural deformations 
may be inaccurately represented in a lumped mass 
formulation. Thus a consistent mass formulation 
is desirable. 

The purpose of this paper is to present the 
development and application of a consistent mass 
finite element formulation for convected transient 
analysis of two-dimensional beam-frame configura- 
tions for use in space. The computational conse- 
quences of a consistent mass formulation are dis- 
cussed with emphasis on reducing computations 
through appropriate approximations. The analysis 
is applied to rotational maneuver and translation 
of long slender booms, deployment of an unfolding 
multi-flexiblemember axi symmetric hoop and deploy- 
ment of a planar unfolding multi -flexible-member 
structure. 



Figure 2. Finite Element Deformation and Convec- 
ted Coordinate System 

(a) Rotation of Convected Axes 

(b) One, Two and Four Finite Element Modeling 
of a Single Beam Member 

Coordinate system 


Convected Analysis 

The purpose in using a convected analysis is 
to separate rigid body rotations from structural 
deformations for in each finite element. The need 
for this separation arises since any truncated 
form of the physical strain measure in a fixed coor- 
dinate system is inaccurate for large rotations. 
Separation may be accomplished by setting up a 
collection of convected coordinate systems which 
follow the rotation of each finite element. The 
convected analysis used herein deviates from that 
of references 1 thru 3 in that the development of 
the inertial terms uses a consistent mass formula- 
tion rather than a lumped mass formulation. The 
consistent formulation is important when rigid 
body rotations must be accurately accounted for. 
To illustrate this consider the rotation of a 
single rigid beam member. Figure 1 illustrates 
the error from the exact rotational Inertia asso- 
ciated with a lumped mass or diagonal mass matrix 
formulation. This simple Illustration demonstrates 
that when a lumped mass formulation is used, five 
rigid finite elements are required to approximate 
the exact rotational Inertia to within eight per- 
cent. On the other hand, a consistent mass formu- 
lation requires only one rigid finite element for 
the exact value. 



Consider the deformed beam finite element 
shown in figure 2a. The convected motion of the 
element is defined by the straight line which 
joins the ends of the element. Deviation from 
this lines constitutes the deformation of the ele- 
ment. A convected coordinate system, local to 
each element is thus defined which rotates with 
the element, but does not translate with it. (The 
convected coordinate system could be defined to 
translate with the element, but this is not neces- 
sary.) The convected system in general differs for 
each element. Coordinates in the global stationary 
coordinate system are related to those in the i^” 
convected system by 



where barred and unbarred quantities denote mea- 
surements in the convected and stationary coordi- 
nate systems respectively and 


E osS^ - sinSi 
1n6^ cos^i 

The convection angle 6^ for the 
defined by 


( 2 ) 

i^^ element is 


tanei = (Li,y + ACi)/(Li,x ♦ (3) 

where L^x and Liy are the projected lengths of the 
ith undeformed element on the x and y axes respec- 
tively and 

Ani = ni^2 - Hi j 

^^i = ^1,1 - 

in which n, i , ni^2* ^1,1 ^i 2 * are the values 

of grid point displacements at’ the ends of the 
i^" element and parallel to the x and y axes as 
shown in figure 2a. 


The displacements in the convected system at any 
point on the 1^^ element transform as 


Figure 1. Percent Error in Centroidal Rotational 
Inertia For Lumped Mass Modeling of Beam Member 
[Error=(Approx. value)/(Exact value) - ij 



(4) 
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However, the velocities transform as 


Virtual Work 


ORIGINAL PAGE fs 
POOR QUALITY 


f I • 

> [T,-]U U [Tj][i*]L 
V / \v Vv 


where. 


[1*3 


p -1 

1 0 


(5a) 


and the accelerations transform as. 







In equation (5b) the second and fourth terms on the 
right hand side are recognized as the Coriolis and 
centrifugal contributions, respectively , to the 
absolute acceleration vector. Converse relation- 
ships for equations (5a) and (5b) may readily be 
derived. ' * ' ’ 


Rotation and Strain Rates 


The virtual work contribution from the inter- 
nal forces of the i^*^ element Is given by, 

I (10) 

where Pj is the average axial load over the length 
of the 1t;h element and M is the bending moment in 
the i^^ eluent given byi 

Mj = -(EI)^ 

From eqs. (7) and (9) , 

®^i (11) 

Substituting eqs. (11) and (8b) into eq. (10) gives, 

■ 1 

6u^=p.(aAn^cose^+6AC^sin6^)+(l/Lf) /M^ 6h^ds 
1 0 

+ (1/Li)/Pih{6h{d? (12) 

0 


The rotational rate of the straight line join- 
ing the end points of the i^^ finite element may 
be found by differentiating eq. (3) with respect 
to time and the mean axial strain rate may be 
found by differentiating with respect to time the 


rel dtion 

^i ,m 

,x 


The result 

is 


/ 

) 

= (l/Lj)[T^]7 . 


)0i/ 



For small strains, on the left hand side of eq. 
(7) may be neglected compared to unity. 


Herein the flexural displacement is defined as, 

-Cai,i(s)+ai^2(s)](0i-“i) (13) 

where 4>i,l and ^^^2 rotations of nodes 1 and 2 
at the ends of the i^^ element. Substituting into 
eq. (12) yields, 

«U, = (fi)(«di) (U) 

where 

= (’^i ,1 ,1 *^i J *’^1 ,2 ^i,2»*i,2^^ 

and the expressions for entries of the internal 
force vector (f-j) are given In Appendix A. 


Finite Element Displacements 

The finite element displacements may be approx- 
imated in the convected system as, 

ui(s) = Hi j + sAni (8a) 

Vi(s) = €j,i + sA5j + hj(s) (8b) 

where, h^ is the flexural displacement and 

hi(0) = hi(l) = 0 
and 

0 <1 <) 

The local strain at the neutral axis of the i^^ ele- 
ment may be found from by adding on to p, the 
usual first-order neutral axis strain due to flexure 
to yield, 

S.o = " (’/2L?)(h!)2 (9) 

and the axial strain at any point through the 
thickness of the element is, 

= s.i - (z/Li)Vi" 

where z is a thickness coordinate measured from the 
neutral axis. 


For a consistent mass formulation, the contri- 
bution of the 1 th element to the total virtual 
work done by the inertial forces is expressed as, 

«w, = m^^’(«u^, 6v^)(V^,V^)Tds (15) 

where m^ is the uniform mass of the 1^^ element and 
the expressions for 6 and 6 v^ may be found by 
substituting the variations of eqs. (8a) and (8b) . 
into eq. (5) to yield 



+ 

+ 



66^cos0i-6h^s1neA 

6B^sin0i+ah^cos0i/(16d) 


The expressions for Uj and v^ are found by differ- 
entiating eqs. (8a) and (8b) and substituting 
into eq.(6) to yield, 



h^ si nS^ -20^ h^ cosO^^-h^ 5?s1 n0^ - 0 ^ h^ cosbA 
^ h|Cos0^ -20^ h^ si n0^ -h^ 0?cos0^ -0^ si n0^ / (16b) 
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In each of equations (16a) and (16b) the first 
column vector on the right hand side arises from 
the rigid body motion of the finite element while 
the second term arises from the flexural deforma- 
tions of the element, Coriolis and centrifugal 
terms are explicit in the the second term of eq. 
(16b) and implicit in the first term. The implicit 
effect may be observed by transforming the first 
term of eq. (16b) to convected coordinates vi a 
eq. ( 6 ) and then evaluating the result at the ele- 
ment end points. 

When the product of eqs. (16a) and (16b) is 
taken to produce the integrand of eq. (18) and the 
indicated integration is performed, rigid body ro- 
tation and flexural deformation terms are coupled. 
The result for the i^^ element may be expressed as 

6W^ = m^(6d* + m^(6di)T(G^) (17) 

where [M^] is the mass matrix of the i^*’ finite 
element, and (Gi) is a vector of nonlinear kine- 
matic terms. 

The mass matrix contains terms due to rigid 
body motion and the predominant inertial coupling 
between rigid body motion and flexural deforma- 
tion. The treatment of the mass matri x i s sim- 
plified by observing that it may be separated into 
three distinct parts, namely, 

3 

[M^] = 

j = l 

The entries yiven in Appendix 

B. Entries of [M^ ' * ' ] are i ndependent of i and 
arise from rigid body motions and stretching of 
tbe.K finite element while those of and 

arise from flexural motions Jha.i^ fi- 
nite element. The entries of both and 

depend on the convected rotation and 
hence would generally require up-dating during the 
analysis. Furthermore, those entries requiring up- 
dating contain h^j . 

As the size of the finite element is reduced, 
the value of h^ decreases as^iSv noted in figure 
2 b and then the entries of also decrease. 

Thus for a sufficiently ,ttne finite element grid, 
the contribution of to the total virtual 

work becomes negligible and if retained would re- 
quire few up-dates during the analysis. Further- 
more, if an element is rigid, then from eq. (13), 

♦i ,1 = = ®i (18) 

It may^then shown that the contributions of 
( 6 d^ )[Mi ](d^ ) to the virtual work of the 
inertial forces vanishes for a rigid element and 
furthermore for a stiff element, approaches zero 
as h approaches zero. It follows, that lor^Cela- 
tively stiff elements, few up-dates of will 

be required. For the cases considered herein, it 
was found that even for booms with slenderness 
ratios up to 1000 . still required few 

up-dates and could be entirely neglected 

when two finite elements were used per half-wave 
of beam member deformation. 

Similarly the vector (Gi) may be separated as, 

j=l 


whose entries are given in Appendix B, These vec- 
tors arise from coupling between element flexure 
and nonlinear centrifugal and Coriolis effects. 
All five vectors vanish as h^ approaches zero. 
Furthermore, for small hj , and 

(G^v^l) behave as h^, while (G^^^O behaves as 
h?. (This does not imply that all centrifugal 
and Coriolis effects vanish as h-j approaches 
zero, for these effects are implicit in the piece- 
wise representation of the entire member). 

As a consequence of the foregoing discussion, 
it is clear that many terms in the governing equa- 
tions will be quite small in many applications and 
either be negligible or require few up-dates. 


Equations of Motion and Constraint 

The equations of motion may be derived from 
the principle of virtual work using eqs. (14) and 
(17) 

Z (f,)} =«We ( 19 ) 

i=l 

where 6 We is the virtual work of applied external 
forces and N is the total number of beam members. 
Lumped masses are accomodated by appropriate addi- 
tions to the diagonal entries of [ni]. Equation 
(19) yields the equations of motion if the entries 
of ( 6 oi) are independent. However, they are not 
independent since members are connected together 
and since multi-point constraints or other cons- 
traints may exist. An independent set of displace- 
ments can be formed by utilizing the constraints 
on the system. 


Connectivity Constraints 

Connectivity relations may be expressed by 
the matrix equation, 

((^i)(? 2 )j = CB](d) (20) 

where (d) is a vector of Nq degrees of freedom. 
The connectivity matrix [B] will usually contain 
many zeros and ones. Rather than substitute eq. 
(20) into eq. (19), connectivity is usually account- 
ed for in finite element computer implementation 
when the mass matrix and internal force vectors 
are assembled. This is accomolished by appropriate 
superposition of entries in (G-j), (r^) and [frj] as 
directed by the connectivity array and computer 
storage addresses. Thus after assembly, eq. (19) 
effectively takes on the form, 

( 6 d)T{[M](d) + (G) +{f)} = ( 6 d)T(F) (21) 

where (F) is a vector of external forces. 


Linear Multi-point Constraints 

For N(. (less than Np) multi-point constraints 
one has. 


t^DD ! ^ 
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where the constraint matrix and displacement 
vector (d) have been partitioned Into dependent 
and independent sets. Ihe selected set of depen- 
dent degrees-of-freedom be expressed In terms of 
the independent set of degrees-of-freedom as 

(do) = [C](di) 


where. 


[c] = -[Coo]-’[Coi] (22) 

Eq. (22) could be substituted into a partitioned 
form of eqs. (21). However, in some applications, 
the multi-point constraints can change many times 
during the course of analysis. Tfiis. can be due to 
lock-up of certain members in a deployment analysis, 
contact conditions as in assembly, docking and ber- 
thing or articulation of structural components. 
In such cases it may be more computationally effi- 
cient to eliminate the j th dependent degree-of- 
freedom by multiplying the j th row of [M], (G), 
(f) and (F) by Cj)^, (where cj^ is an entry of the 
constraint matrix}, and adding the result to all 
the other rows of [M], (G), (f) and (F) respective- 
ly, for all k and j in the ranges 1< k< Nq and l£ 
j_< Nc* This effectively eliminates A simi- 
lar procedure is followed for the columns of [M] 
thereby effectively eliminating (dp). 


Constraints Due to Rigid Members 
Axial rigidity 

The iTb member is defined as axially rigid if 
the neutral axis strain, its variation and its 
derivatives vanish. Then from eqs (11), 

6nij = fiACitan 0 i (23a) 

^i ,1 “ ^^i,2 (23b) 

where, for simplicity, the first-order nonlinear 
term of eq. (11) caused by local finite element 
flexure has been dropped. Differentiation with 
respect to time of eq. (23b) yields, 

n- 2 + AC^tane^ + AC^e^sec^e^ (23c) 

and from eqs. (7) and (23b) 

§1 = ACi/(LiCOsei) 


ROOT 

ROTATION 

(RAOIANSI 



Figure 3. Hub Rotation Due To Applied Torque Step 
Load 


Equations (23a) and (23c) ma[y be used to eliminate 
and m i in eq., (21) . by appropriate mani- 
pulation of the rows' of . (G) , (f) and (F), and the 
rows and columns of [M]. The nonlinear term in 
eq. (23c) cannot be dropped as it is not a negli- 
gible term. It may be handled through appropriate 
augmentation of the Internal force vector (f). 


Flexural rigidity 

If the 1^^ member is flexural ly rigid then,. from 
eq. (18), 

64 ) 1,1 = 

4i,l = ^,2 

and from the second of eqs. (7), 

^1 (24a) 

.l=5^,2-^’ii^®^®1“'*-lil,2/cos0i (24b) 

Differentiation of eq. (24b) with respect to time 
yields, 

,1 =5i ,2-Anitan6i -Li?i ,2/cos6i -ei ^„,eiLiCOse^ ^ 

Eqs. (24a) and (24c) may be used to eliminate 
*♦ 1 , 1 . ii,l. «5i,l. ^i.l eq. (21) by appropriate 
manipulation of the rows of (G), (f), and (F) 
and the rows and columns of [M]. The nonlinear 
term in eqs. (24c) cannot be dropped as it is not 
negligible. It may be handled through appropriate 
augmentation of the internal force vector (f). 


Numerical Solution of Equations of Motion 

Equation (24) may be solved numerically by 
selecting from any of a wide choice of time inte- 
gration algorithms. (See for example the survey 
article of reference 6.) An explicit algorithm is 
chosen here because it can readily accommodate a 
host of nonlinearities, damping mechanisms and con- 
trol actuators. In particular a modified central 
difference algorithm is selected. This leads to 
the following recursive procedure, 

(d(t))= -tCM] 4- (l/2)At[D]}-M(f(a(t))) 

+ (G(d(t) ,dp(t) ,dp(t) )+[D](d(t-4t) )-(F) } 

(d(t))=(d(t-4t))+(l/2)At(d(t)+d(t-At)) 

(d(t)) = (d(t-At)) + At(d‘(t)) (25) 

where (dp) and (dp) are velocity and acceleration 
predictors which are used only in conjunction with 
terms on the order of the flexure of a finite 
element. The predictors are chosen herein as, 

.* • .. 

dp(t) = d(t-At)+Atd(t.At) (26) 

dp(t) = d(t-At) (27) 

and [D] Is a linear viscous damping matrix. 
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Applications 


The convected transient analysis presented 
herein is applicable to multi-body flexible truss 
and frame problems which may involve maneuver 
through large angles, translation through space or 
deployment. In this section the following problems 
are considered: slewing maneuver of a long slender 
boom, translation of two slender beams hinged to- 
gether, deployment of an axisymmetric multi-flexi- 
ble-member hoop and deployment of a planar multi- 
flexible member unfolding structure. In each 
application cubic flexural shape functions are 
chosen. 


1.5 

AVERAGE AXIAI 
BOOM FORCE. , 
NEWTONS 


P TORQUEATOTAI ROTATIONAL INERTIA) = 0.0025 
^HUB goOM 

SLENDERNESS = 

50 m ' 50 

TIP MASS/BOOM MASS« 10 
~ HUB ROTATIONAL INERTIA 
TOTAL ROTATIONAL INERTIA 



^CENTRIFUGAL 
FORCE BASED 
ON RIGID 
BODY MOTION 


Figure 5, Axial Load Due To Applied Torque Step 
Load 


Slewing Maneuver 

Figure 3 displays the time varying response 
of a slender boom with a rigid central hub subject 
to an applied torque. The hub rotation for a 
connected flexible boom with slenderness ratio 1000 
is compared with that of a connected rigid boom. 
The hub rotation with a connected flexible boom is 
always equal to or greater than that of a connected 
rigid boom. This occurs because the inertial mass 
accelerated by the applied torque is the sum of the 
hub inertial mass plus the inertial mass of the 
boom from the hub out to the flexural wave front 
and the latter mass is less than or equal to the 
total boom inertial mass. Hence, the variation of 
the hub motion with connected rigid and flexible 
booms depends on the relative rotational inertial 
masses of the hub and boom as well as boom slender- 
ness ratio. For the case displayed in figure 3, 
the inertial mass of the hub is considerably greater 
than that of the boom. Hence the deviation between 
the rigid and flexible responses is small. 

As shown in figure 4, the tip motions of the 
flexible and rigid booms are quite different from 
that of the hub. Sped fi cal ly , the tip motion 
with the flexible boom lags that with the rigid 
boom. This may be understood by considering the 
time delay caused by flexural wave propagation 
along the boom. As might be expected, the diff- 
erence between rigid and flexible boom tip motion 
i s relatively large compared to the di fferences 
at the hub. Thus boom flexiblity can greatly 
influence pointing without significantly influenc- 
ing hub rotation. Finally, in figure 5, the axial 
load experienced by the flexible and rigid booms 
is compared. Axial loads are due to centrifugal 



Figure 4, Boom Tip Displacement Due To Applied 
Torque Step Load 


force. The flexible boom axial load is seen to 
oscillate about the rigid value with increasing 
amplitude and frequency as the boom gains angular 
velocity. 


Translation of two hinged flexible beams by a 
^n^ant-di recti on step load 


Figures 6a and 6b display the calculated mo- 
tion of two identical beams hinged together at 
their ends and subject to constant-direction step 
load at one end. The problem is kinematically 
nonlinear because the torque caused by the con- 
stant-directional load and the load components in 
the beam modes vary with the kinematics. Displace- 
ments are shown in figure 6a for a flexible beam 
and the trajectory of the load application point 
is shown in. figure 6b for both a flexible and a 
rigid beam. Examination of the trajectories of 
the load application point indicate that flexure 
has little effect on trajectory during the early 
phase of motion, but becomes more important as the 
motion proceeds. 


5r 
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Figure 6. Large Distortion and Motion of Two 

Pin-Connected Beams Subject to a Vertical 
Tip Step Load 

(a) Member Deformations 

(b) Load Point Trajectory 


Hoop Deployment 


The hoop investigated here is similar in 
some respects to that of the Hoop/Column Antenna 
concept of reference 7, In the present study 
the complex control linkage between adjacent mem- 
bers of the hoop in reference 7 is replaced by 




Figure 7. Hoop Deployment (N=16) 

(a) Deployment Sequence 

(b) Top View of Two Adjacent Members 

(c) Front View of Two Adjacent Members 

simple linear rotational springs at each joint 
which drive the deployment as shown in figures 7a 
and b. The hoop is assumed to be composed of N, 
(an even integer), identical flexible beam members 
of length L. The members are assumed to be con- 
nected through universal joints which do not trans- 
mit any moments. At any time during deployment, 
the members project on the horizontal plane as an 
N sided polygon with the joints between adjacent 
members lying on a circle of radius R. As deploy- 
ment proceeds, R increases to a fully deployed 
hoop values of R^. 

Due to the symmetry of the hoop, only one mem- 
ber of the hoop, say the member, needs to be 
analyzed provided the effects of the other members 
on it are accounted for. The effect of the other 
members on the member can be systematically 
derived by establishing the virtual work of all 
the members in terms of the motions of the i^^ 
member. Effectively, the other members constrain 
the motion of the i^^ member and it is shown 
herein that this constraint may be modeled as an 
additional mass matrix superimposed on the mass 
matrix of the i^^ member. 

As shown in figure 7b, the if^ member lies in 
a plane normal to and translating parallel to the 
local radius r. Thus its virtual work has a 
contribution from the motion of the member in the 
plane and a contribution from the motion of the 
plane, namely, 


The first term is associVted with member motion in 
the plane and has already been established in eq. 
(n)» It remains to establish the second term. 

The second term is derived from the geometric 
relationship between the translations of the plane 
and motions of the i^*^ member in the plane. For 
the purpose of deriving this geometric relationship, 
the members are assumed to behave rigidly. However, 
this assumption is relaxed later in performing the 
dynamic deployment analysis. 

The angle 2y which is formed by the radii of 
two adjacent joints depends on the number of hoop 
members as 

Y = n/N 

The value of a remains constant during deployment, 
but as R increases during the deployment process, 
the angle ^ formed by the member and the vertical 
axis varies since, 

Rsiny = (L/2)sintJ> (27) 

When fully deployed \|> is n/2 and 

r 

Rh = L/(2sinY) (28) 

Since, 

sini|) = (tib-na)/L 

where Og and nb are horizontal displacements of 
the hoop member ends in a plane normal to r, eq. 

(27) yields, 

RsinY = (l/2)(nb-na) (29) 

The contribution to the virtual work due to radial 
motion is, 

= mr6r 

Noting that 

r=RcosY 

and substituting from eq. (29) gives 

= (m/4)(ntj-n3)(6n^-6n3)cot^Y 

The principle of virtual work applied to the entire 
hoop may be expressed as, 

N6W^n)+ N(«d)''’[M^]{(l) + 2N(6d)%3](d)=0 (30) 

where, 

{d)=(n*i j ^2»^i ,2* ^ ,2* • • • 

and p is the number of finite elements employed to 
model the i^h hoop member. The nonzero entires 
of [Mf,] are, 

= ("ipJq.q = (n>/4)cot2Y 

K^l.q = (™h)q,l = 
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The nonzero entries of [K5] are 

( ^s) 1 ,3 “ ( )6p,6p ' ^ 

where k is the rotational spring constant at each 
joint. Dividing eq. (30) by N reveals that the 
equation of motion for axisymmetric deployment is 
identical to that for a single hoop member whose 
mass matrix is incremented by [Mh] and which has 
grounded rotational springs at each end with a 
spring constant equal to twice the nominal spring 
constant value. 

As an example application, consider a twenty 
member hoop (N=20). Figure 8 displays the predict- 
ed motion of two adjacent flexible hoop members 
during the deployment process. Each member has a 
slenderness ratio of 1000, Hoop member deformation 
is shown at one quarter, one-half and near fully 
deployed. During deployment, each member deforms 
into a two half-wave pattern which is consistent 

with the symmetry of the hoop- In figure 9, the 

time for the hoop to reach full deployment is 

examined over a range of N values. Two cases are 

displayed. In one case the individual hoop member 
length is fixed so that as N increases, the fully 
deployed hoop radius and total hoop weight also in- 
crease, while in the other case, the fully deployed 
hoop radius is fixed so that as H increases, the 
individual member lengths and mass decrease, however 
the total hoop weight remains essentially unchanged. 

In each case results are shown for hoops con- 
taining all rigid or all flexible members and flex- 
ibility is seen to have little affect on deployment 
time. 


The trends for the two cases are different and, 
where they cross, the hoop geometries are identical. 
The trend for each case may be approximately estab- 
lished in closed form by considering only the 
influence of the diagonal entries of the added 
mass matrix [M^] on a rigid member. Such a consi- 
deration reveals that for the fixed length case 
the deployment time varies linearly with N while 
for the fixed radius case it varies inversely with 
the square root of N. The deployment analysis 
confirms these trends. 


TOP VIEWS 



Figure 8. Top and Front Views of Two Adjacent 

Hoop Members During Deployment of Twenty 
Member Hoop, k=2.82 m.N/radian 



Figure 9. Deployment Times For Hoops Composed of 
Various Numbers of Members 


Multi-member Unfolding Deployment 

An accordian-like flexible beam deployable 
structure having N pin-connected members is shown 
in figure 10. When fully deployed the structure 
is a continuous flexible boom. Thus, when the 
angle between two adjacent members reaches 180 
degrees, the joint is assumed locked and thereafter 
is a rigid joint. The deployment of the structure 
is driven by rotational springs located at each 
joint. All springs are identical, initially com- 
pressed and fully relaxed when the joint they are 
located at locks-up. Considerable insight into 
unfolding deployable structures is gained by exam- 
ination of this generic structure. 





Figure 10. Geometric Configuration of Planar 

Unfolding Structure Showing Rotational 
Joint Springs and Joint Numbering 
Pattern. (Not An Actual Deployment) 


Even Number of Members (N=2,4,6,.,. ) 

N=2 - Figure 11 displays the calculated deploy- 
ment sequence of a two-member structure in a cumu- 
lative superimposed fashion where each member was 
modeled with two finite elements. Member deforma- 
tions have not been scaled-up, but are consistent 
with the drawn member length and little flexure of 
the members prior to lock-up is visible. Following 
lock-up considerable member flexure occurs and the 
deployed boom vibrates in essentially a first sym- 
metric free-free mode shape. 

N=4 - Figure 12 displays the deployment of a 
four member structure in the same fashion as that 
used in figure 11. Prior to lock-up, members 
appear to exhibit more flexure than in the single- 
member case. Since the outer members accelerate 
less mass, the outer joints lock-up first shortly 
followed by the inner joint. Following lock-up. 
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the deployed boom vibrates with essentially a sec- 
ond symmetric free-free mode shape. However^ the 
post-lock-up shape depends upon member stiffness. 
Though not shown, stiffen members can yield a 
first symmetric free-free vibration shape. 

N=6,8,10,...When the structure has six or 
more members, the time difference between the lock- 
up of the joints becomes greater than that exhibit- 
ed in the two-member case. Lock-up begins at the 
outermost joints and proceeds towards the center 
of the structure as shown, for example, in the pro- 
gressive deployment states of the sixteen-member 
structure of figure 13. Also, lock-up times for 
some of the joints in various multi-member configu- 
rations are provided in the Table. 

Joints are numbered beginning at the outermost 
joint and progressing towards the center so that, 
independent of the number of members, the outermost 
joint is always joint number one. Note that the 
lock-up time for joint 1 is independent of the 
number of members and that other joint lock-up 
times also become independent of the number of 
members as the number of members increases. In 
general, lock-up times for all but the Innermost 
joints are independent of the number of members. 
The result is that the structural joints lock-up 
sequentially even though all members and springs 
are identical. Again the reason for this is traced 
to the quantity of mass the various members must 
accelerate. 




Figure 11. Deployment Sequence of Two-Flexible- 
Member Unfolding Structure 



Figure 12. Deployment Sequence of Four-Flexible- 
Member Unfolding Structure 

This basically sequential pattern is consider- 
ed a desirable quality for deployable space struc- 
tures. ( See references 8 and 9.) For the multi- 
member structure considered here, sequential de- 
ployment appears to occur naturally; that is, with- 
out the aid of external mechanisms, controls or 
dampers. It is anticipated that multi-member de- 
ployment configurations other than the one consi- 
dered here will also exhibit natural deployment 
patterns, though they may not be of the sequential 
type. If this is the case, it would allow the 
designer to select member and spring properties to 
achieve a preferred deployment pattern without the 
expensive addition of controlled or passive mechan- 
isms to enforce a prescribed pattern. 

Odd Number of Members (N=3,5,7, , . . ) 

For the unfolding type structures considered 
here the deployment pattern for an odd number of 
members differs from that of an even number of 
members. For example, consider the five member 
structure of figure 14. Due to the odd nunber of 
members, the deployment pattern tends* to appear 
rotated counter-clockwise. Actually, the two 
outer members are rotating clockwise while the 
three inner members are rotating counter-clockwise. 
The result is that the pattern satisfies the con- 
servation of angular momentum. As in the case of 
an even number of members, the deployment pattern 
tends to become sequential as the number of members 
in the structure increases. However, the pattern 
will be distorted somewhat due to the conservation 
of angular momentum. This distortion should de- 
crease as the number of members Increases. 
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Figure 13. Deployment Sequence of Sixteen- 

Flexible-Member Unfolding Structure 

Concluding Remarks 

Tlie application of convected finite element 
transient analysis to maneuver and deployment of 
two-dimensional truss and frame space structures 
has been presented. This approach permits the 
utilization of a vast storehouse of finite element 
computational technology and, within the finite 
element framework, the convected analysis accounts 
for both large rotations and large elastic deforma- 
tions. 


A consistent mass (rather than a lumped mass) 
formulation of the equations of motion has been 
incorporated so that exact rotational inertial 
properties are retained in the analysis which is 
critical for space structures which can undergo 
large rotations. Within the convected analysis 
framework, a consistent mass approach has been 
shown to allow the resulting inertial terms to be 
ordered in relative magnitude and many of the 
terms which require frequent updating may often 
be safely neglected when the finite element grid 
is a relatively fine one. For the problems con- 
sidered herein, little or no updating was required 
even for booms with slenderness ratios of 1000 when 
two finite elements were used per halfwave of 
deformation. Thus, frequent up-dating, a major 
shortcoming of the consistent mass approach in a 
convected coordinate system can be avoided. 

Results of the convected analysis have been - 
presented for maneuvering of slender booms through 
large rotations, translation of two hinged beams, 
deployment of an unfolding multi -flexible-member 
axi symmetric hoop and deployment of a planar un- 
folding multi -flexible-member structure. For a 
symmetrical arrangement of booms extended from a 





Figure 14. Deployment Sequence of Five-Flexible- 
Member Unfolding Structure 


central hub and subject to a torque, flexibility 
influences boom tip deflection and hence pointing, 
but has little influence on hub rotation. 

Deployment of an axisymmetric multi -flexible- 
member hoop may be analyzed by considering only 
one hoop member when an appropriate derived mass 
matrix Is superimposed on the member mass matrix. 
Applying this procedure some aspects of hoop deploy- 
ment have been studied. It has been found that 
hoop deployment time is essentially independent 
of member flexibility and that flexible members 
deform in a two-half-wave pattern. 

Furthermore, planar unfolding deployment of a 
multi-flexible-member structure has indicated that 
these structures may possess natural deployment 
patterns which are very desirable in practice. 

The structure considered herein displays a natural 
sequential deployment pattern without the use of 
any control mechanisms which are usually quite 
complicated and expensive. 


Number 

of 

Members 

1 

2 

Joint 

3 

Number 

4 

5 

4 

1.448 

1.468 




8 

1.453 

2.308 

2.792 

2.814 


12 

1.452 

2.330 

2.856 

3.631 

3.930 

16 

1.454 

2.401 

2.861 

3.680 

4.643 

20 

1.461 

2.476 

2.877 

3.816 

- 


Table. Joint Lock-Up Times In Seconds For An Even 
Number Of Deploying Members 
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